Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,630 G→C 7.3% 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 7.3% 837.0 / 44.9 257pseudogene (1239/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (148/86);  new base C (9/10);  total (160/97)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

TTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466482‑1466779
                                                                                                                                                    |                                                                                                                                                     
ttGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGc                                                                                                                                                      >  2:4512641/1‑150 (MQ=255)
    tatCACCGGCGACAGCGCATCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAt                                                                                                                                                  <  2:3536594/150‑1 (MQ=255)
    tatCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAt                                                                                                                                                  >  1:4629440/1‑150 (MQ=255)
    tatCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAt                                                                                                                                                  <  2:679137/150‑1 (MQ=255)
       cACCGGCGTGTGCGCACCGGTGGCTATCGCGGGCTCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGTTCGACGGCGACCAGGCGGTCGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTaa                                                                                                                                               <  2:4997384/150‑1 (MQ=14)
       cACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTaa                                                                                                                                               <  1:4176942/150‑1 (MQ=255)
        aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                                                                                              <  2:5212522/150‑1 (MQ=255)
         ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                                                                                             <  2:5172040/150‑1 (MQ=255)
           ggCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCTATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGt                                                                                                                                           <  2:1581/150‑1 (MQ=255)
              gACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgat                                                                                                                                        >  1:1946509/1‑150 (MQ=255)
                          gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCa                                                                                                                            >  1:4680588/1‑150 (MQ=255)
                            ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaaa                                                                                                                          <  2:6666767/150‑1 (MQ=255)
                             gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGGGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCCaac                                                                                                                         >  1:4958084/1‑150 (MQ=17)
                              aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaaca                                                                                                                        >  1:5732277/1‑150 (MQ=255)
                                aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaac                                                                                                                      >  2:1310069/1‑150 (MQ=255)
                                  cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaacaa                                                                                                                    >  2:2959379/1‑150 (MQ=255)
                                    aGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTAAAGACGCCACCGCGaaaaacaaag                                                                                                                  >  1:945223/1‑148 (MQ=14)
                                    aGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACg                                                                                                                  <  2:5131658/150‑1 (MQ=255)
                                    aGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACg                                                                                                                  >  1:1702785/1‑150 (MQ=255)
                                        cACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGTGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaa                                                                                                              >  2:3223825/1‑150 (MQ=255)
                                         aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATAAACGGTGAAGACGCCACGGCAAAAAACAACGGCaaa                                                                                                             >  2:2308703/1‑150 (MQ=255)
                                         aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGGGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCAAAAAACAAAGGCaaa                                                                                                             >  2:3351940/1‑150 (MQ=18)
                                             tGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                                                                                                         >  1:1187228/1‑150 (MQ=255)
                                                 cGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGt                                                                                                     >  1:6568272/1‑150 (MQ=255)
                                                   tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                                                   >  1:6317139/1‑150 (MQ=255)
                                                     agagATCCGGAGTCCATCGGTATCCAGTTTGACGGCGACCAGGCGTTTGTTAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc                                                                                                 <  2:2887445/147‑1 (MQ=17)
                                                     aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCGAACAACAACGGCAAAACCACCGTTGAc                                                                                                 >  1:5191968/1‑150 (MQ=255)
                                                       cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                               <  2:3392465/150‑1 (MQ=255)
                                                          tCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTa                                                                                                                                                >  2:5279029/1‑98 (MQ=255)
                                                          tCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTa                                                                                                                                                <  1:5279029/98‑1 (MQ=255)
                                                          tCCGGAGTCCATCGGTATCCAGGTCGACGGCGACCAGGCGGTTGTTAACAACGACGGCGAGAGCGCCAACACCCACGGTGGCAACTGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                            <  2:447619/150‑1 (MQ=17)
                                                             ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGGTGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGa                                                                                         <  1:3254333/150‑1 (MQ=255)
                                                                    aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                  >  1:5919907/1‑150 (MQ=255)
                                                                    aTCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTTAAAAACGAAGGCGATAGCGACATCACCAACGCTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCTAACAACAACGGAAAAACCACCGTTGAAGGCAAGGATTCCAcc                                                                                  >  1:3920184/1‑150 (MQ=255)
                                                                        gTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCAAAAAACAAAGGCAAAACCACCGTTGAAGGCAAGGATTCCACCGGTa                                                                              >  1:2876951/1‑150 (MQ=18)
                                                                           tCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATGAACGGTGATGACGCCACGGCAAACAACAACGGCACAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                           >  2:3972465/1‑150 (MQ=255)
                                                                            ccAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCGAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                                          >  1:3111053/1‑150 (MQ=255)
                                                                            ccAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                                          >  1:274352/1‑150 (MQ=255)
                                                                              aGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                                                        >  2:1256791/1‑150 (MQ=255)
                                                                                    aCGGTGACAAGGCGTTAGTGAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                                                  <  2:1385851/150‑1 (MQ=11)
                                                                                    aCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                                                  <  2:3904008/150‑1 (MQ=255)
                                                                                    aCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                                                  <  1:2499380/150‑1 (MQ=255)
                                                                                        cGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCaa                                                              <  2:6317139/150‑1 (MQ=255)
                                                                                            cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAAc                                                          >  1:6111425/1‑150 (MQ=255)
                                                                                            cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAAc                                                          >  1:2449526/1‑150 (MQ=255)
                                                                                               gCGGTCGTTAAGAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACggg                                                       <  2:6021739/150‑1 (MQ=255)
                                                                                                 ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCACCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaa                                                     <  1:2840835/150‑1 (MQ=255)
                                                                                                 ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaa                                                     <  2:2500604/150‑1 (MQ=255)
                                                                                                 ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaa                                                     >  2:419539/1‑150 (MQ=255)
                                                                                                  gttgttAACACCGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTACCGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATACCGGGAAg                                                    <  2:5757738/150‑1 (MQ=255)
                                                                                                  gttgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAGCGGAAAAACCAGAGGTGACGGCCAGGATTCCACCGGTCCGGGAGTTGCTGGCAGTAACGCGacg                                                    >  2:5572169/1‑148 (MQ=17)
                                                                                                   ttgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAgg                                                   >  1:3017019/1‑150 (MQ=255)
                                                                                                    tgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGt                                                  <  2:4911959/150‑1 (MQ=255)
                                                                                                    tgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGt                                                  >  1:3894397/1‑150 (MQ=255)
                                                                                                    tgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGt                                                  <  2:6496226/150‑1 (MQ=255)
                                                                                                    tctcaagaacGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACGCAGATTAATGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGt                                                  <  2:359744/143‑1 (MQ=14)
                                                                                                     gttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCAACAGGTACGGAAATTGCTGGCAATAACGGGAAGGTg                                                 >  1:3519443/1‑150 (MQ=255)
                                                                                                      ttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGa                                                >  2:465293/1‑150 (MQ=255)
                                                                                                       tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGCTGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAt                                               >  1:534019/1‑150 (MQ=255)
                                                                                                       tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCACTAACGGGAAGGTGAt                                               >  2:3070422/1‑150 (MQ=255)
                                                                                                       tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAt                                               >  1:3085615/1‑150 (MQ=255)
                                                                                                       tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAt                                               >  2:840179/1‑150 (MQ=255)
                                                                                                        aataacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGa                                                      >  1:1831941/4‑142 (MQ=255)
                                                                                                         acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                             >  1:367078/1‑150 (MQ=255)
                                                                                                         acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACAGTTGACGGAAAGGATTCCACCGGTACCGAAATAGCTGGCAATAACGGCAAGGTGATTc                                             >  2:1531752/1‑150 (MQ=17)
                                                                                                         acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATGAACGGTGAAGACGCCACCGCAAAAAACAGCGGCAAAACCACAGTTGACGGCAAGGATTCCACCGGTACGAAAATAGCTGGCAATAACGGCAAGGAGAtgc                                             >  1:3589374/1‑148 (MQ=17)
                                                                                                          taacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                                                                                         <  2:1124162/56‑1 (MQ=255)
                                                                                                          taacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                                                                                         >  1:1124162/2‑57 (MQ=255)
                                                                                                           aacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAAGACGCCACGGCTAACAACAACGGAAAAACCACAGTTGACGGAAATGATTCCACCGGTACGGAAATAGCTGGCAATAACGGCAAGGAGATACac                                           >  1:5052874/1‑149 (MQ=17)
                                                                                                             cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa                                         <  2:2681188/150‑1 (MQ=255)
                                                                                                              gAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAACGAATCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAc                                        >  1:4574808/1‑150 (MQ=255)
                                                                                                               aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACg                                       >  1:5540173/1‑150 (MQ=255)
                                                                                                                  gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg                                    >  1:5608282/1‑150 (MQ=255)
                                                                                                                        gcgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTgg                              <  1:1310069/150‑1 (MQ=255)
                                                                                                                         cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGa                             <  2:3195246/150‑1 (MQ=255)
                                                                                                                         cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGa                             <  1:2210150/150‑1 (MQ=255)
                                                                                                                         cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAACGGTACGGAAATTGCTTGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGa                             <  1:695473/150‑1 (MQ=255)
                                                                                                                          gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                            >  1:4135695/1‑150 (MQ=255)
                                                                                                                          accATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                            <  1:14700/149‑1 (MQ=255)
                                                                                                                            cATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                   <  1:3745697/93‑1 (MQ=37)
                                                                                                                            cATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                          >  2:373931/1‑150 (MQ=255)
                                                                                                                             aTCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAGCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAAGACCGGGAAGGTGATTCAGGACGGCGATGTGGTTGTc                         >  2:4055281/1‑150 (MQ=17)
                                                                                                                              tCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCa                        <  2:2232989/150‑1 (MQ=255)
                                                                                                                                aCCAACGGTGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCGAAAAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAACGGCAAGGTGATTCAGGACGGCGATCTCGATGtagga                      >  1:4523668/1‑146 (MQ=17)
                                                                                                                                 ccAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                            >  1:4014965/1‑143 (MQ=255)
                                                                                                                                   aaCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggc                   >  1:5748338/1‑150 (MQ=255)
                                                                                                                                       gTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcggcg               <  2:2654733/150‑1 (MQ=255)
                                                                                                                                        tGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg              <  2:6568272/150‑1 (MQ=255)
                                                                                                                                          gCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc            >  1:4441220/1‑150 (MQ=255)
                                                                                                                                          gCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc            >  2:257132/1‑150 (MQ=255)
                                                                                                                                            aCCGGCACGCAGATTAACGGTGATGACGCCACGCCAAACAACAACGGCAAAACCACCGTTGACCGCAAGCATACCACCGGTCCGGACATTGCTGCCAATAACGGGCAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc          >  1:3416984/1‑150 (MQ=17)
                                                                                                                                             ccGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg         >  1:370995/1‑150 (MQ=255)
                                                                                                                                              cGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg        >  1:5815444/1‑150 (MQ=255)
                                                                                                                                                gCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTa      >  1:6174361/1‑150 (MQ=255)
                                                                                                                                                gCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTa      >  1:6536734/1‑150 (MQ=255)
                                                                                                                                                gCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTa      >  1:6604389/1‑150 (MQ=255)
                                                                                                                                                  aCGCAGATTAACGGTGCTGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt    <  2:534019/150‑1 (MQ=255)
                                                                                                                                                  aCGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTg                                                                     >  2:1856867/1‑85 (MQ=40)
                                                                                                                                                    gCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa  >  1:3300175/1‑150 (MQ=255)
                                                                                                                                                    gCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa  >  1:4856183/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466482‑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: