Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4410AC44.2% 131.0 / 23.6 86G384G (GGT→GGGyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (44/4);  new base C (0/38);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 3.00e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGA  >  NC_000913/2984297‑2984588
                                                                                                                                                 |                                                                                                                                                   
tGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                               >  2:172233/1‑151 (MQ=255)
  tCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                                                                             >  1:437429/1‑151 (MQ=255)
     tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccc                                                                                                                                          >  1:423225/1‑151 (MQ=255)
     tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccc                                                                                                                                          >  2:193793/1‑151 (MQ=255)
      gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccg                                                                                                                                         <  1:851448/151‑1 (MQ=255)
         aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgccg                                                                                                                                      >  1:1369750/1‑151 (MQ=255)
         aaTGAAAAAATGAGGGAGCGCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCGCACCCTGGCCCccaccg                                                                                                                                      <  2:912193/151‑1 (MQ=255)
              aaaaaTGAGGGAGCCCAGAAAATGTATGTCCACTTGCCAGCCATCATTATTTAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATAc                                                                                                                                  >  1:629359/1‑150 (MQ=255)
               aaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGTAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACa                                                                                                                                 >  2:73165/1‑150 (MQ=255)
               aaaaTAAGGGAGCCAAGAAAAGTTATGTCGACTTGCCAGCCATCATTATTCAACAGTTTAATGACGGCAGATTTTTATCGTAAAGGAGGAAAGGCAGAGCTTTTGTTCACTTTCAATGGTAAATGCCACCCCCTGCCCCCCGCCGATACaa                                                                                                                                <  2:460913/151‑1 (MQ=255)
                    gAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGt                                                                                                                           >  2:968462/1‑151 (MQ=255)
                    gAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGt                                                                                                                           >  1:344646/1‑151 (MQ=255)
                    gAGGGAGCCCAAAAAATGTTTGTCGACTGGCCACCCATCATTATTCAACAGATTAAGGCCGGCAGATTTTAACGGGAAAGGAGGAAATGCGGAGCTTTTCTCCTCTTTCAAGGGCAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCg                                                                                                                            <  1:1062109/150‑1 (MQ=255)
                         aGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGcc                                                                                                                       >  1:566915/1‑150 (MQ=255)
                           cccAGAAAATGTATGTCGACTTACCAGCCATCATTATTCAACAGATTAATGACGGGATATATTTATCGATAAGGAGAAAATGCAGAGCTATTCGTCAGGTGAAATGGTCATCGCCACACCCGGGCCCCCGCCGATACAAAGCGTTGCCa                                                                                                                      <  2:841940/149‑1 (MQ=255)
                            ccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                   >  2:1313767/1‑151 (MQ=255)
                            ccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                   >  2:1073571/1‑151 (MQ=255)
                             cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGt                                                                                                                    >  2:427365/1‑149 (MQ=255)
                             cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                  >  1:968004/1‑151 (MQ=255)
                                aaaaTGTATGTCGACTGCCCATCAATCATTATAAAAAAGATTAATGCGGTCGGATTTTTATCGGAAGGGAGAAAATGAAGAGCTTTTCGTCACTTTCAATGGTCAACCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCtt                                                                                                                <  1:1119006/150‑1 (MQ=255)
                                   atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                            >  2:1258429/1‑151 (MQ=255)
                                    tgtatgtCGACTTGCCAGTCATGATTATTCAACAGATTAATGACGGCAGATATTTATCGGTACGGCGGAAATGCAGAGCTATTCGTCACGATCAATGGTCAATGCCACACCCTGACCCCCGCCAATACAAAGCGTTGCCACCCCTTTGc                                                                                                             >  2:1001351/1‑149 (MQ=255)
                                     gtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                                                                          >  1:967129/1‑151 (MQ=255)
                                           cGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                    >  2:623540/1‑151 (MQ=255)
                                            gACTGGCAAGCCATAATAATCCAACAAATAAATGACTGCAGATTTTTAACGGAAAGAGGAAAATCAAGATTATTTGGCCACTTTAAATAGTAAATGCCCCCCATTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                    <  2:441058/150‑1 (MQ=255)
                                               ttcccagccaTCATTATTCAAAAGTTAAATGCCGGAAGATATTTATCGGTAAGGGGAAAATCCAGAGTTATTCGTCAGGTTAATTGGCCATTCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                           <  2:1224695/137‑1 (MQ=255)
                                                tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                               >  1:1231794/1‑151 (MQ=255)
                                                tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                               >  1:507779/1‑151 (MQ=255)
                                                   cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                            >  1:1254138/1‑151 (MQ=255)
                                                        atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                       >  2:1208699/1‑151 (MQ=255)
                                                         tcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                       >  1:993196/1‑150 (MQ=255)
                                                            ttattaAACAGATTAATGACGGAATATATTTATGGGAAAGGAGGAAATGCAGAGCTTTTCGTCATGTTAAATGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                     <  1:148745/149‑1 (MQ=255)
                                                            ttattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                     >  2:148704/1‑149 (MQ=255)
                                                            ttatcaaaCAGATTAATGACGGAAGATATTTATCGGTAAGGAGAAAATGAAGATCTTTTCGTCACTTTAAATGGTCAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                       <  1:984017/141‑1 (MQ=255)
                                                                tCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                >  2:323188/1‑150 (MQ=255)
                                                                  aaCAGATTAATGCCGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCACGTTAAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                               <  1:72082/149‑1 (MQ=255)
                                                                  aaCAGATTAATGACGGAAGATATTTATCGGTAGGGAGGAAATGCAGAGCTATTCGTCACGTTCAAGGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                               <  2:425544/149‑1 (MQ=255)
                                                                   aaaGATTAATGACGGAAGATTTTTATGGGAAGGGAGAAAATGCAGAGTTTTTCTCCCGGTTCAATGGTAAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag                                                                            <  2:622390/149‑1 (MQ=255)
                                                                   aCAAATTAATGCCGGCAGATATTTAGCGGTAAGGGGAAAATTCAGAGCTTTTCGACCCGTTCAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag                                                                            <  2:1188315/151‑1 (MQ=255)
                                                                       aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGAAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTAGCCCCCGCCGATACAAAGCGTTGCCAGTCATTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                        >  1:1165984/1‑151 (MQ=255)
                                                                        ttAATGACGGCAGATATTTATGGATAAGGAGGAAATGAAGGGTTTTTCGTCAGGTTTAATTGTCATTGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaaa                                                                         <  2:1370160/149‑1 (MQ=255)
                                                                        ttAATGACGGCAGATATTTAGGGGTAAGGAGAAAAGGCAGAGTTATTCGCCACGTTAAATGGTCAAGCCCCCCCCCTGCCCTCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaaa                                                                         <  2:74526/149‑1 (MQ=255)
                                                                         tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                       >  1:706487/1‑150 (MQ=255)
                                                                         tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                       >  1:620780/1‑150 (MQ=255)
                                                                          aaTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTAAAGGGTCAACGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                      <  1:1259153/150‑1 (MQ=255)
                                                                          aaTGACGGCAGATATTTATCGGTAAGGAGAAAAGGCAGAGCTATTCGTCACGTTAAATGGTAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAg                                                                     <  1:999380/151‑1 (MQ=255)
                                                                           atgacgtgaattatTTATCGGTAGAGAGAAAAGACAGCGCTTTTCTTAACTTTTTTTGGTCATTGCCAACCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                       <  2:1288998/137‑1 (MQ=255)
                                                                                ggCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                               >  1:832563/1‑151 (MQ=255)
                                                                                  catatatTTATCGAAAGGGAGAAAGTGCAGGTCTTTTCGTCACTTTAAATGGTAACTCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                               <  2:1143112/146‑1 (MQ=255)
                                                                                  cAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                             >  1:554793/1‑151 (MQ=255)
                                                                                   aGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCCCACCCTGGCCCCCGCCGATACAAAGCGTTGCCCGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                             >  2:768824/1‑150 (MQ=255)
                                                                                     atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc                                                          >  2:111889/1‑151 (MQ=255)
                                                                                        tttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                       >  1:169524/1‑151 (MQ=255)
                                                                                        tttATCGGTAAGGAGAAAAGGAAGAGCTTTTCGTCACGTTAAAGGGTAAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcag                                                        <  1:603404/150‑3 (MQ=255)
                                                                                        tttATCGATAGGGAGAAAGTGCAGAGTTTTTCTTACCGTTAAAGGGTAACTGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc                                                          <  2:447269/148‑1 (MQ=255)
                                                                                         ttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGTCATTACGTTTCACCATTTCATGAACCAGCGAAACCAGGATTCGGCAACCgg                                                       >  2:747593/1‑150 (MQ=255)
                                                                                         ttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                       >  1:1344264/1‑150 (MQ=255)
                                                                                           atgggaaaGGAGGAATGGCAGATCTTTTCGTCACTTTCAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                      <  2:799744/143‑1 (MQ=255)
                                                                                              ggaaaGGAGAAAAGGCAGTTCTCTTCGTAAGTTTAATTGTCCAATCCCCCACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCGTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                  <  2:1251199/147‑1 (MQ=255)
                                                                                                tAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                               >  1:1213011/1‑151 (MQ=255)
                                                                                                  aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                             >  1:1230140/1‑151 (MQ=255)
                                                                                                     agaaaaTGCAGAGCTTTTCGTCCCTTTAAATGGTAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                           <  1:1373701/147‑1 (MQ=255)
                                                                                                        aaaTGCAGAGTTTTTCGTCCCTTTAAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt                                        <  1:160275/150‑1 (MQ=255)
                                                                                                        aaaTGAAGAGCTATTCGTCACGTTCAATGGTCAAGGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                          <  1:11002/148‑1 (MQ=255)
                                                                                                         aaTGAAGAGCTTTTCGTCACGTTCAATGGTCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGa                                      <  1:1237488/151‑1 (MQ=255)
                                                                                                               gaccttttCGTCACGTTAAAGGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAg                                 <  2:1368996/144‑1 (MQ=255)
                                                                                                                agCTATTCGTCACTTTAAATGGTCACCCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                               <  2:1336662/151‑1 (MQ=255)
                                                                                                                 gttATTCGCCTCTTTCAAGGATCAAGCCCCCCCCCTCGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                               <  1:427464/148‑1 (MQ=255)
                                                                                                                   ttttCGTCACGTTCAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                              <  2:441323/147‑1 (MQ=255)
                                                                                                                   ttttCGTCACGTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                              <  1:745546/147‑1 (MQ=255)
                                                                                                                   ttttCGTCACGTTAAATGGTCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                              <  1:628955/147‑1 (MQ=255)
                                                                                                                    aTTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                           >  1:809140/1‑151 (MQ=255)
                                                                                                                     ttCGTCACTTTCAATGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                              <  1:455021/147‑1 (MQ=255)
                                                                                                                     ttCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                          >  2:1144204/1‑151 (MQ=255)
                                                                                                                        gTCACTTTCAATGGTCAGCGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                         <  2:422460/149‑1 (MQ=255)
                                                                                                                        gTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg                       >  1:450335/1‑151 (MQ=255)
                                                                                                                        gTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                        >  2:891585/1‑150 (MQ=255)
                                                                                                                          cACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                     >  1:70135/1‑151 (MQ=255)
                                                                                                                              ttaaaTGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcca                    <  1:1062833/145‑1 (MQ=255)
                                                                                                                                 attGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccaca                  <  2:361788/145‑2 (MQ=255)
                                                                                                                                 aagggTCATTCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg              <  1:1084341/148‑1 (MQ=255)
                                                                                                                                 aaTGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt               <  1:961384/150‑1 (MQ=255)
                                                                                                                                  aTGGTCAATGCCACACCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg              >  2:961591/1‑150 (MQ=255)
                                                                                                                                    gttCAGTGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat                <  2:344555/144‑1 (MQ=255)
                                                                                                                                         aaTGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattga       <  1:1337787/151‑1 (MQ=255)
                                                                                                                                           tGCCACACCCTGCCACCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAATACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGCTCGCGCCACCattgacattgac      >  2:906215/1‑150 (MQ=255)
                                                                                                                                              cacaCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGa  <  1:705803/151‑1 (MQ=255)
                                                                                                                                               acacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGa  >  2:929632/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGA  >  NC_000913/2984297‑2984588

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: