Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,504 A→G 5.2% pseudogene (1113/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,5040AG 5.2% 492.9 / 9.0 156pseudogene (1113/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (62/82);  new base G (5/3);  total (67/89)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 162 total aligned reads displayed.

CCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCAC  >  NC_000913/1466356‑1466653
                                                                                                                                                    |                                                                                                                                                      
ccACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCaa                                                                                                                                                       <  1:3571173/150‑1 (MQ=255)
ccACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCaa                                                                                                                                                       >  1:3946797/1‑150 (MQ=255)
     gtaaaCAACAATGGCAAAACCACCGTTGACGGTAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTg                                                                                                                                                  <  1:4468999/148‑1 (MQ=255)
     gCAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGTGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTg                                                                                                                                                  <  1:1556117/150‑1 (MQ=255)
       aaaCAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGa                                                                                                                                                >  1:4049089/1‑150 (MQ=255)
       aaaCAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGa                                                                                                                                                >  1:454368/1‑150 (MQ=255)
        aacaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGAt                                                                                                                                               >  1:4776762/1‑150 (MQ=255)
            acaccGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACa                                                                                                                                           <  2:3944115/150‑1 (MQ=255)
            acaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACa                                                                                                                                           >  1:4221247/1‑150 (MQ=255)
                tggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAggg                                                                                                                                       <  2:2296797/149‑1 (MQ=18)
                cGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAggg                                                                                                                                       >  1:2495613/1‑150 (MQ=255)
                 ggCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGc                                                                                                                                      <  1:2766085/150‑1 (MQ=255)
                  gCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCa                                                                                                                                     >  2:4052662/1‑150 (MQ=255)
                    aaaaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCAcc                                                                                                                                   <  1:3656653/150‑1 (MQ=255)
                     aaaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCa                                                                                                                                  >  1:899282/1‑150 (MQ=255)
                        ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCACCGTAGGCGGTCACGGTATTGATAACACCGGCGACAGCGCAACGGAGCATAACAAGGGCACCAT‑Ga                                                                                                                              >  2:1967184/1‑150 (MQ=14)
                         caccGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACaa                                                                                                                                          <  1:4797920/138‑1 (MQ=255)
                         caccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GAc                                                                                                                             >  1:3836584/1‑150 (MQ=255)
                          accGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GAcc                                                                                                                            <  1:4052662/150‑1 (MQ=255)
                          accGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GAcc                                                                                                                            <  1:3353624/150‑1 (MQ=255)
                              ttGACGGCAAGGATTCCACCGGTACGGAAATCTACGGTAATAACGGTAATGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCa                                                                                                                        <  2:3836584/150‑1 (MQ=255)
                              ttGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCa                                                                                                                        <  2:883539/150‑1 (MQ=255)
                              ttGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCa                                                                                                                        <  1:2744599/150‑1 (MQ=255)
                              ttGACGGAAACGATTCCACCGGTACGGAAATCAATGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCa                                                                                                                        <  1:2148079/150‑1 (MQ=255)
                               tGACGGCAAGGATTCCACCGGTACGGAAATCCCCGGTAATAACGGCAATGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCAc                                                                                                                       <  2:2495613/150‑1 (MQ=11)
                                   ggCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGAt                                                                                                                   >  2:4527998/1‑150 (MQ=255)
                                   ggCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGAt                                                                                                                   >  2:2682762/1‑150 (MQ=255)
                                     cAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATcc                                                                                                                 >  2:355635/1‑150 (MQ=255)
                                     cAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATcc                                                                                                                 >  1:922672/1‑150 (MQ=255)
                                     cAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATcc                                                                                                                 <  2:1750988/150‑1 (MQ=255)
                                      aaGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCg                                                                                                                <  2:2227117/150‑1 (MQ=255)
                                       aagATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCgg                                                                                                               <  1:4351554/148‑1 (MQ=25)
                                       aGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCgg                                                                                                               >  1:5659984/1‑150 (MQ=40)
                                        ggATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGa                                                                                                              <  1:5353273/150‑1 (MQ=40)
                                              cACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCAt                                                                                                        <  1:2443307/150‑1 (MQ=39)
                                               aCCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATc                                                                                                       <  2:646389/150‑1 (MQ=40)
                                               aCCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATc                                                                                                       >  2:4717075/1‑150 (MQ=40)
                                               aCCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATc                                                                                                       >  1:3724605/1‑150 (MQ=40)
                                                      cGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATcc                                                                                                >  2:1710002/1‑150 (MQ=255)
                                                       ggAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCa                                                                                               <  1:4728108/150‑1 (MQ=17)
                                                       ggAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCa                                                                                               <  2:4776808/150‑1 (MQ=255)
                                                         aaaTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAgg                                                                                             <  2:1887887/150‑1 (MQ=255)
                                                           aTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGtt                                                                                           <  2:4328167/150‑1 (MQ=255)
                                                               aCGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACg                                                                                       <  2:1756422/150‑1 (MQ=38)
                                                                 ggTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGc                                                                                     <  1:3370067/150‑1 (MQ=255)
                                                                 ggTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGc                                                                                     <  1:3463478/150‑1 (MQ=255)
                                                                 ggTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGc                                                                                     >  1:838399/1‑150 (MQ=255)
                                                                 ggTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGc                                                                                     >  1:5797826/1‑150 (MQ=37)
                                                                     ataaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGAcc                                                                                 >  1:5138875/1‑150 (MQ=39)
                                                                     ataaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGAcc                                                                                 >  1:3517934/1‑150 (MQ=21)
                                                                      taaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCa                                                                                >  1:3925807/1‑150 (MQ=39)
                                                                       aaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAg                                                                               >  2:4632396/1‑150 (MQ=39)
                                                                       aaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAg                                                                               >  1:2408679/1‑150 (MQ=39)
                                                                          ggAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCg                                                                            >  1:1269393/1‑150 (MQ=40)
                                                                          ggAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGATGACGGCGACCAGGCg                                                                            >  2:4572837/1‑150 (MQ=255)
                                                                            aaaaGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgt                                                                          <  2:3435134/150‑1 (MQ=40)
                                                                               aGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgt                                                                       <  2:5315670/150‑1 (MQ=40)
                                                                               aGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgt                                                                       <  2:1118906/150‑1 (MQ=40)
                                                                               aGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgt                                                                       >  2:3712374/1‑150 (MQ=255)
                                                                                 tgaTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCg                                                                            >  1:3736635/3‑143 (MQ=9)
                                                                                 tgaTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGGGACAGCGCTACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCg                                                                            <  2:3736635/141‑1 (MQ=9)
                                                                                  tATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCa                                                                                                                        <  2:2026372/98‑1 (MQ=31)
                                                                                  tATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCa                                                                                                                        >  1:2026372/1‑98 (MQ=31)
                                                                                  tATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCgg                                                                                                               <  1:2140896/107‑1 (MQ=11)
                                                                                   aTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaac                                                                   <  2:1107807/150‑1 (MQ=40)
                                                                                    ttCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaaca                                                                  <  2:514758/150‑1 (MQ=255)
                                                                                    ttCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaaca                                                                  >  1:5002282/1‑150 (MQ=40)
                                                                                     tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaa                                                                 <  2:3343957/150‑1 (MQ=255)
                                                                                     tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaa                                                                 <  2:3949042/150‑1 (MQ=255)
                                                                                        ggACGGTGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGa                                                              >  1:1556855/1‑150 (MQ=255)
                                                                                         gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGaa                                                             <  2:4804616/150‑1 (MQ=255)
                                                                                               gATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCTACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGa                                                                                   <  2:1176719/122‑1 (MQ=17)
                                                                                                aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCgaga                                                      <  2:4049089/150‑1 (MQ=255)
                                                                                                  cTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgc                                                    <  1:5528981/150‑1 (MQ=255)
                                                                                                     gATGTGAGCGGCGGCGGTCACGGATCTCATCACCCCGGCGACATCGGACCGGTGGCTACGCAGGGC‑CGATCGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGcc                                                 <  2:2805586/150‑1 (MQ=9)
                                                                                                      aTGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCa                                                >  1:1632670/1‑150 (MQ=255)
                                                                                                         tCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCa                                             >  1:556475/1‑150 (MQ=255)
                                                                                                         tCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCa                                             <  1:2265581/150‑1 (MQ=255)
                                                                                                            gcggcggcgGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCAGCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCa                                          <  2:4411599/150‑1 (MQ=255)
                                                                                                             cggcggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCaa                                         <  2:899282/150‑1 (MQ=255)
                                                                                                                 ggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGt                                     <  2:884003/150‑1 (MQ=255)
                                                                                                                   cggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTgg                                   >  1:2612443/1‑150 (MQ=255)
                                                                                                                             aTTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGAGACACCGGCAct                         >  2:3938101/1‑149 (MQ=255)
                                                                                                                                gATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAg                      >  1:1065269/1‑150 (MQ=255)
                                                                                                                                    tCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTa                  >  1:101613/1‑150 (MQ=255)
                                                                                                                                      aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                >  2:1367026/1‑150 (MQ=255)
                                                                                                                                      aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                >  1:5501548/1‑150 (MQ=255)
                                                                                                                                       ccGGTGACAGCGCGACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg               >  1:767796/1‑150 (MQ=255)
                                                                                                                                       ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg               <  2:4664034/150‑1 (MQ=255)
                                                                                                                                       ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg               <  1:319867/150‑1 (MQ=255)
                                                                                                                                       ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg               <  2:5659984/150‑1 (MQ=255)
                                                                                                                                          gCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtg            >  1:3434674/1‑150 (MQ=255)
                                                                                                                                           cGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga           >  1:1523205/1‑150 (MQ=255)
                                                                                                                                             aCAGCGCGACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGAc                                                                                        <  1:2498416/71‑1 (MQ=21)
                                                                                                                                              cAGCGCAACGTGTGAGAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgatga        <  2:5675872/150‑1 (MQ=255)
                                                                                                                                                 cgcAACTGTTGATAACAAGGGCACCAT‑GACCGTCACAGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGTTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGc     <  2:4846934/150‑1 (MQ=255)
                                                                                                                                                 cgcAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGTTGGTTAACAACGAAGGCgaga                                                      <  2:5262694/101‑1 (MQ=255)
                                                                                                                                                 cgcAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCgaga                                                      >  1:5262694/1‑101 (MQ=255)
                                                                                                                                                    aaCGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCAc  >  1:1776357/1‑150 (MQ=255)
                                                                                                                                                    aaCGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCAc  >  1:1776331/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT‑GACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCAC  >  NC_000913/1466356‑1466653

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