Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 418,803 A→G 5.1% I305M (ATA→ATG phoR → sensory histidine kinase PhoR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913418,8030AG 5.1% 240.4 / 2.1 78I305M (ATA→ATGphoRsensory histidine kinase PhoR
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base G (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTG  >  NC_000913/418654‑418951
                                                                                                                                                     |                                                                                                                                                    
cTGTTGAAGCAATTGCTGACGCTGTCGAAAATAGAATCCGCACCGAAGGATTTGCTCAATGAAAAGGTTGATGTGACGATGATGCTGAGCGTTGTTGAGCTAGAGGGTCAGACTATGAGTCATAAAAAACAGACATTTACCTTTGagata                                                                                                                                                      <  2:706945/150‑1 (MQ=255)
cTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGagata                                                                                                                                                      <  1:747243/150‑1 (MQ=255)
 tGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGagatag                                                                                                                                                     >  2:890374/1‑150 (MQ=255)
  ggTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGagataga                                                                                                                                                    >  2:1014167/1‑150 (MQ=255)
   gTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGagatagat                                                                                                                                                   >  1:446615/1‑150 (MQ=255)
       aGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACg                                                                                                                                               >  1:621045/1‑150 (MQ=255)
             tGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCa                                                                                                                                         <  2:959799/150‑1 (MQ=255)
             tGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCa                                                                                                                                         <  1:189117/150‑1 (MQ=255)
             tGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCa                                                                                                                                         <  2:782037/150‑1 (MQ=255)
               cTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAg                                                                                                                                       <  2:1433749/150‑1 (MQ=255)
                tGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAgg                                                                                                                                      >  1:412238/1‑150 (MQ=255)
                      tGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTg                                                                                                                                <  1:1407277/150‑1 (MQ=255)
                        tCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGTCTCAAGGTGTCTGGc                                                                                                                              <  1:1221211/150‑1 (MQ=255)
                        tCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGc                                                                                                                              <  1:582409/150‑1 (MQ=255)
                        tCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGc                                                                                                                              >  1:161120/1‑150 (MQ=255)
                         cGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCa                                                                                                                             >  2:1143298/1‑150 (MQ=255)
                           aaaaTAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAAc                                                                                                                           <  2:1476735/150‑1 (MQ=255)
                              aTAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGaa                                                                                                                        <  2:1514786/150‑1 (MQ=255)
                                    gCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAg                                                                                                                  <  1:1143298/150‑1 (MQ=255)
                                     ccGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGc                                                                                                                 >  1:1074602/1‑150 (MQ=255)
                                     ccGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGc                                                                                                                 >  2:396457/1‑150 (MQ=255)
                                               gCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGc                                                                                                       <  2:446615/150‑1 (MQ=255)
                                                    tGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGAttt                                                                                                  >  2:1261418/1‑150 (MQ=255)
                                                     gCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTc                                                                                                 >  2:185537/1‑150 (MQ=255)
                                                     gCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTc                                                                                                 >  2:191475/1‑150 (MQ=255)
                                                      cTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCg                                                                                                >  1:82596/1‑150 (MQ=255)
                                                        cAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGaa                                                                                              >  1:1233369/1‑150 (MQ=255)
                                                        cAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGaa                                                                                              >  2:110807/1‑150 (MQ=255)
                                                         aaTGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGAGCGATTTCGaa                                                                                              >  2:1459757/1‑149 (MQ=255)
                                                         aaTGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGAGCGATTTCGaa                                                                                              <  1:1459757/149‑1 (MQ=255)
                                                         aaTGAAAAGGTTGAGGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAAc                                                                                             <  1:941700/150‑1 (MQ=255)
                                                          aTGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAAcc                                                                                            >  1:369108/1‑150 (MQ=255)
                                                          aTGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAAcc                                                                                            >  1:1251761/1‑150 (MQ=255)
                                                             aaaaGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTgg                                                                                         <  2:105358/150‑1 (MQ=255)
                                                             aaaaGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTgg                                                                                         <  1:694187/150‑1 (MQ=255)
                                                               aaGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTc                                                                                       <  2:1074602/150‑1 (MQ=255)
                                                               aaGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTc                                                                                       <  1:185537/150‑1 (MQ=255)
                                                               aaGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTc                                                                                       <  2:866875/150‑1 (MQ=255)
                                                               aaGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTc                                                                                       <  1:191475/150‑1 (MQ=255)
                                                                aGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCt                                                                                      <  2:82596/150‑1 (MQ=255)
                                                                 ggTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCta                                                                                     >  2:1528946/1‑150 (MQ=255)
                                                                     gATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAAc                                                                                 <  1:672570/150‑1 (MQ=255)
                                                                       tgtgCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGTTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGc                                                                               <  2:61647/150‑1 (MQ=255)
                                                                         tgCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCg                                                                             >  1:1129241/1‑150 (MQ=255)
                                                                         tgCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGAGCGATTTCGAACCTGGTCTATAACGCCg                                                                             <  1:659832/150‑1 (MQ=255)
                                                                            cGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGa                                                                          >  2:673562/1‑150 (MQ=255)
                                                                             gatgatGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGaa                                                                         <  1:1278094/150‑1 (MQ=255)
                                                                               tgatgCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATc                                                                       >  1:717410/1‑150 (MQ=255)
                                                                                gatgCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCa                                                                      >  1:180348/1‑150 (MQ=255)
                                                                                        gcgTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCgg                                                              >  1:866997/1‑150 (MQ=255)
                                                                                        gcgTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCgg                                                              <  2:7383/150‑1 (MQ=255)
                                                                                           ttgttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAg                                                           <  1:532218/150‑1 (MQ=255)
                                                                                            tgttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAgg                                                          >  1:1014810/1‑150 (MQ=255)
                                                                                            tgttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAgg                                                          >  1:684627/1‑150 (MQ=255)
                                                                                             gttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCACCGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGc                                                         >  1:873661/1‑150 (MQ=255)
                                                                                                     cgAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCATGCatat                                                 >  2:791261/1‑150 (MQ=255)
                                                                                                     cgAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCatat                                                 >  2:856622/1‑150 (MQ=255)
                                                                                                          cTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCg                                            <  1:1014167/150‑1 (MQ=255)
                                                                                                          cTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCg                                            >  1:1102739/1‑150 (MQ=255)
                                                                                                           gcAGACTCTGAGGCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTGCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGt                                           <  2:1097080/149‑1 (MQ=255)
                                                                                                             aGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTAc                                         >  2:957996/1‑150 (MQ=255)
                                                                                                                 tctGAGTCAGAAAAAACAGACATTTACCTTTGAGATGGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTg                                     >  2:143160/1‑150 (MQ=255)
                                                                                                                 tctGAGTCAGAAAAAACAGACATTTACCTTTGAGATGGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTg                                     >  2:148944/1‑150 (MQ=255)
                                                                                                                 tctGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTg                                     >  1:725585/1‑150 (MQ=255)
                                                                                                                   tGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGc                                   <  1:1261418/150‑1 (MQ=255)
                                                                                                                   tGAGTCAGAAAAAACAGACATTTACCTATGAGATAGATAACGGCCTCAAGGTCTCTGGCAACGAAAATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGc                                   >  2:1017355/1‑150 (MQ=255)
                                                                                                                    gAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGGCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCa                                  <  2:1106266/150‑1 (MQ=255)
                                                                                                                     aGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAg                                 <  2:204406/150‑1 (MQ=255)
                                                                                                                       tCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCg                               <  2:870711/150‑1 (MQ=255)
                                                                                                                       tCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCg                               <  2:1233369/150‑1 (MQ=255)
                                                                                                                       tCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCg                               <  2:161120/150‑1 (MQ=255)
                                                                                                                        cAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGa                              >  1:177530/1‑150 (MQ=255)
                                                                                                                        cAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGa                              >  2:954743/1‑150 (MQ=255)
                                                                                                                         aGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAg                             >  2:549697/1‑150 (MQ=255)
                                                                                                                          gAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGt                            >  1:733659/1‑150 (MQ=255)
                                                                                                                                 cAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCAc                     <  2:1465932/150‑1 (MQ=255)
                                                                                                                                   gACATTTACCTTTGAGATGGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACgg                   <  1:148944/150‑1 (MQ=255)
                                                                                                                                   gACATTTACCTTTGAGATGGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACgg                   <  1:143160/150‑1 (MQ=255)
                                                                                                                                   gACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACgg                   >  2:52860/1‑150 (MQ=255)
                                                                                                                                      aTTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGc                >  1:1157246/1‑150 (MQ=255)
                                                                                                                                             tttGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAAttt         >  1:297034/1‑150 (MQ=255)
                                                                                                                                                gagaTAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGc      >  2:1529294/1‑150 (MQ=255)
                                                                                                                                                gagaTAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGc      >  2:1523637/1‑150 (MQ=255)
                                                                                                                                                  gatagataACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGt    <  2:866997/150‑1 (MQ=255)
                                                                                                                                                    tagataACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTg  <  2:1102739/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTGCGATTTCGAACCTGGTCTATAACGCCGTGAATCATACGCCGGAAGGCACGCATATCACCGTACGCTGGCAGCGAGTGCCGCACGGTGCCGAATTTAGCGTTG  >  NC_000913/418654‑418951

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