Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,630 G→C 5.2% pseudogene (1239/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,6300GC 5.2% 273.0 / 6.5 89pseudogene (1239/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (57/25);  new base C (3/2);  total (62/27)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466487‑1466779
                                                                                                                                               |                                                                                                                                                     
atCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAtt                                                                                                                                                 >  2:2683795/1‑150 (MQ=255)
    cgggCGACAGCGGAGCGGTGGATATGTAGGGCACCTTGTCGGTCACGGTTGCGGAGTCCATGGGTATCCAGTTTGACGGCGACCAGGCGGTGGTTAAGAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                                                                                             <  2:1696739/148‑1 (MQ=17)
     cGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACgg                                                                                                                                            >  2:2318074/1‑150 (MQ=255)
      ggCGACAGCGCAACCGTTGATATCAAGGGCACCATGACCGTCTCCGATCCTGTTCGCATCTGTCTACTGTTTGACGGTGACAAGGCCTTAGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGt                                                                                                                                           <  2:982283/150‑1 (MQ=2)
          aCAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTGAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACTCAGATCAACGGtgagg                                                                                                                                       >  2:1524379/1‑148 (MQ=14)
                   cGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGTTTGACGGAGACAAGGCGTTTGTTAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACgg                                                                                                                              >  1:1100829/1‑150 (MQ=11)
                     gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCa                                                                                                                            >  2:1579581/1‑150 (MQ=255)
                     gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCa                                                                                                                            >  2:2883763/1‑150 (MQ=255)
                          tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAGGACGCCACGGCAaacaa                                                                                                                       >  1:2325316/1‑150 (MQ=255)
                          gcgatctggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTTGACGGCGACAAGGCGGTTGTTAAGAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaa                                                                                                                       <  2:2310068/143‑1 (MQ=14)
                             cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTGAACAACGAAGGCGAGAGCGCGATCACCAAAGGGGGCACCGGCACGCAGATGAACGGTGAGGACGCCACGGCGaagaacag                                                                                                                    >  1:1636334/1‑149 (MQ=14)
                              aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCAaacaacaac                                                                                                                   >  1:1692624/1‑150 (MQ=255)
                                     ccATGACCCTCACCGATCCGGAGTCCATCGGTATCCAGCTAGACGGCGACCAGGCGGTAGATAACAACGAAGGCGAGAGCGCCAACACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaaaa                                                                                                            <  2:1100829/150‑1 (MQ=17)
                                       aTGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAGGGAGGCACCGGCACTCAGATAAACGGTGAAGACGCCACCGCGAAAAACAGAGGAAAAacc                                                                                                          >  2:2343124/1‑150 (MQ=14)
                                          aaccacACCGATCATCAGACCATCGGTATCCAGCTAGACGGAGACCAGGCCTTTGTAAAAAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacc                                                                                                       <  2:642647/145‑1 (MQ=14)
                                             gTCACCGATCCGTAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                                                    <  2:1827275/150‑1 (MQ=255)
                                             gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTGAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATAAACGGTGAGGACGCCACCGCCAAAAACACCGGCAAAACCACCGtt                                                                                                    >  1:602579/1‑150 (MQ=17)
                                             gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTGAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCCAAAAACACAGGCAAAACCACAGtt                                                                                                    >  1:322413/1‑150 (MQ=17)
                                              tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCGATCACCAAAGGCGGCACCGGCACGCAGATAAACGGTGAGGACGCCACGGCGAATAACAGAGGAAAAACCACCGTTg                                                                                                   >  2:2645233/1‑150 (MQ=14)
                                              tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                                                   >  2:1036046/1‑150 (MQ=255)
                                               cACCGATCCGGAGTCCATCGGTATCCAGTTTGACGGTGACAAGGCCTTTGTAAAAAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGTTGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                  <  1:1185760/150‑1 (MQ=11)
                                                  cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                               >  2:2657669/1‑150 (MQ=255)
                                                  cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                               >  2:2657670/1‑150 (MQ=255)
                                                   gATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc                                                                                              >  2:2364066/1‑150 (MQ=255)
                                                     tCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGAAGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                            >  1:995459/1‑150 (MQ=255)
                                                       cGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                                          <  2:216096/150‑1 (MQ=255)
                                                         gAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTCATGACGCCACCGCAAAAAACAAAGGCAAAACCACCGTTGACGGCAAGGAt                                                                                        >  2:2868328/1‑150 (MQ=18)
                                                            tCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                                                     >  1:1700385/1‑150 (MQ=255)
                                                            tCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCAAAAAACAGAGGCAAAACAACCGCACACGGCAAGCATAcc                                                                                     >  1:2521866/1‑150 (MQ=18)
                                                               aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                  <  2:2903502/150‑1 (MQ=255)
                                                               aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                  >  1:1528133/1‑150 (MQ=255)
                                                               aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                  >  1:2669041/1‑150 (MQ=255)
                                                                tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                                                 <  2:2765882/150‑1 (MQ=255)
                                                                tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                                                 <  2:2765889/150‑1 (MQ=255)
                                                                           gTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAAtt                                                                      >  1:1514794/1‑150 (MQ=255)
                                                                           gTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGACAtt                                                                      >  2:2589965/1‑150 (MQ=255)
                                                                           gTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAtt                                                                      >  1:2524243/1‑150 (MQ=255)
                                                                              gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCt                                                                   >  1:1833455/1‑150 (MQ=255)
                                                                               tcGGCGACCAGGCGGTTGTTAACCACGCAGGCGCGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTACCGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                                                  <  2:2613710/149‑1 (MQ=17)
                                                                               aCGGAGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                                                  <  2:2088247/150‑1 (MQ=255)
                                                                                 ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGc                                                                >  1:2459433/1‑150 (MQ=255)
                                                                                 ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAAAAACAAAGGAAAAACCACCGTTGAAGGCAAGGATTCCACCGGTACCGAAATAGCTGGc                                                                >  2:1213190/1‑150 (MQ=17)
                                                                                   cGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCaa                                                              >  1:1775783/1‑150 (MQ=255)
                                                                                    gACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGGGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCGAAAAACATCGGAAAAACCACCGTTGACGGCAAGGAATCCACCGGTACGGAAATAGCTggacat                                                             >  2:675019/1‑146 (MQ=17)
                                                                                       cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAAc                                                          >  2:1293272/1‑150 (MQ=255)
                                                                                       cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCTAAAAACAACGGCAAAACCACCGTTGAAGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAAc                                                          >  2:652748/1‑150 (MQ=17)
                                                                                                    acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                             >  1:172616/1‑150 (MQ=255)
                                                                                                    acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                             <  2:2572853/150‑1 (MQ=255)
                                                                                                     taacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCa                                            <  2:1061607/149‑1 (MQ=255)
                                                                                                     caacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCa                                            <  2:1061442/150‑1 (MQ=255)
                                                                                                          aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACg                                       >  1:1944512/1‑150 (MQ=255)
                                                                                                          aaGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                   >  1:1800022/1‑122 (MQ=25)
                                                                                                          aaGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                   <  2:1800022/122‑1 (MQ=25)
                                                                                                             gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                                          >  2:250271/1‑112 (MQ=25)
                                                                                                             gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                                          <  1:250271/112‑1 (MQ=39)
                                                                                                               gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt                                  >  1:1071426/1‑150 (MQ=255)
                                                                                                               gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt                                  >  1:1442281/1‑150 (MQ=255)
                                                                                                                 gagCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCGGGCAATAACGGGAAGGTGATTCAGGACGGCGATCt                                >  2:2660364/1‑150 (MQ=255)
                                                                                                                    cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGa                             >  1:677/1‑150 (MQ=255)
                                                                                                                     gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACACAGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAACGGGAAGGTGATTCAGGACGGAGAGCTGGAt                            >  1:2286912/1‑150 (MQ=17)
                                                                                                                     gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                            <  1:2657669/150‑1 (MQ=255)
                                                                                                                     gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                            <  1:2657670/150‑1 (MQ=255)
                                                                                                                     gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                            >  1:2152304/1‑150 (MQ=255)
                                                                                                                     gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                            >  1:2152370/1‑150 (MQ=255)
                                                                                                                     gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAGCTGGAt                            >  1:593270/1‑150 (MQ=255)
                                                                                                                      ccATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATg                           >  1:2158680/1‑150 (MQ=255)
                                                                                                                      ccATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATg                           >  1:1670877/1‑150 (MQ=255)
                                                                                                                       cATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                          >  1:348465/1‑150 (MQ=255)
                                                                                                                       cATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                          >  2:1221832/1‑150 (MQ=255)
                                                                                                                       cATCACCAACGGTGGCACCGGCACCCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                          >  2:1221647/1‑150 (MQ=255)
                                                                                                                        aTCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTc                         >  1:2549003/1‑150 (MQ=255)
                                                                                                                         tCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCa                        <  1:305008/150‑1 (MQ=255)
                                                                                                                           aCCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgc                      >  1:1382512/1‑150 (MQ=255)
                                                                                                                             cAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcgg                    >  1:136927/1‑150 (MQ=255)
                                                                                                                             cAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcgg                    >  1:131717/1‑150 (MQ=255)
                                                                                                                             cAACGGTGGCACCGGCACGCAGATTAACGGTGAGCACGCCACGGCACACAACAACGGCAAAACCACCGGTGACGGCCCGGAGTCCACCGGTACGGACAGTGCGGGCGCTAACGGGCCGGTGAGTCGGGACGGCGCTCGGGCGGTcggcgg                    >  1:2375744/1‑150 (MQ=9)
                                                                                                                              aaCGGTGGCTCCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGTCGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggc                   <  2:2782581/150‑1 (MQ=255)
                                                                                                                              aaCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggc                   >  2:266816/1‑150 (MQ=255)
                                                                                                                              aaCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggc                   >  2:1832187/1‑150 (MQ=255)
                                                                                                                              aaCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggc                   >  2:872097/1‑150 (MQ=255)
                                                                                                                              aaCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACCGCGAACAACAGAGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACCGAAATGGCGGGCAATAACGGCAGGGTCATTCAGGTCGGCGGTCAGGCGGGCAgcggc                   >  2:1213822/1‑150 (MQ=255)
                                                                                                                              aaCGGGGGCACCGGCACGCAGATTAACGGTGACGACGCCACGGCGAAGAACAGAGGCAAAACCACAGTTGACGGCAAGGATTCCACCGGTACCGAAATGGCGGGCAATATCGGCAAGGTGAATCAGGACGGCGGTCAGGCGGGCAgcggc                   >  2:1832161/1‑150 (MQ=14)
                                                                                                                                   tGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg              <  1:1221647/150‑1 (MQ=255)
                                                                                                                                   tGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg              <  1:1221832/150‑1 (MQ=255)
                                                                                                                                     gCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc            >  2:1616479/1‑150 (MQ=255)
                                                                                                                                      cACCGGCACGCAGATTACCGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAGTCCACCGGTACGGAAATGGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa           <  1:1938136/150‑1 (MQ=255)
                                                                                                                                      cACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCTCCGGTACTGAAATTGCTGGAAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa           <  1:2737085/150‑1 (MQ=255)
                                                                                                                                      cACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa           >  1:845927/1‑150 (MQ=255)
                                                                                                                                        ccGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg         <  1:2318074/150‑1 (MQ=255)
                                                                                                                                               gCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa  >  1:2711868/1‑150 (MQ=255)
                                                                                                                                               gCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa  <  2:1442281/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
ATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466487‑1466779

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: