Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,639 C→T 6.0% pseudogene (1248/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,6390CT 6.0% 1129.0 / 34.0 366pseudogene (1248/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (200/139);  new base T (9/12);  total (212/154)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 376 total aligned reads displayed.

ACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGG  >  NC_000913/1466490‑1466788
                                                                                                                                                     |                                                                                                                                                     
aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                                                                                                       <  2:5413618/150‑1 (MQ=255)
aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                                                                                                       <  1:1333562/150‑1 (MQ=255)
aCCGGCGAAAGCGCAACGGTGCATAACAAGGGCACCAAGACAGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACAAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                                                                                                       <  2:580347/150‑1 (MQ=18)
 ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                                                                                                      >  1:1090004/1‑150 (MQ=255)
  cGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACgg                                                                                                                                                     >  1:8751327/1‑150 (MQ=255)
     cGACAGCGCAACGGTTGATATCAAGGGCACCTTTTCAGTCACGGTTCCGGTGTCCATCGGTATCCAGGTTGACGGAGACAAGGCGTTTGTTAACAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                                                                                                  <  1:1725357/150‑1 (MQ=17)
     cGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                                                                                                  <  2:3048051/150‑1 (MQ=255)
      gACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgat                                                                                                                                                 >  2:8275875/1‑150 (MQ=255)
      gACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAGCGGTCTCCAGGTTGACGGCGCCCAGGCGGGTGGTCACAACGAAGGCGAGAGCGCCATCACCACCGGTGGCCCCGGCACGCCGAGTAACGGggat                                                                                                                                                 >  1:4866425/1‑150 (MQ=18)
         aGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATCAACGGTAACGAc                                                                                                                                              >  1:2001742/1‑150 (MQ=255)
         aGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGAc                                                                                                                                              >  2:2057565/1‑150 (MQ=255)
           cgcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGc                                                                                                                                            <  1:2663431/150‑1 (MQ=255)
                cGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACgg                                                                                                                                       >  1:4569795/1‑150 (MQ=255)
                  gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCa                                                                                                                                     >  2:6144824/1‑150 (MQ=255)
                    ggATAACAAGGGCACCAGGACCGTCACCGATCCGGAGTCCATGGGTATCCAGCTAGACGGCGACAAGGCGCTGGTCAACAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaaa                                                                                                                                   <  1:8223234/150‑1 (MQ=11)
                     gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaac                                                                                                                                  >  2:9646970/1‑150 (MQ=255)
                           aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaacaac                                                                                                                            <  1:6500753/150‑1 (MQ=255)
                              ggCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGc                                                                                                                         >  1:10024533/1‑150 (MQ=255)
                                           tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAAAAACAGAGGAAAAACCACTGTTg                                                                                                            >  2:3491008/1‑150 (MQ=11)
                                           tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTTAAGACGCCACCGCAAAAAACAAAGGCAAAACCACAGTTg                                                                                                            >  2:5146823/1‑150 (MQ=18)
                                             aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc                                                                                                          >  1:1129056/1‑150 (MQ=255)
                                             aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCCAAAAACAGCGGCAAGAACGCCGTAGAc                                                                                                          >  1:4968128/1‑150 (MQ=18)
                                              ccGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                                                                                                  <  1:3695360/109‑1 (MQ=255)
                                               cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                                        >  1:5069589/1‑150 (MQ=255)
                                                 aTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCa                                                                                                      >  2:9507114/1‑150 (MQ=255)
                                                 aTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCGTAGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAAGACGCCACAGCCAAAAACAAAGGCAAAACCACCGTTGAAggaa                                                                                                      >  2:9355798/1‑148 (MQ=17)
                                                  tccgttgTCCATGGGTATCCAGCTCGACGGCGACCAGGCCTTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                                     <  2:7323497/144‑1 (MQ=14)
                                                  tCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                                     <  2:10109166/150‑1 (MQ=255)
                                                  tCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                                     <  1:5573683/150‑1 (MQ=255)
                                                  tCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                                                     <  2:6243911/150‑1 (MQ=255)
                                                        gTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTc                                                                                               <  1:8898373/150‑1 (MQ=255)
                                                        gTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTc                                                                                               <  2:1539395/150‑1 (MQ=255)
                                                          ccATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                                                             <  1:98898/150‑1 (MQ=255)
                                                           cagcGGGATCCAGGAGGCCGGAGACCAGGCGGTTGCTAACAACGACGGCGAGAGCCCCATCACCTACGTTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCTAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                            <  2:8109520/147‑1 (MQ=17)
                                                            aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAAAGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                           >  2:4793482/1‑150 (MQ=255)
                                                             tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                                                          <  1:5332394/150‑1 (MQ=255)
                                                             tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                                                          <  1:4155286/150‑1 (MQ=255)
                                                             tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAAGACGCCACGGCAAACAACAACGGCAATACCACCTATGACGGCACGGATTCCACCg                                                                                          >  1:5811788/1‑150 (MQ=18)
                                                               ggTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGTGGCACCGGCACGCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGTCAAGGATTCCACCGGt                                                                                        >  1:7623425/1‑150 (MQ=255)
                                                                gTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                                       >  1:210655/1‑150 (MQ=255)
                                                                gTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                                       >  2:7060190/1‑150 (MQ=255)
                                                                gTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAAGACGCCACCGCAAAAAACAAAGGCAAAACCACAGTTGAAGGAAATGATTCCACCGGTa                                                                                       >  2:2839920/1‑150 (MQ=18)
                                                                gTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGAGGCACCGGCACTCAGATAAACGGTCAAGACGCCACCGCGAAAAACAAAGGAAAAACCACAGTTGAAGGAAAGGATTCCACCGGTa                                                                                       >  2:3811705/1‑150 (MQ=255)
                                                                   tCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCGAAAAACAACGGCAAAACCACCGTTGACGGACAGGATTCCACCGGTACgg                                                                                    >  1:4596091/1‑150 (MQ=18)
                                                                           gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                            <  2:6758801/150‑1 (MQ=255)
                                                                            aCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgatga                                                                                                                                               <  2:1119217/82‑1 (MQ=255)
                                                                              ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCACCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGc                                                                         <  1:6584794/150‑1 (MQ=255)
                                                                               gCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCa                                                                        >  1:6198854/1‑150 (MQ=255)
                                                                                    cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAAc                                                                   >  1:6938192/1‑150 (MQ=255)
                                                                                    cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAAc                                                                   <  1:145052/150‑1 (MQ=255)
                                                                                         ggTTGTTATCACCTACGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaa                                                              <  2:8714158/150‑1 (MQ=17)
                                                                                         ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaa                                                              >  2:2670387/1‑150 (MQ=255)
                                                                                            tgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATCACGGGAAGGt                                                           >  1:5741047/1‑150 (MQ=255)
                                                                                              ttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGa                                                         >  2:8563863/1‑150 (MQ=255)
                                                                                              tcacccacgacgGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCCGTTTCACGGTGATGTCGTCACGGCTATCGACCACGGGATAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGa                                                         <  1:9577665/144‑1 (MQ=17)
                                                                                                aacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCACGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAtt                                                       >  1:3395290/1‑150 (MQ=255)
                                                                                                 ataacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGACAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                                      >  1:8827919/3‑150 (MQ=255)
                                                                                                 acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCACTAACGGGACGGTGAGTc                                                      >  1:10106677/1‑150 (MQ=255)
                                                                                                 acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                                      >  1:8099894/1‑150 (MQ=255)
                                                                                                   aacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAg                                                    <  2:5069589/150‑1 (MQ=255)
                                                                                                     cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa                                                  <  2:8416027/150‑1 (MQ=255)
                                                                                                      gAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAc                                                 >  1:2728705/1‑150 (MQ=255)
                                                                                                      gAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAc                                                 >  2:4320043/1‑150 (MQ=255)
                                                                                                       aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                                                                  >  1:5496433/1‑116 (MQ=39)
                                                                                                       aaGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAtt                                                       <  1:8308885/143‑1 (MQ=255)
                                                                                                        aGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACgg                                               <  1:94255/150‑1 (MQ=255)
                                                                                                          gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCCCTCACGGGACGGGGATTCAGGACGGCg                                             >  1:3050352/1‑150 (MQ=17)
                                                                                                          gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg                                             >  1:2842877/1‑150 (MQ=255)
                                                                                                          gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg                                             >  1:9272462/1‑150 (MQ=255)
                                                                                                           ctagagCGCCATCACCAACGTTGTCAACGGCACGAAGATTAACGTTGATGACGCCACGGCAAACAAGAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGa                                            <  2:834065/148‑1 (MQ=17)
                                                                                                             agagCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                          >  2:9204837/1‑54 (MQ=255)
                                                                                                                    cATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                                   >  1:4006970/1‑150 (MQ=255)
                                                                                                                    cATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                                   >  2:5625780/1‑150 (MQ=255)
                                                                                                                     aTCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTc                                  >  2:4547183/1‑150 (MQ=255)
                                                                                                                        aCCAACGGTGGCACCGGCACGCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAtt                                                       <  1:9381773/126‑1 (MQ=255)
                                                                                                                        aCCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgc                               >  1:9716556/1‑150 (MQ=255)
                                                                                                                        aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgc                               >  2:9311918/1‑150 (MQ=255)
                                                                                                                        aCCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTgg                                                                          >  2:9068181/1‑107 (MQ=25)
                                                                                                                         ccAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcg                              <  1:9203879/150‑1 (MQ=255)
                                                                                                                         ccAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcg                              >  2:1819572/1‑150 (MQ=255)
                                                                                                                         ccAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcg                              <  1:9863608/150‑1 (MQ=255)
                                                                                                                          cAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcgg                             >  1:2261694/1‑150 (MQ=255)
                                                                                                                          cAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGAACCGAAATTGCTGGCAATAACGGGACGGTGATTCATCACGGCTATCTGGATGTCAgcgg                             >  1:5999584/1‑150 (MQ=17)
                                                                                                                          cAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCAAAAAACAACGTAAAAACCACCAATCACGGCAAGTAAACAAACTTTACGGAAATTGCTGGCAATAACGGGAAGTTGATTCAGGACGGCGATCTGGATGTCAgcgg                             >  1:7254243/1‑150 (MQ=17)
                                                                                                                            aCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcg                           >  1:3794204/1‑150 (MQ=255)
                                                                                                                             cGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                          <  1:7806543/150‑1 (MQ=255)
                                                                                                                             cGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                          <  2:2948582/150‑1 (MQ=255)
                                                                                                                               gTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcggcg                        >  2:2098406/1‑150 (MQ=255)
                                                                                                                                   cACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa                    <  1:8563863/150‑1 (MQ=255)
                                                                                                                                    aCCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc                   <  2:7871208/150‑1 (MQ=255)
                                                                                                                                    aCCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCCTTGACGGCAAGGAGTCCAACGGTACGGAAATTGCCGGCAATAACGGGAATGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc                   <  2:1406983/150‑1 (MQ=14)
                                                                                                                                      cGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg                 >  2:5535982/1‑150 (MQ=255)
                                                                                                                                      cGGCACCCAGATTAATGGTGATGACGCCACGTCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg                 <  2:4242075/150‑1 (MQ=255)
                                                                                                                                          aCGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt             >  1:8032580/1‑150 (MQ=255)
                                                                                                                                          aCGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt             >  1:1514370/1‑150 (MQ=255)
                                                                                                                                            gCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa           >  1:4674685/1‑150 (MQ=255)
                                                                                                                                            ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                          >  2:1195860/2‑135 (MQ=255)
                                                                                                                                                 ttAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCa      >  1:9532284/1‑150 (MQ=255)
                                                                                                                                                  tAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAc     >  2:6611639/1‑150 (MQ=255)
                                                                                                                                                     cGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg  <  2:4031125/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGG  >  NC_000913/1466490‑1466788

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: