Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,955 A→T 17.5% intergenic (+8/+17) tqsA → / ← pntB pheromone AI‑2 transporter/pyridine nucleotide transhydrogenase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,674,9550AT17.5% 232.2 / 41.9 97intergenic (+8/+17)tqsA/pntBpheromone AI‑2 transporter/pyridine nucleotide transhydrogenase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (47/33);  new base T (11/6);  total (58/39)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

GCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGG  >  NC_000913/1674857‑1675052
                                                                                                  |                                                                                                 
gCCGTTGACAATTCCTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCCTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTcag                                                                                                  >  1:221434/1‑100 (MQ=255)
gCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGACTca                                                                                                   <  1:635855/99‑1 (MQ=255)
 ccGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTcagc                                                                                                 <  2:270126/100‑1 (MQ=255)
     tGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCagag                                                                                             <  2:6878/100‑1 (MQ=255)
     tGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCagag                                                                                             <  2:635885/100‑1 (MQ=255)
       aCAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGc                                                                                           <  1:908420/100‑1 (MQ=255)
        cAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGcc                                                                                          >  1:1966238/1‑100 (MQ=255)
        cAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGcc                                                                                          >  2:1384111/1‑100 (MQ=255)
        cAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGcc                                                                                          >  2:527592/1‑100 (MQ=255)
         aaTTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCg                                                                                         <  2:465417/100‑1 (MQ=255)
           ttattGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCGt                                                                                        <  1:283323/99‑1 (MQ=255)
           ttattGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCGTc                                                                                       <  1:658458/100‑1 (MQ=255)
            tattGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCa                                                                                      <  1:1455294/100‑1 (MQ=255)
             attGTCAAGATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCGTCac                                                                                     <  1:2359502/100‑2 (MQ=255)
             attGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCGTCac                                                                                     <  2:2113911/100‑2 (MQ=255)
             attGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAg                                                                                     <  2:876403/100‑1 (MQ=255)
               tGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCGTCactc                                                                                   <  2:1999789/100‑4 (MQ=39)
               tGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAggg                                                                                   <  2:1984991/100‑1 (MQ=255)
                gTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAgcg                                                                                   >  1:1709162/1‑97 (MQ=255)
                gTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGt                                                                                  <  1:2019947/100‑1 (MQ=255)
                gTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGt                                                                                  <  1:446950/100‑1 (MQ=255)
                gTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGt                                                                                  >  2:972883/1‑100 (MQ=255)
                gTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGt                                                                                  >  1:1230654/1‑100 (MQ=255)
                  cAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTa                                                                                <  1:458597/100‑1 (MQ=255)
                          cgcTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCtt                                                                        >  2:1117072/1‑100 (MQ=255)
                          cgcTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCtt                                                                        >  1:1303689/1‑100 (MQ=255)
                          cgcTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCtt                                                                        >  1:498769/1‑100 (MQ=255)
                            cTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTc                                                                      >  2:842846/1‑100 (MQ=255)
                             ttGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCa                                                                     <  1:320307/100‑1 (MQ=255)
                             ttGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTCTCa                                                                     >  1:252277/1‑100 (MQ=255)
                              tGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCan                                                                    <  2:1787312/100‑2 (MQ=255)
                              tGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAg                                                                    >  2:495175/1‑100 (MQ=255)
                               gAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAgg                                                                   <  2:1706785/100‑1 (MQ=255)
                                 acaaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGAt                                                                 >  2:1709666/1‑100 (MQ=255)
                                 acaaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGAt                                                                 >  1:1045010/1‑100 (MQ=255)
                                  caaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGAtt                                                                >  1:1929643/1‑100 (MQ=255)
                                  caaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGAtt                                                                <  1:2006121/100‑1 (MQ=255)
                                   aaccaGCGGGTGCTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTg                                                               <  1:1645582/97‑1 (MQ=255)
                                   aaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTg                                                               <  2:682039/100‑1 (MQ=255)
                                    aacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGc                                                              >  2:2321124/1‑100 (MQ=255)
                                      caGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCAt                                                            >  1:2557010/1‑100 (MQ=255)
                                         cGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCa                                                         <  2:2425932/100‑1 (MQ=255)
                                           ggAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACg                                                       <  1:2321110/100‑1 (MQ=255)
                                              ggTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTg                                                    >  2:661625/1‑100 (MQ=255)
                                              ggTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACAGCCTCAGCAGAGGCCGTCAGGGTTACAGCGCTTTCAGGATTGCATCCAAGCTg                                                    >  1:297035/1‑100 (MQ=255)
                                               gTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTgg                                                   <  2:1890752/100‑1 (MQ=255)
                                               cttccgatctCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGTGTTACAGAGCTTTCAGGATTGCATCCACGCTgg                                                   <  1:1703023/91‑1 (MQ=255)
                                                tCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGc                                                  >  1:1806176/1‑100 (MQ=255)
                                                tCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGc                                                  >  2:803527/1‑100 (MQ=255)
                                                tCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGc                                                  >  2:422009/1‑100 (MQ=255)
                                                tCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCATCAGCAGAGGCCGCCAGAGTTACAGAGCTTTCAGGATTACATcc                                                          >  1:2021605/1‑92 (MQ=255)
                                                 cAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCt                                                 >  1:1688355/1‑100 (MQ=255)
                                                    aGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTg                                              <  1:1372787/100‑1 (MQ=255)
                                                      cATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGc                                            >  1:2202502/1‑100 (MQ=255)
                                                       aTCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCg                                           >  2:1177123/1‑100 (MQ=255)
                                                         cGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTggagatcggaag                                         >  2:1703032/1‑90 (MQ=255)
                                                         cGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTc                                         >  1:107312/1‑100 (MQ=255)
                                                          gCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCa                                        >  2:2203572/1‑100 (MQ=255)
                                                           ccGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCAc                                       <  2:163940/100‑1 (MQ=255)
                                                           ccGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCAc                                       <  1:2579741/100‑1 (MQ=255)
                                                            cGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCAcc                                      <  1:1298445/100‑1 (MQ=255)
                                                               tCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCaaa                                   >  1:2451052/1‑100 (MQ=255)
                                                               tCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCaa                                    <  1:731084/99‑1 (MQ=255)
                                                                cTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAAc                                  <  2:1953823/100‑1 (MQ=255)
                                                                 tGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACa                                 >  2:1728063/1‑100 (MQ=255)
                                                                  gTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCCCCAAACAg                                >  2:1624823/1‑100 (MQ=255)
                                                                    tAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCa                              >  2:2369524/1‑100 (MQ=255)
                                                                    tAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCa                              <  2:793674/100‑1 (MQ=255)
                                                                    tAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCa                              >  1:669642/1‑100 (MQ=255)
                                                                    tAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCa                              <  2:107313/100‑1 (MQ=255)
                                                                    tAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCa                              <  2:1929653/100‑1 (MQ=255)
                                                                    tAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGACGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCa                              >  1:274140/1‑100 (MQ=255)
                                                                     aaGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCAt                             >  2:153467/1‑100 (MQ=255)
                                                                            ctcAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGgtg                      >  1:169251/1‑100 (MQ=255)
                                                                             tcAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGgtgt                     >  1:2061642/1‑100 (MQ=255)
                                                                             tcAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGgtgt                     >  2:740658/1‑100 (MQ=255)
                                                                               aaTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGttt                   >  2:1922477/1‑100 (MQ=255)
                                                                                aTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGtttt                  <  1:1060969/100‑1 (MQ=255)
                                                                                aTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGtttt                  <  2:1335296/100‑1 (MQ=255)
                                                                                   aaccctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCGGCGCTGGGTTTGGCGTGACCAAACAGCATGTGGGTGTTTTctt               >  1:150014/6‑98 (MQ=37)
                                                                                   aaGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCt               >  1:337430/1‑100 (MQ=255)
                                                                                    agagTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCtt              <  2:669645/100‑1 (MQ=255)
                                                                                     ccctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTg             >  1:1766829/4‑100 (MQ=39)
                                                                                      cctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTga            >  2:402510/3‑100 (MQ=255)
                                                                                      cctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTga            >  2:1979784/3‑100 (MQ=255)
                                                                                      cctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTga            >  2:2327613/3‑100 (MQ=255)
                                                                                      cctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTg             >  1:1572821/3‑99 (MQ=255)
                                                                                      cctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGCGTTTTCCTTga            >  2:2216129/3‑100 (MQ=39)
                                                                                       gTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaa           <  1:1177119/100‑1 (MQ=255)
                                                                                       ctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaa           >  2:1246824/2‑100 (MQ=255)
                                                                                       ctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaa           >  2:1960915/2‑100 (MQ=255)
                                                                                       ctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaa           >  2:293460/2‑100 (MQ=255)
                                                                                        tGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaag          >  2:2164046/1‑100 (MQ=255)
                                                                                        tGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaag          <  2:1966248/100‑1 (MQ=255)
                                                                                        tGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaag          >  2:718058/1‑100 (MQ=255)
                                                                                         gACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaaga         >  1:2457233/1‑100 (MQ=255)
                                                                                           cGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAAc       >  2:1230313/1‑100 (MQ=255)
                                                                                            ggCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACa      >  1:2550985/1‑100 (MQ=255)
                                                                                             gCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAg     >  2:2443302/1‑100 (MQ=255)
                                                                                                tCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCgg  >  2:2240209/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
GCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGG  >  NC_000913/1674857‑1675052

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