Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,178 A→C 6.2% E89A (GAG→GCG) ‡ tqsA → pheromone AI‑2 transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,674,1780AC 6.2% 286.4 / 7.1 96E89A (GAG→GCG) ‡tqsApheromone AI‑2 transporter
Reads supporting (aligned to +/- strand):  ref base A (45/45);  new base C (3/3);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCACC  >  NC_000913/1674029‑1674326
                                                                                                                                                     |                                                                                                                                                     
cGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGa                                                                                                                                                       <  2:305111/150‑1 (MQ=255)
 gCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAg                                                                                                                                                      >  1:591073/1‑150 (MQ=255)
 gCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTAGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCTCGCTCATCGat                                                                                                                                                      >  1:1700388/1‑149 (MQ=255)
  cATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGt                                                                                                                                                     <  2:1672566/150‑1 (MQ=255)
   attattTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGtt                                                                                                                                                    <  1:1205785/150‑1 (MQ=255)
       tttattTCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGAcg                                                                                                                                               <  1:1019347/150‑1 (MQ=255)
       tttattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGAcg                                                                                                                                               <  1:395512/150‑1 (MQ=255)
       tttattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGAcg                                                                                                                                               >  1:1462308/1‑150 (MQ=255)
         tattGCTGTTTTTCTTAACGCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCTTCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGTCAACTCGTTGagcgcg                                                                                                                                             <  1:793909/150‑6 (MQ=255)
         tattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGAcgcg                                                                                                                                             <  2:1511568/150‑1 (MQ=255)
         tattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGTCAACTCGTTGagcgcg                                                                                                                                             <  1:133938/150‑6 (MQ=255)
          attGCTGTTATTCTTAACCCGCTGGGGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCgg                                                                                                                                            <  1:49560/150‑1 (MQ=255)
           ttGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGa                                                                                                                                           >  2:1579502/1‑150 (MQ=255)
             gCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACg                                                                                                                                         >  1:42086/1‑150 (MQ=255)
             gCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACg                                                                                                                                         >  1:1517311/1‑150 (MQ=255)
             gCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGTCAACTCGTTGAGCGCGGaac                                                                                                                                         <  2:1250006/150‑3 (MQ=255)
             gCTGTTATTCTTAACCCGCTGGTACAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACg                                                                                                                                         >  2:1809263/1‑150 (MQ=255)
                gTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTa                                                                                                                                      >  2:28261/1‑150 (MQ=255)
                   aTTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCg                                                                                                                                   <  2:353150/150‑1 (MQ=255)
                    ttCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGc                                                                                                                                  >  1:117674/1‑150 (MQ=255)
                    ttCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGc                                                                                                                                  <  2:1080379/150‑1 (MQ=255)
                    ttCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGc                                                                                                                                  <  1:166950/150‑1 (MQ=255)
                      cTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCaa                                                                                                                                >  2:257811/1‑150 (MQ=255)
                          aCCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATc                                                                                                                            >  1:1344362/1‑150 (MQ=255)
                           cccGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCg                                                                                                                           <  2:1064101/150‑1 (MQ=255)
                           accGCTGGTGCAACAAATGGTCCGCTGTCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCg                                                                                                                           <  2:157751/149‑1 (MQ=255)
                                  gTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGGGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAActct                                                                                                                    <  1:574505/150‑1 (MQ=255)
                                       acacATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGTGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCtattat                                                                                                               >  2:899157/1‑150 (MQ=255)
                                          caTGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAc                                                                                                            >  1:1256974/1‑150 (MQ=255)
                                                 cgacggCTTGCGCCGCGTGTACTGGCGGAGGCGATTTTTATGACCATCATCGTGATGGCTATTGTTTTGCTAGTAGCTTTTCGGGGTTCCGCGCTCAAAGAGTTGACGCGGACGTTACCGCAATATCCCAACTCTATTATGACGCCGCTg                                                                                                     <  2:743873/146‑1 (MQ=255)
                                                   cTGGCGGGTGCCGCGTGTACTGGCGGTGGCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCACAATATCGCAACTCTATTATGACGCCGCTGCa                                                                                                   <  1:1611509/150‑1 (MQ=255)
                                                     ggCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAg                                                                                                 <  1:566129/150‑1 (MQ=255)
                                                         tgtgCCGCGTGTACTGTCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCt                                                                                             <  2:1288688/150‑1 (MQ=255)
                                                         tgtgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCt                                                                                             >  2:22965/1‑150 (MQ=255)
                                                          gtgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCGTCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCtt                                                                                            <  2:1282278/150‑1 (MQ=255)
                                                            gCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGa                                                                                          <  1:813814/150‑1 (MQ=255)
                                                            gCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGGGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGa                                                                                          >  2:1154710/1‑150 (MQ=255)
                                                                cgTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATAAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATAGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCg                                                                                      >  2:277956/1‑150 (MQ=255)
                                                                    tACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgt                                                                                  <  2:714586/150‑1 (MQ=255)
                                                                     aCTGGCGGTGTGGATTTTGATGACCTTCATCGTGATGGCGATGGGGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                                 <  2:1369595/150‑1 (MQ=255)
                                                                      cTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGttgttg                                                                                <  1:31197/150‑1 (MQ=255)
                                                                      cTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCAttgttg                                                                                <  2:579168/150‑1 (MQ=255)
                                                                        ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                              >  2:817582/1‑150 (MQ=255)
                                                                        ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                              <  1:857525/150‑1 (MQ=255)
                                                                         gCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCaa                                                                             <  1:1391294/150‑1 (MQ=255)
                                                                           ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAAcg                                                                           >  2:1575930/1‑150 (MQ=255)
                                                                             tgtCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAAcgcg                                                                         >  1:868352/1‑150 (MQ=255)
                                                                             tgtCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGCGTTCCGCGCTAAACGAGTTGACGCCGACGATACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAAcgcg                                                                         >  2:32400/1‑150 (MQ=255)
                                                                              gtCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGt                                                                        >  1:1201877/1‑150 (MQ=255)
                                                                                cGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                      >  2:1493307/1‑150 (MQ=255)
                                                                                cGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                      >  1:1793773/1‑150 (MQ=255)
                                                                                 gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGAGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAgg                                                                     <  1:299901/150‑1 (MQ=255)
                                                                                  aTTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAggg                                                                    >  2:49449/1‑150 (MQ=255)
                                                                                  aTTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAggg                                                                    <  1:599336/150‑1 (MQ=255)
                                                                                       gatgaCCATCATCGTGAGGGCGATGGTGGTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGa                                                               <  2:1667794/150‑1 (MQ=255)
                                                                                          gaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGATGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGt                                                            >  2:1239653/1‑150 (MQ=255)
                                                                                           aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtc                                                           <  2:42086/150‑1 (MQ=255)
                                                                                            ccATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAGCCGTTGTTGCAACGCGTAGGGATTGACGtct                                                          >  1:1259752/1‑150 (MQ=255)
                                                                                             catcatCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtctc                                                         >  2:700961/1‑150 (MQ=255)
                                                                                             catcatCGTGAGGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtctc                                                         <  1:1597732/150‑1 (MQ=255)
                                                                                              atcatcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCa                                                        <  2:191696/150‑1 (MQ=255)
                                                                                                catcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                      <  1:1497592/150‑1 (MQ=255)
                                                                                                catcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                      <  2:1517311/150‑1 (MQ=255)
                                                                                                         ggCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCt                                             >  1:1765466/1‑150 (MQ=255)
                                                                                                          gCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTg                                            <  1:1224268/150‑1 (MQ=255)
                                                                                                           cGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTgg                                           >  1:845593/1‑150 (MQ=255)
                                                                                                                 tgtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCAtt                                     <  2:1384015/150‑1 (MQ=255)
                                                                                                                  gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGATGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                    >  2:1310895/1‑150 (MQ=255)
                                                                                                                  gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGTTTGCCGTCTCAGTTGACCAGCTGGCGCATta                                    >  1:360047/1‑150 (MQ=255)
                                                                                                                  gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                    <  2:1256974/150‑1 (MQ=255)
                                                                                                                  gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                    >  2:435560/1‑150 (MQ=255)
                                                                                                                  gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGAGCATta                                    <  1:22965/150‑1 (MQ=255)
                                                                                                                   ttGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATtat                                   <  1:1646406/150‑1 (MQ=255)
                                                                                                                       tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTg                               >  1:1354180/1‑150 (MQ=255)
                                                                                                                       tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTg                               >  1:1049793/1‑150 (MQ=255)
                                                                                                                        aTTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGa                              >  2:907081/1‑150 (MQ=255)
                                                                                                                        aTTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGa                              >  1:1561269/1‑150 (MQ=255)
                                                                                                                          tAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATc                            <  2:1344362/150‑1 (MQ=255)
                                                                                                                          tAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATc                            <  1:233021/150‑1 (MQ=255)
                                                                                                                           aGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTCTTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATcc                           <  2:117674/150‑1 (MQ=255)
                                                                                                                           aGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATcc                           >  1:1150626/1‑150 (MQ=255)
                                                                                                                           aGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATcc                           >  1:944013/1‑150 (MQ=255)
                                                                                                                           aGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGGTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGCACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGGTGACCAGCTGGCGCCTTATTTTGAtgc                           >  2:303065/1‑148 (MQ=255)
                                                                                                                              ttATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGCCGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGaa                        >  1:935656/1‑150 (MQ=255)
                                                                                                                                 tCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGTAcgc                     >  1:1648971/1‑150 (MQ=255)
                                                                                                                                     acgtCCGCGTCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGACCgcggcg                 >  1:730310/4‑150 (MQ=255)
                                                                                                                                      gTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGa                <  1:257811/150‑1 (MQ=255)
                                                                                                                                       gtCCGCGTCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt               >  2:1042095/2‑150 (MQ=255)
                                                                                                                                       gtCCGCGTCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt               >  2:549137/2‑150 (MQ=255)
                                                                                                                                        tCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTCTTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGGCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATg              <  1:104878/150‑1 (MQ=255)
                                                                                                                                        tCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATg              >  2:320884/1‑150 (MQ=255)
                                                                                                                                         ccGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGATCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGa             >  1:1688578/1‑150 (MQ=255)
                                                                                                                                          cgagcTCAACGAGATGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCCTTGTTCCAACGCGTATGGATTGACGACTCAGATGACCAGCTGGCGCATTATATTGAACAGAACGCGGCGATGAc            >  2:872088/4‑150 (MQ=255)
                                                                                                                                               tgtaCGAGTTGACGCGGACGTTACCTCAATATCGCAACTCAATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGAAAAGCTTAGTGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGc       <  2:327976/147‑1 (MQ=255)
                                                                                                                                                cAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt      >  2:955609/1‑150 (MQ=255)
                                                                                                                                                  aCGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCa    <  1:700961/150‑1 (MQ=255)
                                                                                                                                                  aCGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCa    <  1:1575930/150‑1 (MQ=255)
                                                                                                                                                    gAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCACCGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGGGGCGATGACGTTGCTcccc  >  1:944096/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCACC  >  NC_000913/1674029‑1674326

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