Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,177 G→T 5.1% 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 5.1% 278.9 / 4.0 78E89* (GAG→TAG) tqsApheromone AI‑2 transporter
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base T (2/2);  total (37/41)
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.78e-01

CATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCAC  >  NC_000913/1674031‑1674325
                                                                                                                                                  |                                                                                                                                                     
cATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAg                                                                                                                                                     >  2:111700/1‑149 (MQ=255)
             tGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGt                                                                                                                                       >  2:469195/1‑149 (MQ=255)
             tGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGt                                                                                                                                       >  3:466981/1‑149 (MQ=255)
               ttATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTa                                                                                                                                     >  4:252398/1‑149 (MQ=255)
               ttATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTa                                                                                                                                     >  4:557694/1‑149 (MQ=255)
                 aTTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTAcc                                                                                                                                   <  3:118486/149‑1 (MQ=255)
                     ttAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCaa                                                                                                                               >  2:77134/1‑149 (MQ=255)
                      tAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAat                                                                                                                              <  3:166798/149‑1 (MQ=255)
                       aaCCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAata                                                                                                                             <  2:531808/149‑1 (MQ=255)
                        aCCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAatat                                                                                                                            >  1:338480/1‑149 (MQ=255)
                        aCCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAatat                                                                                                                            >  2:151275/1‑149 (MQ=255)
                            gCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCa                                                                                                                        >  2:56225/1‑149 (MQ=255)
                                 tggaaCACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTAGATAACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGAGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAActct                                                                                                                   >  4:364161/4‑149 (MQ=255)
                                  gCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCta                                                                                                                  >  3:128018/1‑149 (MQ=255)
                                   cAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCtat                                                                                                                 <  4:216245/149‑1 (MQ=255)
                                     acacATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCtatta                                                                                                               <  1:56225/149‑1 (MQ=255)
                                       acaTGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATg                                                                                                             >  3:282504/1‑149 (MQ=255)
                                        caTGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGa                                                                                                            <  2:589940/149‑1 (MQ=255)
                                           ggTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgc                                                                                                         <  2:210418/149‑1 (MQ=255)
                                                       tgtgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAgcgc                                                                                             <  1:43241/149‑1 (MQ=255)
                                                             gcgTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAAc                                                                                       >  3:186702/1‑149 (MQ=255)
                                                                 gtACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgtt                                                                                   <  2:573988/149‑1 (MQ=255)
                                                                 gtACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgtt                                                                                   <  3:95813/149‑1 (MQ=255)
                                                                     tGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGttgttg                                                                               >  1:24084/1‑149 (MQ=255)
                                                                     tGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGttgttg                                                                               <  3:417632/149‑1 (MQ=255)
                                                                      ggCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGc                                                                              <  2:163553/149‑1 (MQ=255)
                                                                          gtgtCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                                >  4:513875/1‑143 (MQ=255)
                                                                          gtgtCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTAGAACCgttgtt                                                                                <  3:513875/143‑1 (MQ=255)
                                                                             tCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGt                                                                       <  3:407570/149‑1 (MQ=255)
                                                                             tCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGt                                                                       >  1:314600/1‑149 (MQ=255)
                                                                              cGATTTTTATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTa                                                                      <  4:466981/149‑1 (MQ=255)
                                                                               gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                     >  4:76432/1‑149 (MQ=255)
                                                                               gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                     >  2:520611/1‑149 (MQ=255)
                                                                               gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                     <  1:520611/149‑1 (MQ=255)
                                                                               gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                     <  3:76432/149‑1 (MQ=255)
                                                                                 ttttGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAggg                                                                   >  4:61212/1‑149 (MQ=255)
                                                                                 ttttGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAggg                                                                   >  3:106208/1‑149 (MQ=255)
                                                                                     gatgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTg                                                               <  3:61212/149‑1 (MQ=255)
                                                                                      atgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGa                                                              >  1:330524/1‑149 (MQ=255)
                                                                                       tgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATg                                                                                                             <  3:261356/101‑1 (MQ=255)
                                                                                       tgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATg                                                                                                             >  4:261356/1‑101 (MQ=255)
                                                                                         aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGt                                                           <  2:463280/149‑1 (MQ=255)
                                                                                          ccATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGtc                                                          <  4:399874/149‑1 (MQ=255)
                                                                                                tcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGtt                                                    <  4:143910/149‑1 (MQ=255)
                                                                                                     aTGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCa                                               <  1:469195/149‑1 (MQ=255)
                                                                                                      tGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAg                                              <  1:440174/149‑1 (MQ=255)
                                                                                                          gATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGAcgc                                                                                                                                             <  2:272363/50‑1 (MQ=255)
                                                                                                          gATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGAcgc                                                                                                                                             >  1:272363/1‑50 (MQ=255)
                                                                                                           aTGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGgc                                         <  4:268480/149‑1 (MQ=255)
                                                                                                               tgtTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCAt                                     <  1:111700/149‑1 (MQ=255)
                                                                                                                gtTGGTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCAtt                                    >  2:300748/1‑149 (MQ=255)
                                                                                                                     tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATAtt                               >  4:76505/1‑149 (MQ=255)
                                                                                                                     tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATAtt                               >  1:519168/1‑149 (MQ=255)
                                                                                                                     tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATAtt                               >  2:42409/1‑149 (MQ=255)
                                                                                                                       ttAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGa                             >  1:245177/1‑149 (MQ=255)
                                                                                                                        tAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGGATTATATTGAt                            <  2:79454/149‑1 (MQ=255)
                                                                                                                        tAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGAt                            <  4:478897/149‑1 (MQ=255)
                                                                                                                           cTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCg                         >  2:298502/1‑149 (MQ=255)
                                                                                                                            ttATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGa                        >  4:371329/1‑149 (MQ=255)
                                                                                                                             tATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGaa                       <  1:298502/149‑1 (MQ=255)
                                                                                                                              aTCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGttgttg                                                                               <  4:406463/92‑1 (MQ=255)
                                                                                                                              aTCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGttgttg                                                                               >  3:406463/1‑92 (MQ=255)
                                                                                                                               tCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcg                     <  1:241680/149‑1 (MQ=255)
                                                                                                                               tCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcg                     >  4:179902/1‑149 (MQ=255)
                                                                                                                                 tGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGAt                            >  3:378347/1‑140 (MQ=255)
                                                                                                                                 tGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGAt                            <  4:378347/140‑1 (MQ=255)
                                                                                                                                  gggTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcgg                  <  4:282504/149‑1 (MQ=255)
                                                                                                                                  gggTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcgg                  <  2:245177/149‑1 (MQ=255)
                                                                                                                                  gggTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcgg                  >  3:531506/1‑149 (MQ=255)
                                                                                                                                   ggTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTAGTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACgcggc                 <  2:297403/149‑1 (MQ=255)
                                                                                                                                     ttCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTACTATGACGCCGCTGCAAGCGCTTGAACCGTTGCTGCAACGCGTAGGGATTGACGTCTCAGCTGACCAGCTGGCGCATTATATCGATCCGAACGCGGCGa               <  3:364161/149‑1 (MQ=255)
                                                                                                                                     gtCCGCGTCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                                >  4:491210/2‑84 (MQ=37)
                                                                                                                                     gtCCGCGTCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                                <  3:491210/83‑1 (MQ=37)
                                                                                                                                     gtCCGCGTCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGAt                            <  1:92833/135‑1 (MQ=255)
                                                                                                                                     gtCCGCGTCAACTCGTTGAGCGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGAt                            >  2:92833/2‑136 (MQ=255)
                                                                                                                                      tCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt              <  2:24084/149‑1 (MQ=255)
                                                                                                                                        cgcgCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGa            <  1:42409/149‑1 (MQ=255)
                                                                                                                                         gcgcTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGAc           >  1:189735/1‑149 (MQ=255)
                                                                                                                                                  gAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCAc  <  1:77134/149‑1 (MQ=255)
                                                                                                                                                  gAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCAc  >  3:127954/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCTCAC  >  NC_000913/1674031‑1674325

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: