Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,604 G→A 5.7% pseudogene (1213/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,6040GA 5.7% 597.9 / 15.2 193pseudogene (1213/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (114/65);  new base A (5/6);  total (120/73)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

TGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAT  >  NC_000913/1466455‑1466753
                                                                                                                                                     |                                                                                                                                                     
      tCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTAGACGGCGACCAGGCCATGGTTAAAAACGAAGGCGAGAGCGCCATCa                                                                                                                                                 >  1:1290931/1‑150 (MQ=11)
      tCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCa                                                                                                                                                 >  1:1973845/1‑150 (MQ=255)
        aGCGGCGGCGGTCACGGTATTTATAACAACGTCGACAGCGAAACGCAGTATAAAAAGGGCACCATGACAGTCACCGATCCGGAGTCCATCGGTATCCAGGTGGACGGCGACCAGGCGGTTGTTGACAACGAAGGCGAGTGCGCCCTCAcc                                                                                                                                               >  1:5961833/1‑150 (MQ=17)
        aGCGGCGGCGGTCACGGTATTTATAAAACAGGAGACAGCGCAACGGAGGATAACAACGGCACCATTAACGTCACCGATCCGCAGACCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAcc                                                                                                                                               >  1:3166635/1‑150 (MQ=17)
                           ttGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGACCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGc                                                                                                                            >  1:5457811/1‑150 (MQ=255)
                           ttGATATCACCGGCGACAGCGCTACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGc                                                                                                                            <  2:2684097/150‑1 (MQ=255)
                             gATATCACCGGCGACAGCGCTACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAg                                                                                                                          <  1:5527154/150‑1 (MQ=255)
                             gATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGTACGCAg                                                                                                                          <  2:2111070/150‑1 (MQ=255)
                              atatCACCGGCGACAGCTCATCTGTTGATAACTAGGTCACCTTTAACGTCTCAGATCATGAGTCCTTCGGTATCCATGTTGACTGCGACCTGGCGTTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                                                                                                                         <  2:1123260/150‑1 (MQ=18)
                              atatCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTCACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                                                                                                                         >  1:5418173/1‑150 (MQ=255)
                               tatCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAt                                                                                                                        >  1:4578501/1‑150 (MQ=255)
                                atCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAtt                                                                                                                       >  2:6555508/1‑150 (MQ=255)
                                  cACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTaa                                                                                                                     >  2:3312788/1‑150 (MQ=255)
                                  cACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTaa                                                                                                                     >  1:4228866/1‑150 (MQ=255)
                                   cccGGCGACCGCGCAACGGTGGATTACAAGGGCCCCAAGCCCGTCACCGAGCCGGAGTCCATCGGTTTCCAGGCTGTCGGCGCCCCGGCGGTCGTCAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                                                                    <  2:726801/149‑1 (MQ=17)
                                   aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                                                                    <  2:3556722/150‑1 (MQ=255)
                                    ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCCTGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGGTGACGGCGACCAGGCGGGTGTTAACCACGAGGGCGAGAGCGCCAGCCCCACCGGGGGCACCGGCACGCGGATTAACg                                                                                                                   >  2:1734954/1‑150 (MQ=17)
                                    ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                                                                   <  2:826266/150‑1 (MQ=255)
                                    ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                                                                   >  2:4788855/1‑150 (MQ=255)
                                          aCAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgatg                                                                                                             >  1:1410648/1‑150 (MQ=255)
                                                     gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGGCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCa                                                                                                  >  1:1736692/1‑150 (MQ=255)
                                                      tGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCTGATTAACGGTGGTGACGCCGCGGCaa                                                                                                 >  2:2457098/1‑150 (MQ=255)
                                                      tGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaa                                                                                                 >  1:428380/1‑150 (MQ=255)
                                                       ggATAACCAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaaa                                                                                                <  1:1100806/150‑1 (MQ=255)
                                                       ggATAACAAGGGCACCAGGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCTACGTTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaaa                                                                                                <  1:300224/150‑1 (MQ=255)
                                                        gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCTCCATCACCAACGGAGGCACCGGCACGCAGATAAACGGTGATGACGCCACCGCaaaa                                                                                               >  1:3160097/1‑149 (MQ=255)
                                                          tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaa                                                                                             >  2:5329513/1‑150 (MQ=255)
                                                          tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaa                                                                                             >  2:3962697/1‑150 (MQ=255)
                                                           aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaac                                                                                            <  2:1255566/150‑1 (MQ=255)
                                                            aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCAaaaaaca                                                                                           >  1:491728/1‑150 (MQ=17)
                                                                    aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGTTGATGACGCCACGGCAAACAACAACGGCaaa                                                                                   <  2:3458309/150‑1 (MQ=255)
                                                                    aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaaa                                                                                   <  2:108245/150‑1 (MQ=255)
                                                                      cATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                 >  1:3542577/1‑150 (MQ=255)
                                                                      cATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                 >  1:6472180/1‑150 (MQ=255)
                                                                       aTGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGGGGCACCGGCACTCAGATAAACGGTAACGACGCCACCGCGAAGAACAGAGGCAAAacc                                                                                >  2:3134979/1‑150 (MQ=11)
                                                                             gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                          <  2:5490777/150‑1 (MQ=255)
                                                                             gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                          <  2:6227829/150‑1 (MQ=255)
                                                                                aCCGATCCGGAGTCCATCGGTATCCAGGTCGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc                                                                       <  2:4880992/150‑1 (MQ=255)
                                                                                  cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCATGCGGTTGTTAACAACGAAGTCGAGAGCGCAATCACCAAAGTAGGCACCGGCACTCAGATAAACGGTGAAGACGCCACCGCAAAAAACAAAGGAAAAACCACAGATGACgg                                                                     >  1:5056488/1‑150 (MQ=17)
                                                                                  cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                     >  1:5160020/1‑150 (MQ=255)
                                                                                  cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                     >  1:4019679/1‑150 (MQ=255)
                                                                                    aTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCa                                                                   >  2:2064353/1‑150 (MQ=255)
                                                                                      ccGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                 <  2:3878580/150‑1 (MQ=255)
                                                                                       cGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                >  2:5772308/1‑150 (MQ=255)
                                                                                        ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGa                                                               >  1:5435828/1‑150 (MQ=255)
                                                                                        ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACTCAGATAAACGGTGAAGAAACAACGTCAAGAAACAACGGAAAAACCACATTAGACTTCAAGGa                                                               >  2:1962777/1‑150 (MQ=18)
                                                                                             ccATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                          <  2:4621245/150‑1 (MQ=255)
                                                                                               aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                        >  1:5689645/1‑150 (MQ=255)
                                                                                                tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACAGCAAGGATTCCACCg                                                       >  1:2203538/1‑150 (MQ=255)
                                                                                                   gTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                    >  1:3849941/1‑150 (MQ=255)
                                                                                                   gTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAAAAACAAAGGCAAAACCACCGTTGACGGCAAGGATACCACCGGTa                                                    >  1:3971592/1‑150 (MQ=255)
                                                                                                      tCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                 >  1:3930191/1‑150 (MQ=255)
                                                                                                         aGGTTGACGGCGACCAGGCGGTTGTTAACAACGATGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                              >  1:4089380/1‑150 (MQ=255)
                                                                                                         aGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                              <  2:3106290/150‑1 (MQ=255)
                                                                                                         aGGTAGACGGCGACCAGGCGTTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                              <  2:6206006/150‑1 (MQ=255)
                                                                                                             tGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGc                                          <  2:4205146/150‑1 (MQ=255)
                                                                                                             tGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGc                                          <  2:3111202/150‑1 (MQ=255)
                                                                                                             tGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGc                                          <  2:2442309/150‑1 (MQ=255)
                                                                                                              gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                         <  1:4681154/150‑1 (MQ=255)
                                                                                                                 ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCAATCACCAAAGGTGGCACCGGCACGCAGATTAACGGTGAGGACGCCACCGCAAAAAACAAAGGCAAAACCACAGTTGAAGGAAAGGATACCACCGTAACGGAAAGTGCTGGc                                      >  1:1624781/1‑150 (MQ=255)
                                                                                                                 ggCGACCAGGCGGTCGTCAAGAACGACGGCGAGAGCGCCATCACCACCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAACCAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGc                                      <  2:6589833/150‑1 (MQ=14)
                                                                                                                       cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAAc                                >  1:3623690/1‑150 (MQ=255)
                                                                                                                       cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAAc                                >  1:2517180/1‑150 (MQ=255)
                                                                                                                        aGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACg                               >  2:2465488/1‑150 (MQ=255)
                                                                                                                         ggCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                              >  1:6394736/1‑150 (MQ=255)
                                                                                                                              ttgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAgg                         >  1:2283728/1‑150 (MQ=255)
                                                                                                                               tgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCCAACAACAACGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACCCAAATCGCTGGCAATAACGGCCAGGt                        >  2:4005845/1‑150 (MQ=17)
                                                                                                                                gttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTg                       >  1:4342893/1‑150 (MQ=255)
                                                                                                                                gttAACAACGAAGGCGAGAGCGCAATCACCAACGGGGGCACCGGCACGCAGATAAACGGTGAGGACGCCACCGCGAAAAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACCGAAATAGCGGGCAATAACGGGAAGGTg                       >  2:3671036/1‑150 (MQ=17)
                                                                                                                                gttAACAACCAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGATGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTg                       >  1:5853682/1‑150 (MQ=255)
                                                                                                                                 ttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACTGGAAGGTGa                      >  1:3337096/1‑150 (MQ=255)
                                                                                                                                 ttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGa                      >  2:552907/1‑150 (MQ=255)
                                                                                                                                  tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAt                     >  1:537188/1‑150 (MQ=255)
                                                                                                                                   aataacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                              >  1:480152/4‑140 (MQ=255)
                                                                                                                                   aataacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                              >  1:2892717/4‑140 (MQ=255)
                                                                                                                                   aataacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                              <  2:2892717/137‑1 (MQ=255)
                                                                                                                                   aacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAtt                    <  1:4971356/150‑1 (MQ=255)
                                                                                                                                      aacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAAGACGCCACGGCGAACAACAAAGGCAAAACCAACGCTGACGGCAAGGATTCCACAGCTCCGGAAATTGCTGGCAAGAACGGGAAGGTGATTCAg                 >  1:1890700/1‑150 (MQ=17)
                                                                                                                                       acGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAgg                <  1:6153023/150‑1 (MQ=255)
                                                                                                                                        cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCTAACAACAAAGGCAAAACCACCGTTGAAGGCAAGGATTCCACCGGTACGGAAATTGCTGGCTGTAACGGGAAGGTGATTCAGGa               >  1:73135/1‑150 (MQ=17)
                                                                                                                                        cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa               <  2:5740430/150‑1 (MQ=255)
                                                                                                                                        cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa               <  2:2607787/150‑1 (MQ=255)
                                                                                                                                        cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa               >  2:4486808/1‑150 (MQ=255)
                                                                                                                                         gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGAGGGGAATAACGGGAAGGTGATTCAGGAc              >  2:6060492/1‑150 (MQ=11)
                                                                                                                                         gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAc              <  2:5899908/150‑1 (MQ=255)
                                                                                                                                         gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGAGGGCAATAACGGGAAGGTGATTCAGGAc              >  2:4802482/1‑150 (MQ=11)
                                                                                                                                          aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACg             <  2:3036460/150‑1 (MQ=255)
                                                                                                                                             gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGACGGTGATTCAGGACGGCg          >  1:84288/1‑150 (MQ=255)
                                                                                                                                             gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg          >  1:5004704/1‑150 (MQ=255)
                                                                                                                                               gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt        >  2:6093203/1‑150 (MQ=255)
                                                                                                                                               gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt        >  1:4643953/1‑150 (MQ=255)
                                                                                                                                               gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGATGACGGCAAGGAATCCACCGGTACGGAACTTGCTGGCAATAACGGGCAGGTGAGTCAGGACGGCGAt        >  2:6460989/1‑150 (MQ=17)
                                                                                                                                               gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGCGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATCAGGTGAAGGTGATTCAGGACGGCGAt        >  2:1309824/1‑150 (MQ=255)
                                                                                                                                                agagCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATc       >  1:2067820/1‑150 (MQ=255)
                                                                                                                                                 gagCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCt      <  1:4747705/150‑1 (MQ=255)
                                                                                                                                                 gagCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCt      <  2:446644/150‑1 (MQ=255)
                                                                                                                                                  agCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTg     <  2:3711660/150‑1 (MQ=255)
                                                                                                                                                  agCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTg     >  2:6123178/1‑150 (MQ=255)
                                                                                                                                                    cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGa   <  1:2793742/150‑1 (MQ=255)
                                                                                                                                                     gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt  >  1:2991278/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAT  >  NC_000913/1466455‑1466753

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: