Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,317,978 A→C 6.2% intergenic (‑83/+28) phnL ← / ← phnK ribophosphonate triphosphate synthase subunit; putative ABC transporter‑related ATPase/carbon‑phosphorus lyase complex subunit, putative ATP transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,317,9780AC 6.2% 401.9 / 11.2 133intergenic (‑83/+28)phnL/phnKribophosphonate triphosphate synthase subunit; putative ABC transporter‑related ATPase/carbon‑phosphorus lyase complex subunit, putative ATP transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (58/66);  new base C (3/6);  total (61/72)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 142 total aligned reads displayed.

CGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCG  >  NC_000913/4317829‑4318127
                                                                                                                                                     |                                                                                                                                                     
cGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATaa                                                                                                                                                       >  2:585438/1‑150 (MQ=255)
cGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATaa                                                                                                                                                       >  2:528571/1‑150 (MQ=255)
  ggCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAgg                                                                                                                                                     >  2:804663/1‑150 (MQ=255)
   gCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGc                                                                                                                                                    >  2:645765/1‑150 (MQ=255)
     aGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGt                                                                                                                                                  <  1:1004279/150‑1 (MQ=255)
      ggCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGtt                                                                                                                                                 <  2:1308023/150‑1 (MQ=255)
      ggCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGtt                                                                                                                                                 <  1:429950/150‑1 (MQ=255)
        cgcACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTa                                                                                                                                               <  1:641350/150‑1 (MQ=255)
          cACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                             >  2:1305111/1‑150 (MQ=255)
            cGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                           <  2:1359465/150‑1 (MQ=255)
             gCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                          <  1:1410187/150‑1 (MQ=255)
              ccGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTTTACGTTAATCATTTTTCGCCTCGTGGTTCTGGTGTCGGATGCGACGCTGACGAGTCTTATACGACCTAAGGGGAGCGCATTTGTAGGCCGGATAATGCGTTTAcgcctc                                                                                                                                         >  2:1234931/1‑148 (MQ=255)
              ccGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                         >  2:10897/1‑150 (MQ=255)
              ccGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                         >  2:446102/1‑150 (MQ=255)
                gcattGCTGGGGCAGGCTTAAGGTTTTACTGTCGTTTTGTCCCTTAATCATTTTTCGCCTCGTGGATCGGGTGGCGGATGCGAGGTGGACGCGTCTAATCGGACGTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                       <  2:723677/147‑1 (MQ=25)
                   ttGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                    >  2:837067/1‑150 (MQ=255)
                   ttGCTGGGGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGGGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                    <  2:252291/150‑1 (MQ=37)
                    tGCTGGTGCAGGATGCAGGTTTTACTGACGTTTTGTACGTTAATCATTTATCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                   <  1:1269496/150‑1 (MQ=37)
                    tGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                   >  1:891272/1‑150 (MQ=255)
                    tGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGAGGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                   <  2:1446155/150‑1 (MQ=37)
                     gCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                  <  2:800118/150‑1 (MQ=37)
                     gCTGGGGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                  <  2:900536/150‑1 (MQ=37)
                           gCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGc                                                                                                                            <  1:636844/150‑1 (MQ=255)
                              ggATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGTGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGAGTTTACGCCGCATCCGGCACCGGCTCa                                                                                                                         >  2:729704/1‑150 (MQ=255)
                               gATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCaa                                                                                                                        >  1:8348/1‑150 (MQ=255)
                               gATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCaa                                                                                                                        >  1:585599/1‑150 (MQ=255)
                                   aaGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCt                                                                                                                    >  1:1122983/1‑150 (MQ=255)
                                    aGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTg                                                                                                                   <  2:62312/150‑1 (MQ=255)
                                    aGGTTTTACGGACGTTATGTAGGTTAATCATTTTTCGCCTCGGGGTTGGGGTGTCGGATGCGACACTGACGCGTCTTATCCGACCTACGGGGAGCGCAGTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCGCAATTCTg                                                                                                                   <  2:1698211/150‑1 (MQ=25)
                                     ggTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGc                                                                                                                  >  1:1125165/1‑150 (MQ=255)
                                     ggTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGc                                                                                                                  <  2:337718/150‑1 (MQ=255)
                                      gTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATACGACCTACGGGGAGCGCATTTGTAGGCCGGATACGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCa                                                                                                                 >  1:121750/1‑150 (MQ=255)
                                        tttACTGACTTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCaaa                                                                                                               <  2:1214005/150‑1 (MQ=255)
                                            cTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGTCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCg                                                                                                           >  2:1377893/1‑150 (MQ=255)
                                              gACGTTTTGTACGTTAATCATTTTTCGCCTGGGGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTAGGGGGAGCGCATTTGTAGGCCGTATAAGGCGTTTACGCCGCATCGGGCACCGGCTCAATTCTGCAAAACCGAt                                                                                                         <  1:404208/150‑1 (MQ=255)
                                              gACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCGCAATTCTGCAAAGCCGAt                                                                                                         >  2:250214/1‑150 (MQ=255)
                                               aCGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATg                                                                                                        <  1:717238/150‑1 (MQ=255)
                                                 gTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATACGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATgcc                                                                                                      >  2:1145715/1‑148 (MQ=255)
                                                    ttGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGTTAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACAcc                                                                                                   <  1:397488/150‑1 (MQ=255)
                                                    ttGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACAcc                                                                                                   >  2:765871/1‑150 (MQ=255)
                                                       tACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagc                                                                                                >  1:865727/1‑150 (MQ=255)
                                                       tACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGAGAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACCCcagc                                                                                                >  1:906279/1‑150 (MQ=255)
                                                        aCGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagcc                                                                                               >  1:1788410/1‑149 (MQ=255)
                                                        aCGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagca                                                                                               >  2:366872/1‑150 (MQ=255)
                                                               tCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCCCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTgtgt                                                                                        >  1:1036435/1‑150 (MQ=255)
                                                                cATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGta                                                                                       <  2:344047/150‑1 (MQ=255)
                                                                cATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGta                                                                                       <  2:22957/150‑1 (MQ=255)
                                                                    tttCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACg                                                                                   <  2:565384/150‑1 (MQ=255)
                                                                    tttCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACg                                                                                   <  2:1641054/150‑1 (MQ=255)
                                                                    tttCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACg                                                                                   >  2:718885/1‑150 (MQ=255)
                                                                             gTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCg                                                                          <  2:197625/150‑1 (MQ=255)
                                                                             gTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCg                                                                          >  2:390148/1‑150 (MQ=255)
                                                                               ggTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCggg                                                                        <  1:788949/150‑1 (MQ=255)
                                                                               ggTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCggg                                                                        >  1:1215261/1‑150 (MQ=255)
                                                                                   cGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGgtcgt                                                                    >  2:1710957/1‑150 (MQ=255)
                                                                                     ggTGTCGGATGTGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGtcgtcg                                                                  >  2:1714924/1‑150 (MQ=255)
                                                                                     ggTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGtcgtcg                                                                  <  2:707367/150‑1 (MQ=255)
                                                                                       ggtCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAg                                                                <  1:1305111/149‑1 (MQ=255)
                                                                                        gtCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGCTGCGGGTCGTCgcgc                                                               >  2:892258/1‑150 (MQ=255)
                                                                                         tCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCa                                                              <  1:705618/150‑1 (MQ=255)
                                                                                          cGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCAc                                                             >  2:581387/1‑150 (MQ=255)
                                                                                          cGGATGCGACGCTGACGCGTCTAATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCAc                                                             >  1:1668762/1‑150 (MQ=255)
                                                                                            gATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCAcgc                                                           <  1:922523/150‑1 (MQ=255)
                                                                                              tGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGcgg                                                         >  2:677885/1‑150 (MQ=255)
                                                                                                cGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGcggtc                                                       >  2:1115309/1‑150 (MQ=255)
                                                                                                 gACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGcggtcg                                                      >  2:406679/1‑150 (MQ=255)
                                                                                                   cGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGcggtcggt                                                    >  1:217232/1‑150 (MQ=255)
                                                                                                      tGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTaa                                                 <  1:1712845/150‑1 (MQ=255)
                                                                                                       gACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAAc                                                <  1:528571/150‑1 (MQ=255)
                                                                                                       gACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAAc                                                <  1:1271032/150‑1 (MQ=255)
                                                                                                       gACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAAc                                                <  1:133910/150‑1 (MQ=255)
                                                                                                           cgTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCAc                                            <  1:804663/150‑1 (MQ=255)
                                                                                                           cgTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCAc                                            <  1:60690/150‑1 (MQ=255)
                                                                                                           cgTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCAc                                            <  2:1122983/150‑1 (MQ=255)
                                                                                                             tCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCActc                                          <  1:727808/150‑1 (MQ=255)
                                                                                                                   ccGACCTACGGTGAGCGCATTTGTAGGCCGGATAAGGTGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCATCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTccacc                                    <  1:729704/150‑1 (MQ=255)
                                                                                                                        cTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTg                               >  2:449485/1‑150 (MQ=255)
                                                                                                                         tACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGc                              >  2:951800/1‑150 (MQ=255)
                                                                                                                          aCGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGcc                             <  1:918570/150‑1 (MQ=255)
                                                                                                                          aCGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGcc                             >  2:1086070/1‑150 (MQ=255)
                                                                                                                            ggggAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACTCCAGCAGCGGTGTATACGGATGATGCGGGTCGTCGAGCCCGCGGTCGGTTCACCCACTCTCCACCACGTGcccc                           >  1:158231/1‑150 (MQ=39)
                                                                                                                            ggggAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTACCCCACTCTCCCCCACTTGcccc                           >  2:957161/1‑150 (MQ=255)
                                                                                                                            ggggAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGcccc                           >  2:1658262/1‑150 (MQ=255)
                                                                                                                             gggAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCt                          >  2:1470626/1‑150 (MQ=255)
                                                                                                                                aGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCt                       <  1:1377893/150‑1 (MQ=255)
                                                                                                                                     aTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCAGCTtcatc                  <  1:1170646/150‑1 (MQ=255)
                                                                                                                                      tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTtcatca                 <  2:743808/150‑1 (MQ=255)
                                                                                                                                        tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCAcc               >  1:1883762/1‑150 (MQ=255)
                                                                                                                                        tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCAcc               >  2:1033730/1‑150 (MQ=255)
                                                                                                                                         cggtgCCGGATGCGGCGTAAACGCCTTATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCGccc              >  1:984018/5‑149 (MQ=255)
                                                                                                                                              ccGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCaaa         >  1:1370784/1‑150 (MQ=255)
                                                                                                                                               cGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAAc        <  1:677885/150‑1 (MQ=255)
                                                                                                                                               cGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAAc        <  1:1374516/150‑1 (MQ=255)
                                                                                                                                                ggATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCCGCAAACg       >  1:394411/1‑150 (MQ=255)
                                                                                                                                                 gATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACgg      <  2:207411/150‑1 (MQ=255)
                                                                                                                                                    aaGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGCGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCCTCACCAGCAAACGGTcc   >  2:765036/1‑150 (MQ=255)
                                                                                                                                                     aGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCGCCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCg  >  1:1066505/1‑150 (MQ=255)
                                                                                                                                                     aGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCg  >  1:1601706/1‑150 (MQ=255)
                                                                                                                                                     aGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCg  <  1:1434632/150‑1 (MQ=255)
                                                                                                                                                     aGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCg  >  2:1318073/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGGCAGGCGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCG  >  NC_000913/4317829‑4318127

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