Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,443,086 T→A N338Y (AAC→TAC)  secY ← preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,443,0860TA93.3% 300.0 / 8.9 90N338Y (AAC→TAC) secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base A (35/49);  total (38/52)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAG  >  NC_000913/3442989‑3443176
                                                                                                 |                                                                                          
aTTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGa                                                                                          <  1:1101693/100‑1 (MQ=255)
 ttACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGaa                                                                                         >  1:400311/1‑100 (MQ=255)
   aCTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGaaaa                                                                                       >  2:1029256/1‑100 (MQ=255)
   aCTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGaaaa                                                                                       <  1:509301/100‑1 (MQ=255)
    cTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAAc                                                                                      <  1:715250/100‑1 (MQ=255)
    cTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAAc                                                                                      <  1:1558499/100‑1 (MQ=255)
      ttATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCa                                                                                    >  1:1912867/1‑100 (MQ=255)
       tATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCaa                                                                                   >  1:439293/1‑100 (MQ=255)
          cGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGc                                                                                <  1:1371282/100‑1 (MQ=255)
           gATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGcc                                                                               >  1:1171740/1‑100 (MQ=255)
           gATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCACGTTATGTGCTGTTTCACGCGGGTAGAAAACCAACGcc                                                                               >  1:1307423/1‑100 (MQ=255)
              ataCTTCGCCGTTTGGTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtg                                                                            <  2:1915882/100‑1 (MQ=255)
              ataCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtg                                                                            <  1:969311/100‑1 (MQ=255)
              ataCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtg                                                                            <  2:85259/100‑1 (MQ=255)
               taCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtgt                                                                           <  2:269254/100‑1 (MQ=255)
               taCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtgt                                                                           <  1:223046/100‑1 (MQ=255)
               taCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtgt                                                                           <  2:1369542/100‑1 (MQ=255)
                aCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTa                                                                          <  2:1009711/100‑1 (MQ=255)
                aCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTa                                                                          <  2:1168240/100‑1 (MQ=255)
                   tCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTagaa                                                                       <  1:608295/100‑1 (MQ=255)
                   tCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTagaa                                                                       >  1:285157/1‑100 (MQ=255)
                     gCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTagaaga                                                                     >  1:776923/1‑100 (MQ=255)
                     gCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTagaaga                                                                     >  2:1709223/1‑100 (MQ=255)
                     gCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTagaaga                                                                     >  2:740720/1‑100 (MQ=255)
                     gCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTagaaga                                                                     <  2:303084/100‑1 (MQ=255)
                      ccGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAgaagaa                                                                    <  2:1187749/100‑1 (MQ=255)
                      ccGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAgaagaa                                                                    <  2:664809/100‑1 (MQ=255)
                       cGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGaaa                                                                   >  2:1742316/1‑100 (MQ=255)
                        gTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAAc                                                                  <  2:1705360/100‑1 (MQ=255)
                          ttgacctcCCGGACGAATTCCTGGTACAAATGGACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAAc                                                                  <  1:1494498/93‑1 (MQ=255)
                            gCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACagaa                                                              <  2:468160/100‑1 (MQ=255)
                            gCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaa                                                              <  2:214020/100‑1 (MQ=255)
                            gCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaa                                                              >  2:1754758/1‑100 (MQ=255)
                             ctctCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaag                                                             >  1:1615458/1‑100 (MQ=255)
                             ctctCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaag                                                             >  2:1012974/1‑100 (MQ=255)
                             ctctCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaag                                                             <  1:863695/100‑1 (MQ=255)
                              tctcCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaaga                                                            <  2:186759/100‑1 (MQ=255)
                                tcCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaagaag                                                          >  2:772499/1‑100 (MQ=255)
                                 cccGGACGAATTCTGGTTAGAACTGCATCGGACTTCTTCAGGTTATTTCCTGTGTCACGGGGGTTGAAAACCAACGCCGTGTagaaga                                                                     <  1:1608712/88‑1 (MQ=255)
                                 cccGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaagaaga                                                         >  1:560081/1‑100 (MQ=255)
                                   cGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAgatg                                                       >  2:894896/1‑100 (MQ=255)
                                    ggACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAgatga                                                      <  1:313059/100‑1 (MQ=255)
                                      aCGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGAtt                                                    <  1:740719/100‑1 (MQ=255)
                                       cGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTg                                                   >  1:1491113/1‑100 (MQ=255)
                                        gAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGc                                                  >  2:1667014/1‑100 (MQ=255)
                                        gAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGc                                                  >  2:1741867/1‑100 (MQ=255)
                                         aaTTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCa                                                 <  2:873974/100‑1 (MQ=255)
                                          aTTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAg                                                >  2:493420/1‑100 (MQ=255)
                                            tCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGAc                                              >  2:74861/1‑100 (MQ=255)
                                             ccTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACg                                             <  1:748435/100‑1 (MQ=255)
                                               tGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCa                                           <  1:25458/100‑1 (MQ=255)
                                               tGGTACAAATGCACCGGACCACTGCAGGTCATCTGCTGTTTCACGCGGGTAGAAAACCAACTCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCa                                           <  1:396711/100‑1 (MQ=255)
                                                ggTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCAt                                          <  2:1029910/100‑1 (MQ=255)
                                                 gTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATa                                         >  1:1263584/1‑100 (MQ=255)
                                                  tACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATag                                        >  2:1825194/1‑100 (MQ=255)
                                                    cAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATagag                                      <  1:536907/100‑1 (MQ=255)
                                                    cAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATagag                                      <  1:1576154/100‑1 (MQ=255)
                                                     aaaTGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGt                                     <  2:521644/100‑1 (MQ=255)
                                                      aaTGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTa                                    >  1:1102065/1‑100 (MQ=255)
                                                      aaTGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTa                                    >  2:691382/1‑100 (MQ=255)
                                                        tGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAATACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAac                                  <  1:563400/100‑1 (MQ=255)
                                                        tGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAac                                  <  1:1539759/100‑1 (MQ=255)
                                                         gCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAaca                                 <  1:115241/100‑1 (MQ=255)
                                                              ggACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATaa                            <  2:303304/100‑1 (MQ=255)
                                                               gACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATaaa                           >  1:1855922/1‑100 (MQ=255)
                                                                 cttcttCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGc                         <  2:1171742/100‑1 (MQ=255)
                                                                   tcttcAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCgg                       >  1:1286927/1‑100 (MQ=255)
                                                                    cttcAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCgg                       >  1:1827915/1‑99 (MQ=255)
                                                                    cttcAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGt                      <  2:1122240/100‑1 (MQ=255)
                                                                     ttcAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTTTAGAAGAAACAGAAGTAGATGATTGCAGACGCATAGAGTATCACATAAAGCGGt                      <  1:1201340/99‑1 (MQ=255)
                                                                      tcAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTg                    <  2:31307/100‑1 (MQ=255)
                                                                      tcAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTg                    <  2:439295/100‑1 (MQ=255)
                                                                        aGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGcc                  >  2:595290/1‑100 (MQ=255)
                                                                        aGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGcc                  <  2:799212/100‑1 (MQ=255)
                                                                           cttcgatctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAg               <  1:1791833/93‑1 (MQ=255)
                                                                              tCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCt            >  2:1871216/1‑100 (MQ=255)
                                                                               ctgctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAATATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTg           >  2:1670451/1‑100 (MQ=255)
                                                                               ctgctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTg           >  1:1381800/1‑100 (MQ=255)
                                                                               ctgctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTg           >  2:19377/1‑100 (MQ=255)
                                                                                tgctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGc          <  1:772498/100‑1 (MQ=255)
                                                                                tgctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGc          <  1:172055/100‑1 (MQ=255)
                                                                                 gctgTTTCTCGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCa         <  2:317596/100‑1 (MQ=255)
                                                                                 gctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCa         >  2:261307/1‑100 (MQ=255)
                                                                                 gctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCa         <  1:1709212/100‑1 (MQ=255)
                                                                                    gTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAgagatcggaa      >  2:1791844/1‑91 (MQ=255)
                                                                                      ttCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATAc    <  2:618711/100‑1 (MQ=255)
                                                                                        cACGCGGGTAGAAAACCACCGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAg  <  1:1647704/100‑1 (MQ=255)
                                                                                        cACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAg  >  1:1791815/1‑100 (MQ=255)
                                                                                        cACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAg  <  2:1352383/100‑1 (MQ=255)
                                                                                        cACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAg  >  1:793320/1‑100 (MQ=255)
                                                                                                 |                                                                                          
ATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAG  >  NC_000913/3442989‑3443176

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