Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,505 A→C 8.1% R194R (CGT→CGG) ‡ ftsX ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,5050AC 8.1% 300.1 / 14.1 108R194R (CGT→CGG) ‡ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (49/50);  new base C (5/4);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

Only 100 of 109 total aligned reads displayed.

ACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCAC  >  NC_000913/3601356‑3601653
                                                                                                                                                     |                                                                                                                                                    
acacTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATa                                                                                                                                                      >  2:495841/1‑150 (MQ=255)
  acTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACg                                                                                                                                                    <  1:753669/150‑1 (MQ=255)
         ccGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGc                                                                                                                                             >  2:1700618/1‑150 (MQ=255)
           gATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGTCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATCGTATCACGCAg                                                                                                                                           <  2:1598046/150‑1 (MQ=255)
           gATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATCGTATCACGCAg                                                                                                                                           <  2:1726428/150‑1 (MQ=255)
            aTGAGGAGGAACACGGCAGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTAAATCTGCGTGATACGATCACGCAGc                                                                                                                                          <  1:1262415/150‑1 (MQ=255)
             tGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATCGTATCACGCAgat                                                                                                                                         <  2:1973778/150‑3 (MQ=255)
             tGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCg                                                                                                                                         <  1:2000712/150‑1 (MQ=255)
             tGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCg                                                                                                                                         <  1:415410/150‑1 (MQ=255)
             tGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCg                                                                                                                                         <  2:1545429/150‑1 (MQ=255)
                   ggAACACGGCCGCCACCATCAACACGCCGCTCATCGCCGAAACGCGCCCGACCAGACCGGTCAACGCCGCCAGACGCGCAGACCGGCTGTCAGCCATGAGCACTTCGGCAATTCCGTTAATATGCGTGATACGATCTCGCAGCGTATTCa                                                                                                                                   >  2:828445/1‑150 (MQ=255)
                   ggAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAACCCAGCTGGCATCCATCCGCCCATCGTCACTGCCGTTAATCTGCGTGATACGATCACGAAGCGTATTCa                                                                                                                                   >  2:1316347/1‑150 (MQ=255)
                   ggAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTcc                                                                                                                                   >  2:1593672/1‑149 (MQ=255)
                     aaCACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGt                                                                                                                                 <  1:692460/150‑1 (MQ=255)
                      acacGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTg                                                                                                                                <  1:1667831/150‑1 (MQ=255)
                            ccgccACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCGAGCCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCg                                                                                                                          >  1:1238328/1‑150 (MQ=255)
                              gccACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTc                                                                                                                        <  2:171150/150‑1 (MQ=255)
                                caccaTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTccc                                                                                                                      <  1:1015369/150‑1 (MQ=255)
                                  ccaTCAACACGCAGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCCGCTGTCATCCATCCGCACTTCGTCAATGGCGTTAATCTGCGTGATACGATCCCTCAGCGTATTCAGTGATTCCGTCCCCt                                                                                                                    >  2:315950/1‑150 (MQ=255)
                                   caTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTg                                                                                                                   <  2:1600454/150‑1 (MQ=255)
                                   caTCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGCGCAAACCCGCTGTCAGCCATCCGCACTTCGTCACTGCCGATAATCAGCGTGATACGATCACGCCGCGTATTCAGGGATTCCGTccccgg                                                                                                                   >  1:1886973/1‑148 (MQ=255)
                                     tCACCACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCCTCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGTTTCCGTCCCCTGGa                                                                                                                 >  2:153685/1‑150 (MQ=255)
                                       aaCACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGaaa                                                                                                               <  1:598927/150‑1 (MQ=255)
                                       aaCACGCCGATCATCGCCAAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGaaa                                                                                                               >  1:488271/1‑150 (MQ=255)
                                            gCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCACCGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCTGTGATTCCGTCCCCTGGAAATCGAg                                                                                                          >  1:580823/1‑150 (MQ=255)
                                            gCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGTGCAAACGAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAg                                                                                                          <  1:96922/150‑1 (MQ=255)
                                            gCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAg                                                                                                          <  1:1666990/150‑1 (MQ=255)
                                                   atcGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCACACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCggg                                                                                                   <  2:892229/150‑1 (MQ=255)
                                                      gCCGAAACGCGCCCGACCAGCCCGGTCAGCGCCGCCAGCCGGGCACACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATCCCATCACGCAGCGTATTCAGTGATTCCGTCCCCGGGAAATCGAGGTTCGGGATc                                                                                                >  2:982065/1‑150 (MQ=255)
                                                      gCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCCGCGTATTCAGTGATTCCGTCCCCTGGAAGTCGAGTTTCGGGCTc                                                                                                >  1:1053326/1‑150 (MQ=255)
                                                       ccGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTGCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATca                                                                                               <  1:1820616/150‑1 (MQ=255)
                                                        cGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGGCATACATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATcac                                                                                              >  1:457997/1‑150 (MQ=255)
                                                         gAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATcacc                                                                                             <  2:1426348/150‑1 (MQ=255)
                                                           aaCGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCTAACCAGCTGTCATCCATCCGCACTTCGTGAATGCCGTTAATCTGCGTGATACGAGCACGCAGCGTATTCAGTGATTCCGTCCCCTTGAAATCGAGTGTCGGGATcaccac                                                                                           >  2:135026/1‑150 (MQ=255)
                                                            aCGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCCCTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCtccacc                                                                                          >  2:1197865/1‑150 (MQ=255)
                                                                 cccGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGGATTCAGTGAGTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCAc                                                                                     >  1:948200/1‑150 (MQ=255)
                                                                    gACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGc                                                                                  <  2:888354/150‑1 (MQ=255)
                                                                    gACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGc                                                                                  <  1:488088/150‑1 (MQ=255)
                                                                    gACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGATGTCATCCATCCGCACTTCGGCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGc                                                                                  <  2:954882/150‑1 (MQ=255)
                                                                       cAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCgg                                                                               >  1:1457632/1‑150 (MQ=255)
                                                                        aGCCCGGTCAACGCCGCCAGAAGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGa                                                                              <  1:1394945/150‑1 (MQ=255)
                                                                         gCCCGGTCAACGCCGCCAGACGGGCAAACGAGCTGTCATCAATGCGCTCTTCGTCATTGCCGTTAATCTGCGTGTTACGTTCACGCAGCGTATTCAGTGAGTCCGTCCCCTGGAATTAGAGTTTCGGGATCTCCACCGCCACTGCCGGat                                                                             >  2:87244/1‑149 (MQ=255)
                                                                           ccGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGc                                                                           >  2:1159239/1‑150 (MQ=255)
                                                                           ccGGTCAACGCCGCCAGACGGGAAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCATGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGc                                                                           <  2:1496318/150‑1 (MQ=255)
                                                                                 aaCGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGttt                                                                     <  2:413129/150‑1 (MQ=255)
                                                                                 aaCGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGGTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGttt                                                                     <  2:589759/150‑1 (MQ=255)
                                                                                 aaCGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTGCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGGTTCGGGAGCACCACCGCCACTGCCGGAAGCGGGttt                                                                     <  1:1325119/150‑1 (MQ=255)
                                                                                  aCGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGtttt                                                                    <  1:1149817/150‑1 (MQ=255)
                                                                                       gccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGGGTGATACGATCACTCAGCGTATTCCGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCTCCACTGCTGGAAGCGGGTTTTCTTcc                                                               <  2:747472/150‑1 (MQ=255)
                                                                                        ccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCa                                                              <  1:803992/150‑1 (MQ=255)
                                                                                         cAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAg                                                             <  2:204006/150‑1 (MQ=255)
                                                                                          aGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGc                                                            <  1:88720/150‑1 (MQ=255)
                                                                                          aGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGc                                                            >  2:884205/1‑150 (MQ=255)
                                                                                           gACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGTGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCa                                                           >  1:842961/1‑150 (MQ=255)
                                                                                           gACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCa                                                           <  1:151032/150‑1 (MQ=255)
                                                                                           gACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCa                                                           >  1:1832491/1‑150 (MQ=255)
                                                                                             cGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCata                                                         >  2:590504/1‑150 (MQ=255)
                                                                                               ggCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGACCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATc                                                       >  1:1969913/1‑150 (MQ=255)
                                                                                                gCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCTCCACCGCCACTGCCCGAAGCGGGGTTTCTTCAAGCATATcc                                                      >  2:843266/1‑150 (MQ=255)
                                                                                                 cAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCa                                                     >  1:785557/1‑150 (MQ=255)
                                                                                                 cAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGAGATACGATCACGCAGCGTATTCAGTGATTCAGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCa                                                     >  1:1518494/1‑150 (MQ=255)
                                                                                                  aaaCCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAg                                                    <  2:1279984/150‑1 (MQ=255)
                                                                                                   aaCCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGAAACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCCCTGCCGGAAGCGGGTTTTCTTCCAGCATATCCCgc                                                   >  1:200533/1‑150 (MQ=255)
                                                                                                     ccAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAgcgc                                                 <  1:143916/150‑1 (MQ=255)
                                                                                                      cAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGca                                                >  2:1866101/1‑150 (MQ=255)
                                                                                                        gCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcacc                                              >  2:1162437/1‑150 (MQ=255)
                                                                                                         cTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcacca                                             <  1:459985/150‑1 (MQ=255)
                                                                                                          tGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcaccac                                            <  2:52631/150‑1 (MQ=255)
                                                                                                           gTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCaccacc                                           >  1:281546/1‑150 (MQ=255)
                                                                                                            tCATACATCCGCACTTCGTCAATGCAGTTAATCTGCGTGATACGATCACGAAGCGTAATCAGAGACTCCTTCCCCTCGAAATCGAGATTCGGCATCACCACCGCCACAGCCGGAAGCTGGTTTTCTTCCAGCATATCCAGCGCAccacca                                          <  1:315950/150‑1 (MQ=255)
                                                                                                              atccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCaaa                                        >  1:1983591/1‑150 (MQ=255)
                                                                                                                ccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCTGGTTTTCTTCCAGCATATCCAGCGCACCACCACAac                                      >  2:487400/1‑150 (MQ=255)
                                                                                                                ccatccGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAac                                      <  1:1700618/150‑1 (MQ=255)
                                                                                                                     cGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagac                                 >  1:1020008/1‑150 (MQ=255)
                                                                                                                     cGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagac                                 <  1:1269215/150‑1 (MQ=255)
                                                                                                                     cGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagac                                 >  1:549317/1‑150 (MQ=255)
                                                                                                                     cGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagac                                 <  1:495841/150‑1 (MQ=255)
                                                                                                                       cACTACGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTAATCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagacaa                               >  1:1870233/1‑148 (MQ=255)
                                                                                                                        aCTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCCCCGCCACTGCCGGAAGCGGGTGTGCTTCCAGCATATCCAGCGCACCACCAACaccagaccag                              >  2:119616/1‑150 (MQ=255)
                                                                                                                            cGTCAATGCCCTTAATCTGCGTGATACGATCACGCAGCGTAATCAGTGATTCCGTCCCCTGGAAATCGAGTATCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTAc                          <  1:1802034/150‑1 (MQ=255)
                                                                                                                             gTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACg                         >  1:987296/1‑150 (MQ=255)
                                                                                                                              tCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACgg                        <  2:1852480/150‑1 (MQ=255)
                                                                                                                               cAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGa                       >  2:591861/1‑150 (MQ=255)
                                                                                                                                 aTGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAAc                     >  1:1615092/1‑150 (MQ=255)
                                                                                                                                   gCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTc                   >  1:74532/1‑150 (MQ=255)
                                                                                                                                   gCCGTTAATCTGCGTGATACGAGCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTc                   <  1:1118649/150‑1 (MQ=255)
                                                                                                                                    ccGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTcc                  >  1:825032/1‑149 (MQ=255)
                                                                                                                                        tacgcTGCGTGATCGTATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCa              >  2:1664864/5‑150 (MQ=255)
                                                                                                                                        tAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCa              <  1:1297666/150‑1 (MQ=255)
                                                                                                                                         aaTCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGGTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTGTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAg             <  1:1197865/150‑1 (MQ=255)
                                                                                                                                          cgcTGCGTGATCGTATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGAAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGt            >  2:1829011/3‑150 (MQ=255)
                                                                                                                                           tCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTg           <  1:1643739/150‑1 (MQ=255)
                                                                                                                                           gcTGCGTGATCGTATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCCACGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTg           >  1:377660/2‑150 (MQ=255)
                                                                                                                                           gcTGCGTGATCGTATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTg           >  2:896301/2‑150 (MQ=255)
                                                                                                                                             tGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACGACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCg         <  2:1777090/150‑1 (MQ=255)
                                                                                                                                             tGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAACACCAGCCCAGTTACGGAACTCACCCAGTGCg         >  1:1697658/1‑150 (MQ=255)
                                                                                                                                               cgagATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTGTCGGGATCACCCCCGCCACTGACGGAGGCGGGTTTACTTCCAGCATATCAAGCGCGCCACAAAACCGAGACCAGTTTCGGAACGCAGCCACTGCGTc       >  2:73822/4‑150 (MQ=255)
                                                                                                                                                  gATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGCACTCACCCAGTGCGTCTTc    >  2:717621/1‑150 (MQ=255)
                                                                                                                                                   aTACGATCACGCAGCGTAATCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGACCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCa   <  1:989263/150‑1 (MQ=255)
                                                                                                                                                    tACGATCACACAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCAc  <  1:1945646/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACACTGTTACCGATGACGAGGAACACGGCCGCCACCATCAACACGCCGATCATCGCCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCAC  >  NC_000913/3601356‑3601653

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: