Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,317,985 T→A 5.9% intergenic (‑90/+21) 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,9850TA 5.9% 422.4 / 11.2 136intergenic (‑90/+21)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 T (64/64);  new base A (2/6);  total (66/70)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 137 total aligned reads displayed.

CGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCA  >  NC_000913/4317837‑4318133
                                                                                                                                                    |                                                                                                                                                    
cgcACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTa                                                                                                                                                     <  1:641350/150‑1 (MQ=255)
  cACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                                   >  2:1305111/1‑150 (MQ=255)
    cGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                 <  2:747433/150‑1 (MQ=255)
      ccGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                               >  2:446102/1‑150 (MQ=255)
      ccGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                               >  2:788949/1‑150 (MQ=255)
       cGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                              <  1:702085/150‑1 (MQ=255)
       cGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                              <  2:1300220/150‑1 (MQ=255)
        gcattGCTGGGGCAGGCTTAAGGTTTTACTGTCGTTTTGTCCCTTAATCATTTTTCGCCTCGTGGATCGGGTGGCGGATGCGAGGTGGACGCGTCTAATCGGACGTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                             <  2:723677/147‑1 (MQ=25)
            tGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGAGGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                         <  2:1446155/150‑1 (MQ=37)
            tGCTGGTGCAGGATGAAGGTTTTACTGACGATTTGTACCTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                         <  1:1694061/150‑1 (MQ=37)
             gCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                        <  2:800118/150‑1 (MQ=37)
             gCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                        <  2:347462/150‑1 (MQ=255)
             gCTGGGGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                        <  2:900536/150‑1 (MQ=37)
               tGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTACTCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAc                                                                                                                                      <  1:1380996/150‑1 (MQ=255)
                      ggATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCa                                                                                                                               >  1:325014/1‑150 (MQ=255)
                       gATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCaa                                                                                                                              >  1:585599/1‑150 (MQ=255)
                       gATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCaa                                                                                                                              >  1:8348/1‑150 (MQ=255)
                            aGGTTTTACGGACGTTATGTAGGTTAATCATTTTTCGCCTCGGGGTTGGGGTGTCGGATGCGACACTGACGCGTCTTATCCGACCTACGGGGAGCGCAGTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCGCAATTCTg                                                                                                                         <  2:1698211/150‑1 (MQ=25)
                             ggTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGc                                                                                                                        <  2:337718/150‑1 (MQ=255)
                             ggTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGc                                                                                                                        >  1:1125165/1‑150 (MQ=255)
                              gTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATACGACCTACGGGGAGCGCATTTGTAGGCCGGATACGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCa                                                                                                                       >  1:121750/1‑150 (MQ=255)
                                tttACTGACTTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCaaa                                                                                                                     <  2:1214005/150‑1 (MQ=255)
                                    cTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGTCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCg                                                                                                                 >  2:1377893/1‑150 (MQ=255)
                                    cTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCg                                                                                                                 >  2:486992/1‑150 (MQ=255)
                                    cTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCACATCCGGCACCGGCTCAATTCTGCAAAACCg                                                                                                                 >  2:1059195/1‑150 (MQ=255)
                                     tGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGa                                                                                                                <  1:1106475/150‑1 (MQ=255)
                                      gACGTTTTGTACGTTAATCATTTTTCGCCTGGGGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTAGGGGGAGCGCATTTGTAGGCCGTATAAGGCGTTTACGCCGCATCGGGCACCGGCTCAATTCTGCAAAACCGAt                                                                                                               <  1:404208/150‑1 (MQ=255)
                                      gACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGAt                                                                                                               >  1:1147230/1‑150 (MQ=255)
                                      gACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCGCAATTCTGCAAAGCCGAt                                                                                                               >  2:250214/1‑150 (MQ=255)
                                       aCGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATg                                                                                                              <  1:717238/150‑1 (MQ=255)
                                         gTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGac                                                                                                            <  1:844774/150‑1 (MQ=255)
                                            ttGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACAcc                                                                                                         >  2:765871/1‑150 (MQ=255)
                                               tACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagc                                                                                                      >  1:865727/1‑150 (MQ=255)
                                                aCGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagcc                                                                                                     >  1:1788410/1‑149 (MQ=255)
                                                aCGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagca                                                                                                     >  2:366872/1‑150 (MQ=255)
                                                 cGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTCGGCCGGATAAGGCGATTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagcag                                                                                                    >  1:1342238/1‑150 (MQ=255)
                                                    tAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCtg                                                                                                 >  2:1759088/1‑150 (MQ=255)
                                                    tAATCATTTTTCGCCTCGTGGTGCGGGTGGCGGAGGCGACGCGGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCtg                                                                                                 <  2:309401/150‑1 (MQ=255)
                                                      aTCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCtgtg                                                                                               >  1:1817106/1‑150 (MQ=255)
                                                       tCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCCCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTgtgt                                                                                              >  1:1036435/1‑150 (MQ=255)
                                                          tttttCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGtata                                                                                           <  2:696507/150‑1 (MQ=255)
                                                          tttttCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGtata                                                                                           <  2:1090263/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)
                                                                       ggTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCggg                                                                              <  2:572206/150‑1 (MQ=255)
                                                                         ttCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGgtc                                                                            <  1:949288/150‑1 (MQ=255)
                                                                           cGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGgtcgt                                                                          >  2:1710957/1‑150 (MQ=255)
                                                                             ggTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGtcgtcg                                                                        <  2:724341/150‑1 (MQ=255)
                                                                                gtCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGCTGCGGGTCGTCgcgc                                                                     >  2:892258/1‑150 (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)
                                                                                        cGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGcggtc                                                             <  2:1530942/150‑1 (MQ=255)
                                                                                          aCGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGcggtcgg                                                           >  2:248309/1‑150 (MQ=255)
                                                                                          aCGCGGACGCGACTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGAAAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGcggtcgg                                                           <  2:1710246/150‑1 (MQ=255)
                                                                                               gACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAAc                                                      <  1:1271032/150‑1 (MQ=255)
                                                                                               gACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAAc                                                      <  1:133910/150‑1 (MQ=255)
                                                                                                   cgTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCAc                                                  <  2:1122983/150‑1 (MQ=255)
                                                                                                   cgTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCAc                                                  <  1:60690/150‑1 (MQ=255)
                                                                                                    gTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCAct                                                 <  1:1254728/150‑1 (MQ=255)
                                                                                                     tCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCActc                                                <  1:720884/150‑1 (MQ=255)
                                                                                                         aTCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTcca                                            <  1:1539672/150‑1 (MQ=255)
                                                                                                          tCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTccac                                           >  1:894094/1‑150 (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)
                                                                                                                    ggggAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTACCCCACTCTCCCCCACTTGcccc                                 >  2:957161/1‑150 (MQ=255)
                                                                                                                    ggggAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGcccc                                 <  1:421909/150‑1 (MQ=255)
                                                                                                                        aGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCt                             <  1:1377893/150‑1 (MQ=255)
                                                                                                                         gcgcATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCtt                            <  2:1439841/150‑1 (MQ=255)
                                                                                                                           gcATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTtca                          <  2:881421/150‑1 (MQ=255)
                                                                                                                            cATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTtcat                         >  2:645681/1‑150 (MQ=255)
                                                                                                                             aTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCAGCTtcatc                        <  1:1170646/150‑1 (MQ=255)
                                                                                                                              tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTtcatca                       <  1:609031/150‑1 (MQ=255)
                                                                                                                              tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTtcatca                       <  2:743808/150‑1 (MQ=255)
                                                                                                                                tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCAcc                     >  2:1033730/1‑150 (MQ=255)
                                                                                                                                tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCAcc                     >  1:1883762/1‑150 (MQ=255)
                                                                                                                                      ccGGATAAGGCGTTTACGCCGAATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGGCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCaaa               <  1:799482/150‑1 (MQ=255)
                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAAc              <  1:1374516/150‑1 (MQ=255)
                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAAc              <  1:677885/150‑1 (MQ=255)
                                                                                                                                        ggatgcggCGTAAACGCCTTATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCCCTTGCCCCTGCGTCATCCCCAGCAAACg             >  2:337012/7‑150 (MQ=255)
                                                                                                                                        ggATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCCGCAAACg             >  1:394411/1‑150 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACgg            <  2:207411/150‑1 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCACAGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACgg            >  1:1293702/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:1184198/1‑150 (MQ=255)
                                                                                                                                             aGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCg        >  2:1318073/1‑150 (MQ=255)
                                                                                                                                              ggCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGc       >  1:1348697/1‑150 (MQ=255)
                                                                                                                                                 gTTTACGCCGCATCCGGCACCGGCTCTATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCcag    <  1:895295/150‑1 (MQ=255)
                                                                                                                                                   ttACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGTTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCGGCTTCATCAGCTGCAAGCGGTCCGCcagcg  >  1:305007/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCACGCCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCA  >  NC_000913/4317837‑4318133

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