Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,604 G→A 5.2% 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.2% 359.4 / 4.5 118pseudogene (1213/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (75/34);  new base A (2/4);  total (77/41)
Fisher's exact test for biased strand distribution p-value = 9.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

GGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAT  >  NC_000913/1466456‑1466753
                                                                                                                                                    |                                                                                                                                                     
ggATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGCTAGTAAAGAACGAAGGCGAGAgcgc                                                                                                                                                      >  2:2924869/1‑150 (MQ=255)
  aTGTCAGCCGCGGCGGACACGGTATTGATAACACCGGCGACAGCGCAACGGTGGATAACAAGGTCACCATTACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCa                                                                                                                                                    <  2:2702637/150‑1 (MQ=14)
   tGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCAt                                                                                                                                                   <  1:3583044/150‑1 (MQ=255)
    gTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATc                                                                                                                                                  <  2:1675647/150‑1 (MQ=255)
     tCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCa                                                                                                                                                 <  2:535342/150‑1 (MQ=255)
      cAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAc                                                                                                                                                <  2:2441067/150‑1 (MQ=255)
      cAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAc                                                                                                                                                <  2:2665223/150‑1 (MQ=255)
                 gTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCa                                                                                                                                     >  2:2438320/1‑150 (MQ=255)
                    aCGGTATTGATATCACCGTCTACAGCTCAACGGTGTATAACAAGGGCTCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCg                                                                                                                                  <  2:2979584/150‑1 (MQ=17)
                    aCGGTATTGATATCACCGGCGACAGCGCAACGGTGGATTACAAGGTCACCATGACCGTCTCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCg                                                                                                                                  <  2:2966235/150‑1 (MQ=255)
                    aCGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCg                                                                                                                                  >  1:405575/1‑150 (MQ=255)
                      ggTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGc                                                                                                                                >  2:471510/1‑150 (MQ=255)
                          ttGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGc                                                                                                                            >  2:2826843/1‑150 (MQ=255)
                              tatCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAt                                                                                                                        >  1:2081857/1‑150 (MQ=255)
                               atCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAtt                                                                                                                       >  2:19129/1‑150 (MQ=255)
                                tCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTa                                                                                                                      >  2:3227712/1‑150 (MQ=255)
                                  aCCGGCGACAGCGCAACGGTGGATAACAAGGGCTCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                                                                    <  2:1171265/150‑1 (MQ=255)
                                       tgACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                                                               <  2:2081857/149‑1 (MQ=255)
                                       agACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAACCGGAGTCCATCGGTATCCAGATAGACGGAGACAAGGCCATAGTTAAAAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                                                               <  2:2412639/149‑1 (MQ=11)
                                                aaCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCAc                                                                                                      >  1:1164004/1‑150 (MQ=255)
                                                 aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACg                                                                                                     >  2:718617/1‑150 (MQ=255)
                                                  cGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACgg                                                                                                    <  2:2660205/150‑1 (MQ=255)
                                                      ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaaa                                                                                                <  2:302810/150‑1 (MQ=255)
                                                      ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaaa                                                                                                >  2:191776/1‑150 (MQ=255)
                                                       gtaaaCAAGGGCACCACCACCGTCACCGACCCGGTGGGGAGCGGCTTCCACGTAGAGCGGGTCGACGGGTTTGTCAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaac                                                                                               <  2:3384434/147‑1 (MQ=9)
                                                        aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGGGGCACCGGCACGCAGATGAACGGTGAGGACGCCACCGCCaaca                                                                                              >  2:3249106/1‑150 (MQ=17)
                                                            cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTGAAGACGCCACGGCCaaaaacaa                                                                                          >  1:3699962/1‑150 (MQ=17)
                                                            cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTGAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATGAACGGTGAGGACGCCACGGCCaaaaacaa                                                                                          >  1:3699641/1‑150 (MQ=17)
                                                                 gCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAAGACGCCACCGCAAAAAACAACGGCa                                                                                     >  1:2238662/1‑150 (MQ=255)
                                                                     cATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATCAAc                                                                                                                    >  1:363435/1‑115 (MQ=255)
                                                                     cATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                 >  2:881609/1‑150 (MQ=255)
                                                                     cATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGGGGCACCGGCACTCAGATCAAc                                                                                                                    <  2:363435/115‑1 (MQ=255)
                                                                        gACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAaccac                                                                              <  2:2547269/150‑1 (MQ=255)
                                                                        gACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAaccac                                                                              >  1:1211883/1‑150 (MQ=255)
                                                                            gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAAGAACAACGGCAAAACCACCGtt                                                                          <  2:1164004/150‑1 (MQ=255)
                                                                            gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                          >  2:3550987/1‑150 (MQ=255)
                                                                             tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                         >  2:2506913/1‑150 (MQ=255)
                                                                             tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                         >  2:2506844/1‑150 (MQ=255)
                                                                               aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc                                                                       <  1:4122040/150‑1 (MQ=255)
                                                                                 cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                     >  2:1167398/1‑150 (MQ=255)
                                                                                 cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTGAACAACGAAGGCGAGAGAGAAATCAACAACTGTGGAAACGGCACGCATATTAACTGTGATGACGCCAAGGCAAAAAAAAACGGCAAAAACACCGTTGACgg                                                                     >  1:2096872/1‑150 (MQ=25)
                                                                                   aTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCa                                                                   >  1:2187455/1‑150 (MQ=255)
                                                                                   aTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGAGGCACCGGCACGCAGATAAACGGTTAAGACGCCACCGCGAAAAACAGAGGAAAAACCACCGTTGAAGGCa                                                                   >  2:1626746/1‑150 (MQ=14)
                                                                                   aTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGAGGCAACGGCACGCAGAAAAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCa                                                                   >  1:90860/1‑150 (MQ=18)
                                                                                    tCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                  <  2:523304/150‑1 (MQ=255)
                                                                                      cGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCACCATCACCAACGGAGGCACCGGCACTCAGATAAACGGTAAAGACGCCACCGCGAAAAACAGAGGCAAAACCACCGTTGAAGGaaaag                                                                >  2:356192/1‑148 (MQ=11)
                                                                                      cGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTAAACAACGAAGGCGAGAGCGCCATCACCAAAGGCGGCACCGGCACGCAGATAAACGGTGAGGACGCCACGGCGAAAAACAGTGGAAAAACCACAGTTGATGGaaaag                                                                >  2:1646846/1‑148 (MQ=11)
                                                                                       ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATTAACGGTGAAGACGCCACGGCAAACAACAACGGCAAAACAACAGATGACGGAAAGGa                                                               >  2:3731810/1‑150 (MQ=18)
                                                                                         aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCCCCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                             >  1:2049993/1‑150 (MQ=255)
                                                                                           tCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                           >  1:3259329/1‑150 (MQ=255)
                                                                                              aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                        >  1:908435/1‑150 (MQ=255)
                                                                                              aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAAAGGCAAAACCACAGTTGAAGGCAAGGATTCCAcc                                                        >  1:254402/1‑150 (MQ=255)
                                                                                              aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTCACGACGCCACCGCGAAAAACACAGGCAAAACCACCGTTGAAGGCAAGGATTCCAcc                                                        >  1:2167019/1‑150 (MQ=18)
                                                                                                cGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                      >  1:3165385/1‑150 (MQ=255)
                                                                                                cGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                      >  1:1289732/1‑150 (MQ=255)
                                                                                                cGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGGGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCGAAAAACAAAGGCAAAACCACCGTTGACGGCAAGGATTCCACggg                                                      >  1:1289631/1‑150 (MQ=18)
                                                                                                  gTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                    <  1:483322/150‑1 (MQ=255)
                                                                                                   tATCCAGTTTGACGGTGACAAGGCGTTTGTCAAGAACGAAGGCGAGAGCTCCATCACCAACGGTGGCAGCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAc                                                   <  2:804195/150‑1 (MQ=11)
                                                                                                       cAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                               >  1:2384997/1‑150 (MQ=255)
                                                                                                        aGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                              >  2:2013352/1‑150 (MQ=255)
                                                                                                        aGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                              >  2:2634913/1‑150 (MQ=255)
                                                                                                        aGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAAAAACAACGGCAAAACCACCGTTGAAGGAAATGATTCCACCGGTACGTaaa                                              >  1:1716223/1‑150 (MQ=18)
                                                                                                            tGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTgg                                          >  2:1763214/1‑149 (MQ=255)
                                                                                                               cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTgg                                       <  1:19129/150‑1 (MQ=255)
                                                                                                                ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGc                                      >  2:3066118/1‑150 (MQ=255)
                                                                                                                 gCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCa                                     >  2:666928/1‑150 (MQ=255)
                                                                                                                 gCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATAAACGGTGAGGACGCCACCGCAAAAAACAAGGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGGACCGAAATAGCTGCCa                                     >  2:3834408/1‑150 (MQ=255)
                                                                                                                  cGACCTGGCGGTTGTTAACATCGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCaa                                    <  1:1398279/150‑1 (MQ=255)
                                                                                                                    aCCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATa                                  >  2:3714627/1‑150 (MQ=255)
                                                                                                                      cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAAc                                >  2:3165655/1‑150 (MQ=255)
                                                                                                                      cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACCGCAAACAACAACGGAAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGCAAATAGCTGGCAATAAc                                >  2:4009007/1‑150 (MQ=17)
                                                                                                                           ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTGAAGACGCCACGGCAAAAAACAAAGGCAAAACCACAGTTGACGGAAAGGATTCCACCGGTACCAAAATAGCTGGCAATAACGGCaa                           >  2:2583110/1‑150 (MQ=17)
                                                                                                                           ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAACGGGaa                           >  2:2583164/1‑150 (MQ=255)
                                                                                                                           ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGAGGCACCGGCACTCAGATAAACGGTAAAGACGCCACCGCAAAAAACAGAGGAAAAACCACCGTTGAAGGAAAAGATTCCACCGGTACCAAAATAGCTGGCAATAACggcat                           >  1:2090004/1‑147 (MQ=255)
                                                                                                                            gttgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAg                          <  1:162028/150‑1 (MQ=255)
                                                                                                                                ttAAGAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGa                      <  2:2857755/150‑1 (MQ=255)
                                                                                                                                tcaataacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGa                      <  2:2660574/145‑1 (MQ=255)
                                                                                                                                 tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTCAAGACGCCACGGCGAAAAACAACGGAAAAACCACCGTTGACGGAAAGGATTCCACCGGACAGGAAATTGCAGGCAATAACGGGAAGGTGAt                     >  2:1926554/1‑150 (MQ=17)
                                                                                                                                  aacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAtt                    >  1:178548/1‑150 (MQ=255)
                                                                                                                                  aacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAtt                    >  1:1537654/1‑150 (MQ=255)
                                                                                                                                  aacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGAACGCAGATTAACGGTGATGACGCCTCGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAtt                    <  1:2283983/150‑1 (MQ=255)
                                                                                                                                    taacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                    >  2:1403192/2‑116 (MQ=14)
                                                                                                                                     aacGAAGTCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAg                 >  1:831060/1‑150 (MQ=255)
                                                                                                                                     aacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAg                 >  1:3529918/1‑150 (MQ=255)
                                                                                                                                     aacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAg                 >  1:3996116/1‑150 (MQ=255)
                                                                                                                                     aacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAGGACGCCACGGCAAACAACAACGGAAAAACCACCGTTGAAGGCAAGGATTCCACCGGTACCGAAATAGCTGGCAATATCGGCAAGGTGATTCAg                 >  2:2595001/1‑150 (MQ=17)
                                                                                                                                       cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCa                                     >  1:3463959/1‑128 (MQ=255)
                                                                                                                                       cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCa                                     <  2:3463959/128‑1 (MQ=255)
                                                                                                                                       cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa               <  2:3724805/150‑1 (MQ=255)
                                                                                                                                         aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACg             <  1:471510/150‑1 (MQ=255)
                                                                                                                                         aaGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACg             <  1:2054372/150‑1 (MQ=255)
                                                                                                                                            gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg          >  2:1086059/1‑150 (MQ=255)
                                                                                                                                             cGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa         >  2:290260/1‑150 (MQ=255)
                                                                                                                                                gagCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCt      <  1:3741924/150‑1 (MQ=255)
                                                                                                                                                gagCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCt      <  1:835331/150‑1 (MQ=255)
                                                                                                                                                 agCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGATGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTg     >  2:884499/1‑150 (MQ=255)
                                                                                                                                                 agCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGATGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTg     >  2:884504/1‑150 (MQ=255)
                                                                                                                                                 agCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTg     >  2:3218067/1‑150 (MQ=255)
                                                                                                                                                 agCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAAAAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATGGCTGGCAATAACGGCAAGGTGATTCAGGACGGCGATCTg     >  1:2352248/1‑150 (MQ=255)
                                                                                                                                                    gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt  <  1:2438320/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAT  >  NC_000913/1466456‑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: