Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC43.1% 103.5 / 13.9 58G384G (GGT→GGGyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (30/3);  new base C (0/25);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.88e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC  >  NC_000913/2984292‑2984585
                                                                                                                                                      |                                                                                                                                                 
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                  >  1:393000/1‑151 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                              >  2:296365/1‑151 (MQ=255)
           gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccg                                                                                                                                      <  1:680446/151‑1 (MQ=255)
            agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                                                                     <  2:611378/151‑1 (MQ=255)
            agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                                                                     >  2:481501/1‑151 (MQ=255)
               aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgccg                                                                                                                                   >  1:692239/1‑150 (MQ=255)
                   aaaaaTGAGGGAGCCAAAAAAATATATGTCGACTTGCCACCCATCATTATTCAACAGATTAATGACGGCAGATATTTACCGGTAAGGAGAAAATGCAGCGCTATTCGTCACTTTAAATGGTAATTGCCACCCCCGGGCCCCCGCCGATa                                                                                                                                <  1:115397/149‑1 (MQ=255)
                                ccaaaaaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATAAATGACGGCAGATTTTTATCGGTAGGGAGAAAAGGCAGAGCTTTTCGTCCCGTTAAATGGTCAAGGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCa                                                                                                                   <  2:737289/144‑1 (MQ=255)
                                  cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                               >  2:311801/1‑151 (MQ=255)
                                   aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                              >  1:516449/1‑151 (MQ=255)
                                     aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTCAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                            <  2:799830/151‑1 (MQ=255)
                                         tgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                        >  1:828569/1‑151 (MQ=255)
                                           tatgtCGTCTTGCCACCCATCATTATTCAACAGATTAATGACGGAAGATATTTATGGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAATGGTAAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                        <  2:966206/149‑1 (MQ=255)
                                            atgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                      >  2:69650/1‑150 (MQ=255)
                                                 gACTTGCCAGCCATCATTATTCAAAAGATTAAGGACGGCAGATATTTATGGGAAGGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                                  <  2:717673/149‑1 (MQ=255)
                                                  atttgcacgccaTCTTTTTTAAACAGATTAATGACGGAAATTTTTTATCGGTAGGAAGAAAATGCAGAGCTTTTCTCCCTTTCAAATGGTCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCTTTTGCGGGCATTa                                                                                                  <  2:330891/140‑1 (MQ=255)
                                                      gccagccaTCATTATTCAACAAATAAATGACGGCAAATATTTATCGGTAAGGAGAAAATGCAAAGCTTTTCGTCACGTAAAATGGACAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                           <  2:76756/151‑1 (MQ=255)
                                                        cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                         >  2:468793/1‑151 (MQ=255)
                                                            catcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                     >  1:114631/1‑151 (MQ=255)
                                                             atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                    >  2:280543/1‑151 (MQ=255)
                                                                attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                 >  2:517579/1‑151 (MQ=255)
                                                                attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                 >  2:67928/1‑151 (MQ=255)
                                                                attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                 >  2:168581/1‑151 (MQ=255)
                                                                   attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                              >  2:624095/1‑151 (MQ=255)
                                                                    ttCAACAAATTAATGACGGCAGATTTTTATCGGTAAGGAGAAAAGGCAGAGCTTTTCGTCACGTTAAATGGTAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                     <  2:957710/143‑1 (MQ=255)
                                                                       aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAATGGTCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                           <  2:426380/150‑1 (MQ=255)
                                                                          aGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                                                                        >  2:458919/1‑150 (MQ=255)
                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                     >  2:749975/1‑151 (MQ=255)
                                                                            aTTAAGGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACTTTCAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaaa                                                                      <  1:158519/150‑1 (MQ=255)
                                                                              tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                    >  1:393593/1‑150 (MQ=255)
                                                                                aTGACGAGAGATTTTTATCGGGAAGGAGGAACGGCACAGTTTTTCTCCCCTTTCATGGGTAAATCCCTCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                 <  2:711491/151‑1 (MQ=255)
                                                                                  gACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                               >  1:638491/1‑151 (MQ=255)
                                                                                   aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt                                                              >  2:6746/1‑151 (MQ=255)
                                                                                     ggCAGATATTTACGGGTAAGGAGAAAATGCAGAGTTTTTCGTCACTTTCAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                                <  2:49065/147‑1 (MQ=255)
                                                                                        aGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCAACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAGTCGGc                                                          >  1:378174/1‑150 (MQ=255)
                                                                                        aGATAGTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTTGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCCTGAACCACAGAAACCAGGATTCg                                                            >  1:883735/1‑148 (MQ=255)
                                                                                         gATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa                                                         >  2:676504/1‑150 (MQ=255)
                                                                                         gATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa                                                         >  2:732016/1‑150 (MQ=255)
                                                                                          atatTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCaa                                                        <  2:326469/150‑1 (MQ=255)
                                                                                              ttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                    <  2:393174/150‑1 (MQ=255)
                                                                                               tATGGGTAAGGAGAAAATGCAGACCTATTCGTCACGTTCAATGGCCAGTCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                   <  2:955884/150‑1 (MQ=255)
                                                                                               tATCGGAAAGAAGGAAAGGCAGAGTTTTTCGTCACGTTAAATGTTCATTGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                   <  1:295239/150‑1 (MQ=255)
                                                                                                   ggTAAGGAGGAAAGGAAGAGCTTTTCTTCCAGTTAAAGGTTAATTGCCCCCCCCCGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                               <  2:639393/150‑1 (MQ=255)
                                                                                                      aaGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                             <  2:969715/149‑1 (MQ=255)
                                                                                                        ggagaaAATGCAGAGCTTTTCGTAACGTTCAAGGGCCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                                 <  1:927814/143‑1 (MQ=255)
                                                                                                         gagaaAATGCAGAGCTTTTCGTCACGTTAAATGGTAAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                        <  2:875665/151‑1 (MQ=255)
                                                                                                          aggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                       >  2:498751/1‑151 (MQ=255)
                                                                                                           ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                      >  1:891672/1‑151 (MQ=255)
                                                                                                             aagtccagactgtttCGTCACTTTCAATGTAAACCCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                    <  1:904298/138‑1 (MQ=255)
                                                                                                             aaaTGCAGATCTTTTCGTCCCGTTAAATGGTCAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt                                     <  2:148068/150‑1 (MQ=255)
                                                                                                              aaTGCAGATCTTTTCGTCAGGTTCAAGGGTCAGTGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                        <  1:839526/146‑1 (MQ=255)
                                                                                                                     agCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                            >  1:1003057/1‑151 (MQ=255)
                                                                                                                     actttttCGCCAGGTTCAATGTTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                             <  1:511930/145‑1 (MQ=255)
                                                                                                                          ttCGTCACGTTAAAGGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                        <  2:998343/150‑1 (MQ=255)
                                                                                                                                aCGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACCAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCCTGAACCAGAGAAACCAGGAGTCGGCAACCGGAAGCGCCTATCGGGCGACCGAGTGCGATCGCGcc                  >  1:693040/1‑150 (MQ=255)
                                                                                                                                 cGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccac                >  2:526148/1‑151 (MQ=255)
                                                                                                                                     cAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt            >  1:128943/1‑151 (MQ=255)
                                                                                                                                         gggTCACTACAACACCCCGGCCCCCGCCGATACAAAGAGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAACCCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCTAGTGCGATCGCGCCACcattgac         <  1:964875/150‑1 (MQ=255)
                                                                                                                                                tGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAcc  >  2:902469/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC  >  NC_000913/2984292‑2984585

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: