Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC46.0% 97.5 / 21.3 63G384G (GGT→GGG) yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/7);  new base C (0/29);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.12e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAC  >  NC_000913/2984292‑2984589
                                                                                                                                                       |                                                                                                                                                    
gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAACGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCacac                                                                                                                                                     <  2:224382/151‑1 (MQ=255)
  atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACAccc                                                                                                                                                   >  2:662124/1‑151 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTg                                                                                                                                                 >  1:1005158/1‑151 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTg                                                                                                                                                 <  1:607691/151‑1 (MQ=255)
    atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTg                                                                                                                                                 >  1:364614/1‑151 (MQ=255)
      gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGc                                                                                                                                               <  1:1069358/151‑1 (MQ=255)
           gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCccg                                                                                                                                          <  1:410810/151‑1 (MQ=255)
                  gAAACATGAGGGAGCCCAGAAAATGTATGTCTACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGTAGGAAATGCAGTGCTATTCGTCTCGTTCAA‑TTGTCTATGCCTCACCCTGGCCCccg                                                                                                                                          <  1:1022384/144‑1 (MQ=255)
                         gAGGGAGCCAAGAAAATGTATGTCGACTTG‑CCAGCCATAATTATTCAACAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTAAA‑TGGTCAATCCCACCCCCTGGCCCCCGCCGATACAAAGc                                                                                                                              <  2:99819/149‑1 (MQ=255)
                             gAGCCCAAAAAATGTATGTCG‑CTTGCCCGGCCATCATTATTCAACAGATTAATGCCGGCAGATATTTATCGGTAAGGAGAAAAGGAGGAGTTTTTCTTCACTTTCAA‑GGGTAAAGCCCCCACCCGGGCCCCCGCCGATACAAAGCGTTGcc                                                                                                                        <  2:692043/151‑1 (MQ=255)
                                  cAGAAAATGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                   >  2:75165/1‑151 (MQ=255)
                                   aaaaaaTGTATGTCGACTTG‑CCACCCATCATTATACAACAGATTAAGGACGGCAGATTTTTATGGGTAAGGAGAAAATGCAGAGCTATTCGTCCCGTTCAA‑TGGTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                    <  1:389364/147‑1 (MQ=255)
                                   aaaaaaTGTATGTCGACTGG‑CCAGCCATCATTATCAAACAGATAAAGGAGGGAAATTATTTAGCGGGAGGGAGAAAATGAAGAGCTTTTCGTCACGTCCAA‑TGGTAAAGGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAg                                                                                                                      <  1:61610/145‑1 (MQ=255)
                                   aGAAAATGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                  >  2:39391/1‑151 (MQ=255)
                                     aaatTGTATGTCGACTTG‑CCAGCCATAATTATTAAAAAGATAAATGACGGAAGATATTTATGGGAAGGGAGAAAATGCAGAGCTTTTCTCCACGTACAA‑TGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                  <  2:1028793/149‑1 (MQ=255)
                                      aaaTGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                               >  1:850335/1‑151 (MQ=255)
                                      aaaTGTATGCCGACTTG‑CCAGCCTTCATTATTCAAAAAATTAATGAGGGAAGATTTTTAGGGGTAAGGAGAAAAGGAAGAGCTTTTCGTCACTTTAAA‑TGGCAAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                               <  1:1136285/151‑1 (MQ=255)
                                           tatgtCGACTTG‑CCAGCCATCATTATTCAAAAGATTAATGACGGCAGATTTTTAGCGGAAAGGAGAAAAGGCAGACCTTTTCGTCACGTTAAA‑TGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                             <  1:270397/148‑1 (MQ=255)
                                              gtCGACTTG‑CCACCCATCTTTATTAAACAGATTAATGACGGCAGATATTTATCGGAAAGGAGGAAATGCAGAGCTTTTCGTCACTTTAAA‑GGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                                                                        <  2:766226/150‑1 (MQ=255)
                                               tCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                       >  1:1048382/1‑150 (MQ=255)
                                               tCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                       >  2:1120557/1‑150 (MQ=255)
                                                cGACTTG‑CCATCAAGAATTATACAACAAATAAATGGCGGAAAATATTTAGAGGGAGGGAAAAAATCCAGATGTATTCGTCATTTTCAA‑TGTTCACCCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                     <  2:62285/151‑1 (MQ=255)
                                                   atttgccgccaatcatTATTAAAAAAATTAATGACGGCAAATATTTAGCGGAAAGGAGAAAAGGAAGGGTTTTTGGTCACGTTAAA‑GGGTAAAGGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                          <  2:90998/133‑1 (MQ=255)
                                                    tgg‑ccagcaaATATTATTAAACAGATAAATGCCGGAAAATATTTAGCGGAAAGGAGAAAATGAAGGCCTTTTCGTCACGTTAAA‑GGGTAAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                            <  1:130564/138‑1 (MQ=255)
                                                     tG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGTTTTTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCCCTTTCAA‑TGGCCAATCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                  <  1:850485/149‑1 (MQ=255)
                                                     tG‑CCAGCCATCATTATTCAACAGATTAAGGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAA‑TGGTCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                  <  2:225655/149‑1 (MQ=255)
                                                       ccaagccaTATTTATACAACAAATTAATGACGAAAGATATTTAGCGGAAAGGAGGAAATGAAGATCTTTTCGTCCCGTAAAA‑TGGTAAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                               <  1:426346/148‑1 (MQ=255)
                                                         cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                             >  1:694015/1‑151 (MQ=255)
                                                            ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                          >  2:656074/1‑151 (MQ=255)
                                                              atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCCCACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                         <  1:121940/150‑1 (MQ=255)
                                                                catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                      >  2:194365/1‑151 (MQ=255)
                                                                catTATCCAACAGATGATTGTCGAAATATTTTTATCGGAAAGGAGGATTTTCAGGCTTCTTCATAACTTCCAA‑TGTTAATTACCACCCCCGGGCCACCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAt                                                                                           <  2:509830/146‑1 (MQ=255)
                                                                 attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                     >  1:1093310/1‑151 (MQ=255)
                                                                  ttattaAAAAGATTAATGACGGCAATTTTTTATGGGAAAGGAGAAAATGCAGAGTTTTTCGTCACGTTAAA‑GGGTCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                      <  1:1034157/149‑1 (MQ=255)
                                                                    attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                  >  2:610223/1‑151 (MQ=255)
                                                                     ttCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                 >  1:64792/1‑151 (MQ=255)
                                                                       cAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                 >  2:175571/1‑149 (MQ=255)
                                                                        aaCAGTTTAAGGACGGCAGATTTTTATCGGTAAGGAGAAAAGGCAGAGTTTTTCTTCCCGTTCAA‑GGTTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                <  2:952554/149‑1 (MQ=255)
                                                                         aCATATTATTGACGGCAAATATTTAGCGGAAAGGAGAAAATGCAGAGTTTTTCTTCACTTAAAG‑TGGACAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag                                                                             <  1:563402/151‑1 (MQ=255)
                                                                                     cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                                 >  2:383910/1‑151 (MQ=255)
                                                                                      ggaaGATATTTATCGGTAAGAAGAAAATGCAGAGTTATTCGCCACTTTCAA‑TGTTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                                    <  2:355448/144‑1 (MQ=255)
                                                                                         aGATATTTATCGGTAGGGAGGAAATGCAGAGTTTTTCGTACCTTTCAA‑TGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                               <  1:655521/149‑1 (MQ=255)
                                                                                              tttAGCGGAAGGAAGAAAATGGAGAGCTTTTCCTC‑TTTTCAACTGCTCAACCCCACACCCCGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa                                                             <  2:1099083/146‑1 (MQ=255)
                                                                                                      tAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                                >  1:1085378/1‑151 (MQ=255)
                                                                                                       aaGGAGAAAATGCAGAGTTTTTCGTCACTTTAAA‑TGGCAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                                <  2:12389/150‑1 (MQ=255)
                                                                                                          gagaaAATGCGGATCTTTTCGTAACTTTAAA‑TGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                                <  1:1031736/147‑1 (MQ=255)
                                                                                                            ggaAGTGCAGACCTTTTCGTCATTTACAA‑TGGTAAGTGCCGCCCCGTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                           <  1:723120/150‑1 (MQ=255)
                                                                                                                aTGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                                      >  2:576522/1‑151 (MQ=255)
                                                                                                                  gCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                    >  2:115202/1‑151 (MQ=255)
                                                                                                                    agagCTTTTCGTCACGTTAAA‑TGGTAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                    <  2:121709/149‑1 (MQ=255)
                                                                                                                     gagCTTTTCGTCACGTACAA‑TGGTAAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                     <  1:348415/147‑1 (MQ=255)
                                                                                                                       gCTTTTCGTCACGTTCAA‑GGGTAAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                               <  1:1104339/151‑1 (MQ=255)
                                                                                                                       gCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                               >  2:1087843/1‑151 (MQ=255)
                                                                                                                            tCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                          >  1:510828/1‑151 (MQ=255)
                                                                                                                              gTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg                        >  2:1139525/1‑151 (MQ=255)
                                                                                                                                cACGTTCAA‑TGGCCAATCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                         <  2:243971/148‑1 (MQ=255)
                                                                                                                                 aCGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                      >  1:37328/1‑150 (MQ=255)
                                                                                                                                  acgtaaaaggtgaaaTGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                       <  2:1075226/136‑1 (MQ=255)
                                                                                                                                    ttCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                   >  1:566286/1‑150 (MQ=255)
                                                                                                                                     tCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat                 >  1:604885/1‑151 (MQ=255)
                                                                                                                                      cAA‑TGGTCAATGCCACACCCTGTCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt                >  1:784529/1‑151 (MQ=255)
                                                                                                                                            gtcatgcccACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattga              <  2:1121113/141‑1 (MQ=255)
                                                                                                                                               aaTGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg         <  2:136360/150‑1 (MQ=255)
                                                                                                                                               aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg         >  1:280054/1‑150 (MQ=255)
                                                                                                                                                     acacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc  >  2:184896/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTG‑CCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAA‑TGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAC  >  NC_000913/2984292‑2984589

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