Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,340,547 A→G 36.4% E102G (GAA→GGA)  yciS → DUF1049 family inner membrane protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,340,5470AG36.4% 117.0 / 145.1 143E102G (GAA→GGA) yciSDUF1049 family inner membrane protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (39/52);  new base G (20/32);  total (59/84)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 144 total aligned reads displayed.

TCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGA  >  NC_000913/1340450‑1340646
                                                                                                 |                                                                                                   
 cccTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATa                                                                                                  >  2:3205992/1‑100 (MQ=255)
   cgggCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAAc                                                                                                <  1:3135389/98‑1 (MQ=255)
    tGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACt                                                                                               >  2:432073/1‑100 (MQ=255)
         cgcgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTa                                                                                          >  2:2152470/1‑100 (MQ=255)
           cgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATg                                                                                        >  2:2606576/1‑100 (MQ=255)
             cTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGGACTCGTCAGCGGCGAAGGGATAACTTTCTATGCt                                                                                      >  2:1478903/1‑100 (MQ=255)
              tGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTg                                                                                     <  1:2266297/100‑1 (MQ=255)
               gAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTgg                                                                                    <  1:1075043/100‑1 (MQ=255)
                  cGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAgt                                                                                 >  1:1564734/1‑100 (MQ=255)
                   gTAAAATAAAGCGACTGGGAAACCAGCGTTCACCCGCGAGTGACGTGGCTGTAGTGCGGCACCCGTAAGCGGAGAAGGAATAACTTTCTATGCTgn                                                                                    <  1:812494/96‑2 (MQ=255)
                    tAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAgttg                                                                               >  1:633379/1‑100 (MQ=255)
                    tAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttg                                                                               <  1:767764/100‑1 (MQ=255)
                      aaaTAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgtt                                                                             >  1:3469385/1‑100 (MQ=255)
                       aaTAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgttg                                                                            <  2:2158332/100‑1 (MQ=255)
                         tAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGCAGGAATAACTTTCTATGCTGGAgttgttgtt                                                                          >  1:2805936/1‑100 (MQ=255)
                         tAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAgttgttgtt                                                                          <  2:2550449/100‑1 (MQ=255)
                         tAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgttgtt                                                                          <  2:385233/100‑1 (MQ=255)
                           aaGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAGGGAATAACTTTCTATGCTGGAGTTGTTGTTTc                                                                        >  2:2756592/1‑100 (MQ=255)
                           aaGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTc                                                                        <  1:1949381/100‑1 (MQ=255)
                            aGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCt                                                                       <  1:557435/100‑1 (MQ=255)
                            aGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCt                                                                       <  2:594342/100‑1 (MQ=255)
                                 cTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCttt                                                                  <  1:195521/100‑1 (MQ=255)
                                   ggAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTg                                                                <  1:2405628/100‑1 (MQ=255)
                                   ggAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTg                                                                >  2:204747/1‑100 (MQ=255)
                                   ggAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTg                                                                >  2:197189/1‑100 (MQ=255)
                                   ggAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTg                                                                <  2:293421/100‑1 (MQ=255)
                                   ggAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTg                                                                >  1:1739549/1‑100 (MQ=255)
                                    gAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGc                                                               <  2:1036693/100‑1 (MQ=255)
                                    gAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTATGCTTTTGc                                                               >  2:948879/1‑100 (MQ=255)
                                     aaaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGcc                                                              <  2:2693776/100‑1 (MQ=255)
                                     aaaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGcc                                                              <  2:2325334/100‑1 (MQ=255)
                                     aaaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGcc                                                              <  1:676879/100‑1 (MQ=255)
                                     aaaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGcc                                                              <  2:388087/100‑1 (MQ=255)
                                      aaaCCAGCTTTCACCCGCGACTGACGTGGCTGTTGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGcc                                                              <  1:516299/99‑1 (MQ=255)
                                       aaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTg                                                            <  1:3790550/100‑1 (MQ=255)
                                        aCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGt                                                           <  2:3045627/100‑1 (MQ=255)
                                         ccAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTa                                                          >  1:1820273/1‑100 (MQ=255)
                                         ccAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTa                                                          >  2:1070661/1‑100 (MQ=255)
                                         ccAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGt                                                           >  2:4242942/1‑99 (MQ=255)
                                         ccAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTa                                                          >  2:1113671/1‑100 (MQ=255)
                                          cAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAg                                                         <  2:2481378/100‑1 (MQ=255)
                                           aGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGc                                                        >  2:3060859/1‑100 (MQ=255)
                                            gCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGcc                                                       <  2:3354004/100‑1 (MQ=255)
                                               ttCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCt                                                    <  2:4266068/100‑1 (MQ=255)
                                                tCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTg                                                   <  1:3759223/100‑1 (MQ=255)
                                                tCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTg                                                   <  2:1444060/100‑1 (MQ=255)
                                                tCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTg                                                   <  1:3111819/100‑1 (MQ=255)
                                                tCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTg                                                   >  1:3365130/1‑100 (MQ=255)
                                                tCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTg                                                   <  1:1870609/100‑1 (MQ=255)
                                                 cACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGc                                                  <  2:3225113/100‑1 (MQ=255)
                                                  aCCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGcc                                                 <  2:1694947/100‑1 (MQ=255)
                                                  aCCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTggc                                                 <  2:2887505/100‑3 (MQ=255)
                                                   cccGCGGCTGACGGGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGcc                                                 <  2:1893563/99‑1 (MQ=255)
                                                   cccGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCt                                                >  2:186217/1‑100 (MQ=255)
                                                   cccGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCt                                                >  2:57129/1‑100 (MQ=255)
                                                     cgcgACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTAt                                              <  1:2453335/100‑1 (MQ=255)
                                                      gcgACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGcc                                                 <  1:4231280/96‑1 (MQ=255)
                                                       cgACTGACGTGGCTGTAGTTCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTTGATTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATg                                             <  1:2298826/99‑1 (MQ=255)
                                                        gACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGc                                           <  1:1484000/100‑1 (MQ=255)
                                                        gACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGc                                           <  2:3963314/100‑1 (MQ=255)
                                                         aCTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCt                                          <  1:3852096/100‑1 (MQ=255)
                                                          cTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTg                                         >  1:1180736/1‑100 (MQ=255)
                                                             aCGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGta                                      <  2:1933940/100‑1 (MQ=255)
                                                               gTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGtata                                    >  1:1094498/1‑100 (MQ=255)
                                                                tGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTatat                                   <  2:2077673/100‑1 (MQ=255)
                                                                 ggCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATg                                  >  2:1569789/1‑100 (MQ=255)
                                                                  gCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATgg                                 <  2:1993109/100‑1 (MQ=255)
                                                                  gCTGTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATgg                                 <  2:2470999/100‑1 (MQ=255)
                                                                  gCTGTAGTGCCGCAATCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATg                                  >  1:1509171/1‑99 (MQ=255)
                                                                    tGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGc                               <  2:3145157/100‑1 (MQ=255)
                                                                     gTAGTGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGcc                              <  2:198879/100‑1 (MQ=255)
                                                                     gTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGcc                              <  1:3849426/100‑1 (MQ=255)
                                                                       aGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGcctc                            <  2:2903965/100‑3 (MQ=255)
                                                                         tGCCGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAg                          <  1:1040342/100‑1 (MQ=255)
                                                                          gCCGCGCTCGTCAGCGGGGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAg                          <  2:1490120/99‑1 (MQ=255)
                                                                           ccGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGaa                        >  2:4152028/1‑100 (MQ=255)
                                                                            cGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAg                       <  2:4106380/100‑1 (MQ=255)
                                                                            cGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAg                       <  2:2633450/100‑1 (MQ=255)
                                                                            cGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATCTGGGCCGCAGAAg                       >  1:834745/1‑100 (MQ=255)
                                                                            cGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAg                       >  1:3178960/1‑100 (MQ=255)
                                                                             gCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGCCTGGTATATGGGCCGCAGAAGt                      <  1:2478613/100‑1 (MQ=255)
                                                                              cACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTg                     <  2:1094501/100‑1 (MQ=255)
                                                                               aCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTgc                    <  2:4257780/100‑1 (MQ=255)
                                                                                 tCGTCAGCGGCGAAGGAATAACTTTCTATGCTGTAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTg                     <  2:871286/97‑1 (MQ=255)
                                                                                  cGTCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGca                 >  2:2621945/1‑100 (MQ=255)
                                                                                  cGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGca                 >  2:3150028/1‑100 (MQ=255)
                                                                                    tCAGCGGCGAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGcaac               <  2:2335046/100‑1 (MQ=255)
                                                                                    tCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGcaac               <  1:3063369/100‑1 (MQ=255)
                                                                                     cAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGcaaca              >  1:789640/1‑100 (MQ=255)
                                                                                       gcggcgAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGcaacaa             >  1:2657483/1‑99 (MQ=255)
                                                                                        cggcgAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCaacaaaa           <  2:473011/100‑1 (MQ=255)
                                                                                          gcgAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCaacaaaaca         >  1:920241/1‑100 (MQ=255)
                                                                                           cgAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCaacaaaacaa        >  2:4170487/1‑100 (MQ=255)
                                                                                            gAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaag       >  1:905706/1‑100 (MQ=255)
                                                                                            gAAGGGATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaag       <  2:2443959/100‑1 (MQ=255)
                                                                                            gAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCATAAGCGCGCAACAAAAcaag       <  1:1687418/100‑1 (MQ=255)
                                                                                               ggAATAACTTTCTATGCTGGGGTTGTTGTTTCTGCTGTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaagca     <  1:1234652/99‑1 (MQ=255)
                                                                                               ggAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaagcaa    >  2:1398266/1‑100 (MQ=255)
                                                                                                gAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaagcaag   <  2:1155510/100‑1 (MQ=255)
                                                                                                 aaTAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGa  <  1:379122/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
TCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGA  >  NC_000913/1340450‑1340646

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: