Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,177 G→T 6.2% E89* (GAG→TAG) ‡ tqsA → pheromone AI‑2 transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,674,1770GT 6.2% 288.2 / 6.7 97E89* (GAG→TAG) ‡tqsApheromone AI‑2 transporter
Reads supporting (aligned to +/- strand):  ref base G (45/46);  new base T (3/3);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

CGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCACC  >  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)
       tttattTCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑cg                                                                                                                                               <  1:1019347/150‑1 (MQ=255)
       tttattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑cg                                                                                                                                               <  1:395512/150‑1 (MQ=255)
       tttattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑cg                                                                                                                                               >  1:1462308/1‑150 (MQ=255)
         tattGCTGTTTTTCTTAACGCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCTTCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGagcgcg                                                                                                                                             <  1:793909/150‑6 (MQ=255)
         tattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑cgcg                                                                                                                                             <  2:1511568/150‑1 (MQ=255)
         tattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGagcgcg                                                                                                                                             <  1:133938/150‑6 (MQ=255)
          attGCTGTTATTCTTAACCCGCTGGGGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCgg                                                                                                                                            <  1:49560/150‑1 (MQ=255)
           ttGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGa                                                                                                                                           >  2:1579502/1‑150 (MQ=255)
             gCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACg                                                                                                                                         >  1:42086/1‑150 (MQ=255)
             gCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACg                                                                                                                                         >  1:1517311/1‑150 (MQ=255)
             gCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGAGCGCGGaac                                                                                                                                         <  2:1250006/150‑3 (MQ=255)
             gCTGTTATTCTTAACCCGCTGGTACAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACg                                                                                                                                         >  2:1809263/1‑150 (MQ=255)
                gTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTa                                                                                                                                      >  2:28261/1‑150 (MQ=255)
                   aTTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCg                                                                                                                                   <  2:353150/150‑1 (MQ=255)
                    ttCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGc                                                                                                                                  >  1:117674/1‑150 (MQ=255)
                    ttCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGc                                                                                                                                  <  2:1080379/150‑1 (MQ=255)
                    ttCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGc                                                                                                                                  <  1:166950/150‑1 (MQ=255)
                      cTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCaa                                                                                                                                >  2:257811/1‑150 (MQ=255)
                          aCCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATc                                                                                                                            >  1:1344362/1‑150 (MQ=255)
                           cccGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCg                                                                                                                           <  2:1064101/150‑1 (MQ=255)
                           accGCTGGTGCAACAAATGGTCCGCTGTCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCg                                                                                                                           <  2:157751/149‑1 (MQ=255)
                                  gTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGGGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAActct                                                                                                                    <  1:574505/150‑1 (MQ=255)
                                       acacATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGTGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCtattat                                                                                                               >  2:899157/1‑150 (MQ=255)
                                          caTGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGAc                                                                                                            >  1:1256974/1‑150 (MQ=255)
                                                 cgacggCTTGCGCCGCGTGTACTGGCGGAGGCGATTTTTATGACCATCATCGTGATGGCTATTGTTTTGCTAGTAGCTTTTCGGGGTTCCGCGCTCAAAGAGTTGA‑CGCGGACGTTACCGCAATATCCCAACTCTATTATGACGCCGCTg                                                                                                     <  2:743873/146‑1 (MQ=255)
                                                   cTGGCGGGTGCCGCGTGTACTGGCGGTGGCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCACAATATCGCAACTCTATTATGACGCCGCTGCa                                                                                                   <  1:1611509/150‑1 (MQ=255)
                                                     ggCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAg                                                                                                 <  1:566129/150‑1 (MQ=255)
                                                         tgtgCCGCGTGTACTGTCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCt                                                                                             <  2:1288688/150‑1 (MQ=255)
                                                         tgtgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCt                                                                                             >  2:22965/1‑150 (MQ=255)
                                                          gtgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCGTCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCtt                                                                                            <  2:1282278/150‑1 (MQ=255)
                                                            gCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGa                                                                                          <  1:813814/150‑1 (MQ=255)
                                                            gCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGGGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGa                                                                                          >  2:1154710/1‑150 (MQ=255)
                                                                cgTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATAAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATAGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCg                                                                                      >  2:277956/1‑150 (MQ=255)
                                                                    tACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgt                                                                                  <  2:714586/150‑1 (MQ=255)
                                                                     aCTGGCGGTGTGGATTTTGATGACCTTCATCGTGATGGCGATGGGGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                                 <  2:1369595/150‑1 (MQ=255)
                                                                      cTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGttgttg                                                                                <  1:31197/150‑1 (MQ=255)
                                                                      cTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCAttgttg                                                                                <  2:579168/150‑1 (MQ=255)
                                                                        ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                              >  2:817582/1‑150 (MQ=255)
                                                                        ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCa                                                                              <  1:857525/150‑1 (MQ=255)
                                                                         gCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCaa                                                                             <  1:1391294/150‑1 (MQ=255)
                                                                           ggTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAAcg                                                                           >  2:1575930/1‑150 (MQ=255)
                                                                             tgtCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAAcgcg                                                                         >  1:868352/1‑150 (MQ=255)
                                                                             tgtCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGCGTTCCGCGCTAAACGAGTTGA‑CGCCGACGATACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAAcgcg                                                                         >  2:32400/1‑150 (MQ=255)
                                                                              gtCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGt                                                                        >  1:1201877/1‑150 (MQ=255)
                                                                                cGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                      >  2:1493307/1‑150 (MQ=255)
                                                                                cGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                      >  1:1793773/1‑150 (MQ=255)
                                                                                 gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGAGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAgg                                                                     <  1:299901/150‑1 (MQ=255)
                                                                                  aTTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAggg                                                                    >  2:49449/1‑150 (MQ=255)
                                                                                  aTTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAggg                                                                    <  1:599336/150‑1 (MQ=255)
                                                                                       gatgaCCATCATCGTGAGGGCGATGGTGGTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGa                                                               <  2:1667794/150‑1 (MQ=255)
                                                                                          gaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGATGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGt                                                            >  2:1239653/1‑150 (MQ=255)
                                                                                           aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtc                                                           <  2:42086/150‑1 (MQ=255)
                                                                                            ccATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAGCCGTTGTTGCAACGCGTAGGGATTGACGtct                                                          >  1:1259752/1‑150 (MQ=255)
                                                                                             catcatCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtctc                                                         >  2:700961/1‑150 (MQ=255)
                                                                                             catcatCGTGAGGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtctc                                                         <  1:1597732/150‑1 (MQ=255)
                                                                                              atcatcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCa                                                        <  2:191696/150‑1 (MQ=255)
                                                                                                catcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                      <  1:1497592/150‑1 (MQ=255)
                                                                                                catcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                      <  2:1517311/150‑1 (MQ=255)
                                                                                                         ggCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCt                                             >  1:1765466/1‑150 (MQ=255)
                                                                                                          gCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTg                                            <  1:1224268/150‑1 (MQ=255)
                                                                                                           cGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTgg                                           >  1:845593/1‑150 (MQ=255)
                                                                                                                 tgtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCAtt                                     <  2:1384015/150‑1 (MQ=255)
                                                                                                                  gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGATGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                    >  2:1310895/1‑150 (MQ=255)
                                                                                                                  gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGTTTGCCGTCTCAGTTGACCAGCTGGCGCATta                                    >  1:360047/1‑150 (MQ=255)
                                                                                                                  gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                    <  2:1256974/150‑1 (MQ=255)
                                                                                                                  gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                    >  2:435560/1‑150 (MQ=255)
                                                                                                                  gtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGAGCATta                                    <  1:22965/150‑1 (MQ=255)
                                                                                                                   ttGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATtat                                   <  1:1646406/150‑1 (MQ=255)
                                                                                                                       tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTg                               >  1:1354180/1‑150 (MQ=255)
                                                                                                                       tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTg                               >  1:1049793/1‑150 (MQ=255)
                                                                                                                        aTTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGa                              >  2:907081/1‑150 (MQ=255)
                                                                                                                        aTTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGa                              >  1:1561269/1‑150 (MQ=255)
                                                                                                                          tAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATc                            <  2:1344362/150‑1 (MQ=255)
                                                                                                                          tAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATc                            <  1:233021/150‑1 (MQ=255)
                                                                                                                           aGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTCTTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATcc                           <  2:117674/150‑1 (MQ=255)
                                                                                                                           aGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATcc                           >  1:1150626/1‑150 (MQ=255)
                                                                                                                           aGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATcc                           >  1:944013/1‑150 (MQ=255)
                                                                                                                           aGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGGTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGCACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGGTGACCAGCTGGCGCCTTATTTTGAtgc                           >  2:303065/1‑148 (MQ=255)
                                                                                                                              ttATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGCCGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGaa                        >  1:935656/1‑150 (MQ=255)
                                                                                                                                 tCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGTAcgc                     >  1:1648971/1‑150 (MQ=255)
                                                                                                                                     acgtCCGCG‑TCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGACCgcggcg                 >  1:730310/4‑150 (MQ=255)
                                                                                                                                      gTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGa                <  1:257811/150‑1 (MQ=255)
                                                                                                                                       gtCCGCG‑TCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt               >  2:1042095/2‑150 (MQ=255)
                                                                                                                                       gtCCGCG‑TCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt               >  2:549137/2‑150 (MQ=255)
                                                                                                                                        tCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTCTTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGGCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATg              <  1:104878/150‑1 (MQ=255)
                                                                                                                                        tCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATg              >  2:320884/1‑150 (MQ=255)
                                                                                                                                         ccGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGATCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGa             >  1:1688578/1‑150 (MQ=255)
                                                                                                                                          cgagcTCAACGAGATGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCCTTGTTCCAACGCGTATGGATTGACGACTCAGATGACCAGCTGGCGCATTATATTGAACAGAACGCGGCGATGAc            >  2:872088/4‑150 (MQ=255)
                                                                                                                                               tgtaCGAGTTGA‑CGCGGACGTTACCTCAATATCGCAACTCAATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGAAAAGCTTAGTGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGc       <  2:327976/147‑1 (MQ=255)
                                                                                                                                                cAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt      >  2:955609/1‑150 (MQ=255)
                                                                                                                                                  aCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCa    <  1:700961/150‑1 (MQ=255)
                                                                                                                                                  aCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCa    <  1:1575930/150‑1 (MQ=255)
                                                                                                                                                    gAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCACCGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGGGGCGATGACGTTGCTcccc  >  1:944096/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CGCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCACC  >  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: ‑