Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,325,949 A→G 6.0% intergenic (‑207/+451) 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,325,9490AG 6.0% 254.4 / 3.7 107intergenic (‑207/+451)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base A (41/57);  new base G (2/6);  total (43/64)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGG  >  NC_000913/4325801‑4326096
                                                                                                                                                    |                                                                                                                                                   
gagGCGGCCTGATAAGACGCGCGCGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCa                                                                                                                                                     <  1:891915/149‑1 (MQ=255)
 agGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCAt                                                                                                                                                    <  1:1155727/149‑1 (MQ=40)
    cggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGc                                                                                                                                                 <  1:1297207/149‑1 (MQ=39)
    cggcCCGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGc                                                                                                                                                 <  1:1297290/149‑1 (MQ=38)
     tgcCTGAGAAGACGCGCGCGCGTCGCCTCAGGCAGTCGGCCTTCTAGCCTGATGCGGCGTAACCGCCTTATCCGGCCGCCGGAGGGTGCGGGATCTTGTAGGCCTGACAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGcc                                                                                                                                                <  2:164422/148‑1 (MQ=11)
        cTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGa                                                                                                                                             <  2:1167817/149‑1 (MQ=38)
             acggcgcgcgAGCGTCGCATCAGGCAGTCGGCGCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcgg                                                                                                                                        <  2:657103/145‑1 (MQ=255)
             aaGACGCGCGGGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcgg                                                                                                                                        <  2:657382/149‑1 (MQ=38)
              gcccgtgcgAGCGTCGCATCCGGCAGTCGGCACGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGTCACTTTGTCGGCCTGATAAGACGCGCAAGCGGGGGAGCAGGCACTCGGCACCATTGCCGGATgcggc                                                                                                                                       <  2:14189/143‑1 (MQ=11)
              aGACGCGCGGGCGTCGCAGCAGGCGGTCGGCACTGTTGCGGGATGCGGCGTATGCGTCTTATGCGGCGTGCGGAGTGTGGGGGGATTTGTAGGCCTGCTAAGACGCGCAAGCGTCGCCTCAGGCAGTCGGCACCATTGCCGGATggggc                                                                                                                                       <  2:685814/149‑1 (MQ=11)
              aGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcggc                                                                                                                                       >  1:661870/1‑149 (MQ=38)
              aGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcggc                                                                                                                                       <  1:454650/149‑1 (MQ=38)
               gACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcggcg                                                                                                                                      >  1:210804/1‑149 (MQ=17)
               gACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcggcg                                                                                                                                      <  1:410109/149‑1 (MQ=38)
                aCGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGt                                                                                                                                     >  1:1045663/1‑149 (MQ=17)
                    gcgAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAAc                                                                                                                                 <  1:1307287/149‑1 (MQ=38)
                       aGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGcc                                                                                                                              <  1:446686/149‑1 (MQ=35)
                         cGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCtt                                                                                                                            <  1:1078815/149‑1 (MQ=35)
                         cGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCtt                                                                                                                            <  2:1145469/149‑1 (MQ=35)
                             gCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                        >  1:1167122/1‑149 (MQ=35)
                               aTCGGGCAGTCGGCACTGTTGCCGGGTGCGGCGTAAACGCCTTGTCCGGCCTACGGAGCGTGCTTGAAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                      <  1:700488/149‑1 (MQ=11)
                               aTCAGGCAGTCGGCACTGTTGCCGGATGCGGAGCGAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                      >  1:712369/1‑149 (MQ=255)
                                    gCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTGCGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                 >  2:734219/1‑149 (MQ=25)
                                    gCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                 <  1:289220/149‑1 (MQ=35)
                                     cAGTCGGCACTGTTTCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGTAGGGTGCGGGAATTTGTAGGCCTGAGAAGACGCGCAAGCGTCGAATCAGGCTGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg                                                                                                                >  1:1074333/1‑149 (MQ=11)
                                       gTCGGCATTGTGGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGAGGGTGCGGGAATTTGTCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                              <  2:771963/149‑1 (MQ=255)
                                            cgcTGTGGCCGGATGCGGCGTAAACGCCATATGCGGCCTACAGAGGGTGCGGGAATGTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTg                                                                                                         <  2:566731/147‑1 (MQ=11)
                                            cccTGGTGCCGGATGCGGCGGGACCGCGTTGTCCGGCCTACGGGGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTg                                                                                                         <  2:260451/147‑1 (MQ=11)
                                              cTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                       <  1:1040817/149‑1 (MQ=32)
                                              cTGTCGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGGAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                       <  2:389682/149‑1 (MQ=17)
                                               tgctGCCGGATGCGGCGTAAGCGCCTTATCCGGCCCGCGGGTGGTGCGGGCATTTGTCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg                                                                                                      <  2:1393883/146‑1 (MQ=11)
                                               tGTTGCCGGATGCGGCGTAAACGCCTTATCCGGACTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGc                                                                                                                               >  1:1063459/1‑124 (MQ=17)
                                               ggTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg                                                                                                      <  1:159390/148‑1 (MQ=17)
                                                gTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCggg                                                                                                     >  2:1339733/1‑149 (MQ=32)
                                                   gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcggcg                                                                                                                                      >  1:866361/1‑113 (MQ=32)
                                                     cGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCTTCATTGCCGGTTTCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAttt                                                                                                <  1:347277/149‑1 (MQ=255)
                                                           cggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc                                                                                          <  1:669196/149‑1 (MQ=17)
                                                             gcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCt                                                                                        >  1:263193/1‑149 (MQ=32)
                                                             gcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGGAGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCt                                                                                        <  2:594810/149‑1 (MQ=17)
                                                               gaacaCGTCTAATCCGGCCTGCGGAGGGTGCGGGAATGTGGAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                                                                                      <  2:19066/145‑1 (MQ=11)
                                                                   aCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCACCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGAGCGGGGATTTGTAGGCCTGATGAg                                                                                  >  2:1388442/1‑149 (MQ=255)
                                                                     gCCTTATCCGGCCTACGGAGGGTGCGGGAAAATGCTGTCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAc                                                                                >  1:1441525/1‑149 (MQ=11)
                                                                       cTTGTCCGGGCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCACATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                              <  2:711103/149‑1 (MQ=255)
                                                                             cGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGc                                                                        <  2:306331/149‑1 (MQ=32)
                                                                             cGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGc                                                                        >  1:1179941/1‑149 (MQ=32)
                                                                              ggCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg                                                                       >  2:469243/1‑149 (MQ=32)
                                                                               gCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGt                                                                      <  1:363784/149‑1 (MQ=32)
                                                                               gCCTACAGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAGTCGGCCCCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCAACGGAGGGTGCGGGGATTTGTAGGGCTGAGAAAGCGCGCAAGCGt                                                                      >  2:1184505/1‑149 (MQ=11)
                                                                                ccTACGGAGGGCGCGGTAATTTGTAGGCCTGATAAGCCCCGCACGCGTCCGATCAGTCATTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTCCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                                                     <  2:1090122/149‑1 (MQ=11)
                                                                                  tACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGGATTGGTAGGCCTGCTGAGACGGGCAAGGGTGGc                                                                   >  2:692504/1‑149 (MQ=11)
                                                                                   aCGGTGGGTGCGGGAGTGGGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGACTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCa                                                                  <  2:713087/149‑1 (MQ=11)
                                                                                     ggAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                                >  1:255211/1‑149 (MQ=32)
                                                                                         ggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGc                                                            <  1:792546/149‑1 (MQ=32)
                                                                                             cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGcc                                                                                                                                                >  1:1102802/1‑61 (MQ=35)
                                                                                             cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAttt                                                                                                >  1:729000/1‑109 (MQ=33)
                                                                                             cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTc                                                        >  1:908625/1‑149 (MQ=34)
                                                                                              gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCt                                                                                        >  2:1124576/1‑116 (MQ=32)
                                                                                              ggaggtttGGAGGCCTGAGAGGCCGCGCAAGCGTGGCATCAGGTAGTCCGCACCTTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTAGTATGCCTGAGAAGACGCGCGAGCGTCGCATCAGGCAGTCg                                                       <  2:892854/144‑1 (MQ=11)
                                                                                               ggAATTTGTAGGCCTGATAAGACCCGCAAGCGTCGAATCAGGCAGGCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCgg                                                      <  1:6197/149‑1 (MQ=255)
                                                                                                  aTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTCCACGCCTTATCCGGCATATTTAGGGTGCGGGAATTTGTAGGCCTGATCAGACGCGCAAGCGCTGCATCAGGCAGTCGGCAc                                                   <  1:1066450/149‑1 (MQ=11)
                                                                                                  aTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCAc                                                   >  1:851497/1‑149 (MQ=38)
                                                                                                   tgtGTAGGCCTGATAAGACGCGCAAGCCTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCAcc                                                  <  1:829330/147‑1 (MQ=255)
                                                                                                    tggTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCg                                                 <  2:558628/147‑1 (MQ=17)
                                                                                                      ggaGGCCTGCTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATGTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGtt                                               <  1:435139/147‑1 (MQ=255)
                                                                                                      gTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGtt                                               <  1:401782/149‑1 (MQ=37)
                                                                                                          gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCGTCAGGCAGTCGGCGCCGTTGCCg                                           >  1:1109536/1‑149 (MQ=255)
                                                                                                          gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGtt                                                        >  2:641955/1‑135 (MQ=18)
                                                                                                           ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                    <  1:37847/139‑1 (MQ=14)
                                                                                                           ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                    >  2:37847/1‑139 (MQ=14)
                                                                                                           ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCgg                                          <  1:857030/149‑1 (MQ=38)
                                                                                                            cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                                                                             >  2:908808/1‑97 (MQ=33)
                                                                                                             tGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCTAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGAt                                        <  2:908625/149‑1 (MQ=14)
                                                                                                             tGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                     >  2:1008065/1‑72 (MQ=34)
                                                                                                             tGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCGCCGTTGCCGGAt                                        >  1:137030/1‑149 (MQ=255)
                                                                                                              gATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATg                                       <  1:462905/149‑1 (MQ=17)
                                                                                                               aTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgc                                      <  2:743116/149‑1 (MQ=37)
                                                                                                                 aaGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGCTgcgg                                    >  1:448909/1‑149 (MQ=255)
                                                                                                                   gACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAGTCCGCACCGTTGCCGGATgcggcg                                  >  1:187447/1‑149 (MQ=255)
                                                                                                                    aCGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGt                                 <  1:734219/149‑1 (MQ=38)
                                                                                                                           aGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGcc                          <  2:506106/149‑1 (MQ=38)
                                                                                                                           aGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGcc                          <  1:1339733/149‑1 (MQ=39)
                                                                                                                             cGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGTCCGAACAACAGAGCTGTAATGTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGGCGGGTGCGGCGGAAACGCCtt                        >  2:1251685/1‑149 (MQ=11)
                                                                                                                              gTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAACTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTCGCCGGATGCGGCGCAAACGccaaa                       >  2:910617/1‑146 (MQ=255)
                                                                                                                                  cagcAGCCAGTCGGCCCCGTTGCCGGATGCGGCGTAACCGCCTTGTCCGGCCCCAGGAGGTTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGGCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCg                   <  2:456636/146‑1 (MQ=11)
                                                                                                                                  cATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCg                   <  1:1452189/149‑1 (MQ=38)
                                                                                                                                  cATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCAGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAGAAGAGCCGCAAGCGGGGCATGAGGCAATGGGCACCGTTAGCGCGCGCGGCGCAAAGGACCTATCCg                   >  2:914992/1‑149 (MQ=11)
                                                                                                                                   aTCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                  <  1:469243/149‑1 (MQ=38)
                                                                                                                                    tCAGGCAGTCGGCACCGTGGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                 <  2:281286/149‑1 (MQ=17)
                                                                                                                                    tCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtt                                                                                                 >  1:1412780/1‑69 (MQ=37)
                                                                                                                                    tCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtt                                                                                                 <  2:1412780/69‑1 (MQ=35)
                                                                                                                                        gCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGc                                             <  1:1083191/117‑1 (MQ=39)
                                                                                                                                        gCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGc                                             >  2:1083191/1‑117 (MQ=255)
                                                                                                                                          aGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCg                                           <  2:117187/117‑1 (MQ=17)
                                                                                                                                          aGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCg                                           >  1:117187/1‑117 (MQ=40)
                                                                                                                                              ggCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGc                                             >  1:748188/1‑111 (MQ=37)
                                                                                                                                                 aCCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                                                                                      >  2:394361/1‑67 (MQ=35)
                                                                                                                                                  ccATTGCCGGATGCGGCGTAAACGCCTTACCCGGCCGACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCg   <  1:630034/149‑1 (MQ=25)
                                                                                                                                                   cATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGAt                                        <  2:537222/111‑1 (MQ=31)
                                                                                                                                                   cATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGAt                                        >  1:537222/1‑111 (MQ=31)
                                                                                                                                                   cATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTGCCCGGAGGCGGCGGAAACGCCTTATCCGGCCTACGGATCGGGCgg  >  2:1284451/1‑149 (MQ=11)
                                                                                                                                                    |                                                                                                                                                   
GAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGG  >  NC_000913/4325801‑4326096

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: