Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,242 C→T 5.8% intergenic (‑500/+158) yjdN ← / ← yjdM metalloprotein superfamily protein/zinc‑ribbon family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,326,2420CT 5.8% 225.3 / 4.5 100intergenic (‑500/+158)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base C (50/44);  new base T (2/4);  total (52/48)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

Only 100 of 101 total aligned reads displayed.

CGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGAT  >  NC_000913/4326094‑4326389
                                                                                                                                                    |                                                                                                                                                   
cGGGCATTTGTAGGCCTGATAAGAATCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGTCGTAAACGCCTTATCCGGCCTACGCAGAGTGCAGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTc                                                                                                                                                     >  1:1000764/1‑149 (MQ=11)
  ggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACCGTTGCCGGATGCGGCGCAAACGCCTTATCCGGCCTACGGAGGGCGCGGGGAAAGGGCGGCCTGATAAGAAGCGGAAGCGGCCCGACAGGCAGCCgg                                                                                                                                                   >  2:145190/1‑149 (MQ=11)
   gAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTGGGCCGGAGAAGACGCGGAAGCGTCGGATCAGGCAGCCggg                                                                                                                                                  >  2:801677/1‑148 (MQ=11)
        tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGt                                                                                                                                             <  2:399482/149‑1 (MQ=37)
          tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCCCCGTTg                                                                                                                                           >  1:420214/1‑149 (MQ=255)
            ggCCTGATAAGGCGCTCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGcc                                                                                                                                         <  1:1004660/149‑1 (MQ=255)
               cTGATAAGACGCGGAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGa                                                                                                                                      <  1:696196/149‑1 (MQ=38)
               cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATGAGACGCGCAAGCGTCGCATCAGGGAGCCGGCACCGTTGCCGGa                                                                                                                                      >  2:1018823/1‑149 (MQ=255)
                 gATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGCCGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGCTg                                                                                                                                    >  1:680024/1‑149 (MQ=255)
                     aGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcggc                                                                                                                                >  1:783528/1‑149 (MQ=39)
                       ccgcgCACGCGGCGCAGCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGt                                                                                                                              <  2:855145/148‑1 (MQ=255)
                         gcgcAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGtga                                                                                                                            >  2:858248/1‑147 (MQ=38)
                          cgcAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTaaa                                                                                                                           <  1:1296442/149‑1 (MQ=38)
                          cgcAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCAGCACCGTTGCCGGATGCGGCGTgaa                                                                                                                           >  1:829585/1‑149 (MQ=255)
                            cAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACg                                                                                                                         >  1:554915/1‑149 (MQ=37)
                            cAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACg                                                                                                                         >  1:38100/1‑149 (MQ=38)
                                      aTCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCGAGCGTCGCCTCAGGCAGTCGGCACCGTTGCCGGATGCGGCGACCGCGCCTTAGCCgt                                                                                                               >  2:1149108/1‑148 (MQ=11)
                                       tCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTTCGGGAATGTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                              >  1:523902/1‑149 (MQ=25)
                                       tCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGGAAACGCCTTATCCGGc                                                                                                              >  2:717602/1‑149 (MQ=25)
                                          ggCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCCCCTGGGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                           <  1:1337228/149‑1 (MQ=255)
                                          ggCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGTGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGGCAACGGCTTATCCGGccga                                                                                                           >  2:1158145/1‑147 (MQ=11)
                                                  gCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGt                                                                                                   >  2:608113/1‑149 (MQ=32)
                                                      tGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGTTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg                                                                                               >  1:391971/1‑149 (MQ=17)
                                                      tGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg                                                                                               <  2:38100/149‑1 (MQ=32)
                                                              gATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTgg                                                                                       >  2:820643/1‑148 (MQ=17)
                                                               aTGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATGCGTCGCCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGCGCGGGAATTTGTa                                                                                      <  1:924917/149‑1 (MQ=11)
                                                               aTGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                                                                      <  1:358195/149‑1 (MQ=32)
                                                                 gcggcgTAACCGCCTAGTGCGGCTTGCGGAGGGTGCGGGAATGTGGAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAgg                                                                                    <  2:652848/149‑1 (MQ=11)
                                                                  cggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGCCGGCACCGTTCCCGGGAGCGGCGTAAACGCCTTAACCGGCCTACGGAGGGTCCGGGCATGTGGAGGc                                                                                   >  2:855349/1‑149 (MQ=11)
                                                                   ggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGcc                                                                                  >  1:697070/1‑149 (MQ=17)
                                                                    gcgTTATCGCCTTACCCGGCCTACGGGGGGTGCGGGAATTTGTAGGCCTGATACGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTAGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTTCGGGAAGTTGTAGGCCt                                                                                 <  2:1000764/149‑1 (MQ=11)
                                                                    gcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAt                                                                                           <  2:255086/139‑1 (MQ=32)
                                                                        aaaCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATa                                                                             <  2:539424/149‑1 (MQ=32)
                                                                          aCGCCTTATCCGGCCTACGGAGGGCGCGGGGATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGGCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCGACGGAGGGTGCGGGGATTTGTTGGCCTGATAAg                                                                           >  2:41439/1‑149 (MQ=11)
                                                                              cTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCGGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACCGGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                       <  2:131409/149‑1 (MQ=12)
                                                                               ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                                                                       >  2:1324140/1‑132 (MQ=32)
                                                                               ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                                                                       <  1:1324140/132‑1 (MQ=32)
                                                                                tATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGGCGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACGGAGGGCGCGGCAATTTGTAGGCCTGCTAAGACgcgc                                                                     >  2:1013066/1‑149 (MQ=11)
                                                                                 aTCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCa                                                                    >  1:792256/1‑149 (MQ=32)
                                                                                     tgCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg                                                                <  2:1267709/148‑1 (MQ=17)
                                                                                     ggCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg                                                                >  2:1085279/1‑149 (MQ=17)
                                                                                      gCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATaa                                                                            <  1:268171/136‑1 (MQ=33)
                                                                                      gCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATaa                                                                            >  2:268171/1‑136 (MQ=32)
                                                                                      gCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGCAAACGCCTTATCCGGCCTACCGAGGGCGCGGGAATCTGTAGGCCTTCTAAGAcgcggatccgt                                                               >  2:522474/1‑142 (MQ=11)
                                                                                       ccTGCGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCACCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                                              <  2:379672/149‑1 (MQ=17)
                                                                                       ccTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                                                                      <  2:1092479/125‑1 (MQ=32)
                                                                                        ctgcGGAGGGTGCGGGAATTGGTGGGGCGCATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCg                                                             <  2:1039443/146‑1 (MQ=11)
                                                                                        cTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg                                                                                               <  1:306079/115‑1 (MQ=33)
                                                                                          atgtaGGCTGGGGGAATTTGTAGGTCTGATAAGACGCGCAAGCGTCTCATCAGGCAGTCGGCACCGTTGCCGGATGCTGCGTGCACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCa                                                           <  1:1297359/145‑1 (MQ=11)
                                                                                           cGGATGGTGCGGGAATTTGTAGGCCTGATAAGGCGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCAt                                                          <  2:709159/149‑1 (MQ=12)
                                                                                               gggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAg                                                                                     <  1:1115737/118‑1 (MQ=33)
                                                                                               gggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAg                                                                                     >  2:1115737/1‑118 (MQ=32)
                                                                                                   gCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATc                                                  >  1:382548/1‑149 (MQ=34)
                                                                                                    cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCt                                                 <  1:448307/149‑1 (MQ=17)
                                                                                                    cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCCTCAGGCATCt                                                 >  1:1133991/1‑149 (MQ=21)
                                                                                                    cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCt                                                 >  1:1346978/1‑149 (MQ=35)
                                                                                                    cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAGTTGGTAGGCCTGATAAGACACGCAAGCGGATCATGAGGCATCt                                                 >  1:685277/1‑149 (MQ=11)
                                                                                                     gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGc                                                            <  1:735137/137‑1 (MQ=32)
                                                                                                       gAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                                              <  1:448655/117‑1 (MQ=32)
                                                                                                        aaTTTGTAGGCCTGATAAGACGTGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCTGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGCCGCGCAAGCGTCTCATCAGGCATCTGGCa                                             <  1:855349/149‑1 (MQ=255)
                                                                                                        aaTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCa                                             <  2:420214/149‑1 (MQ=38)
                                                                                                        aaTTTGTAGGCCTGATAAGACGCGCAAGCAACACATCCGGTAGTCGTCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCa                                             <  1:1043388/149‑1 (MQ=11)
                                                                                                         aTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATTTGGCat                                            >  1:1229626/1‑148 (MQ=17)
                                                                                                         aTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCAc                                            >  1:572956/1‑149 (MQ=38)
                                                                                                         aTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGGCGGCACCGTTGCCGGATGCGGCGTAAACGCCGTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCAc                                            <  1:154233/149‑1 (MQ=255)
                                                                                                          tttGTAGGCCGAACAAGACGCGCCAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCAcc                                           <  1:6224/149‑1 (MQ=255)
                                                                                                           ttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGCAAGGGGCTTATCGGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTGGGATTCGGCATCTGGCACCg                                          >  2:1158300/1‑149 (MQ=11)
                                                                                                            tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGt                                         <  1:695716/149‑1 (MQ=40)
                                                                                                              tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGGATTTGTAGGCAGGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTg                                       >  1:591657/1‑149 (MQ=255)
                                                                                                              tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTCGGCCTGATAAGACGCGCCAGCGTCGCATCGGGCATCTGCCCCCGtgg                                       >  2:767787/1‑147 (MQ=11)
                                                                                                               aGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGc                                      >  1:1349497/1‑149 (MQ=255)
                                                                                                                 gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCg                                    >  2:937858/1‑149 (MQ=38)
                                                                                                                 gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCg                                    <  2:146777/149‑1 (MQ=255)
                                                                                                                 gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGGATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGCTGGCACCGTTGCCg                                    >  2:342514/1‑149 (MQ=25)
                                                                                                                 gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCt                                                                                 <  2:1163325/104‑1 (MQ=32)
                                                                                                                 gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCg                                    >  1:348695/1‑149 (MQ=39)
                                                                                                                  ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAg                                                                                                      <  2:445940/82‑1 (MQ=31)
                                                                                                                   cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGc                                                            <  1:85506/123‑1 (MQ=32)
                                                                                                                    tGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                                              <  2:278690/104‑1 (MQ=32)
                                                                                                                    tGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGAt                                 >  1:1310698/1‑149 (MQ=39)
                                                                                                                     gttacggcgcgCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATg                                <  1:1066756/142‑1 (MQ=255)
                                                                                                                     gATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATg                                <  1:841297/149‑1 (MQ=40)
                                                                                                                     gATAAGACGCGCAAGCGTCGCATCAGGCAGCCGGACCCGTTGCCGGAGGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATg                                <  1:1116438/149‑1 (MQ=255)
                                                                                                                      aTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATgc                               >  1:740866/1‑149 (MQ=40)
                                                                                                                        aaGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATgcgg                             <  1:1285756/149‑1 (MQ=40)
                                                                                                                         aGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATgccgc                            >  1:72804/1‑149 (MQ=255)
                                                                                                                           aCGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCCAGCATCTGGCACCCTTGCCGGAtactgggc                          >  2:360242/1‑142 (MQ=14)
                                                                                                                                cAAGCGTCGCATCAGACAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGa                        <  1:783301/146‑1 (MQ=255)
                                                                                                                                cAAGCGTCGCATCAGACAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGa                        >  2:783301/1‑146 (MQ=255)
                                                                                                                                 aaGCGTCGCATCAGGCAGTCGGCTCCGTTGCGGGATGTGGTGTAAAGGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATTTGGCACCGATGCCGGATGCGGCGCGAgcgc                    <  1:704715/149‑1 (MQ=14)
                                                                                                                                  aGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCACCTGGCCCCGTTCCCGGATGCGGCGCGAGCGcc                   >  2:246590/1‑149 (MQ=255)
                                                                                                                                     gTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCTAGCGCCTTa                >  1:1168613/1‑149 (MQ=255)
                                                                                                                                       cGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATc              <  1:937858/149‑1 (MQ=255)
                                                                                                                                        gCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATcc             <  1:456089/149‑1 (MQ=255)
                                                                                                                                          aTCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCgg           >  2:1307425/1‑149 (MQ=255)
                                                                                                                                           tCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCCCCGTTGCCGGACGCGGCGCGACCGCCCTATCCGGc          >  2:589354/1‑149 (MQ=14)
                                                                                                                                           tCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                       <  2:1003643/104‑1 (MQ=32)
                                                                                                                                           tCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGc          >  2:915585/1‑149 (MQ=255)
                                                                                                                                              ggCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCAt                                                          <  2:179030/98‑1 (MQ=33)
                                                                                                                                                   tCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGAt  >  1:129464/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGAT  >  NC_000913/4326094‑4326389

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