Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,443,086 T→A 100% 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,0860TA96.0% 300.8 / 0.8 98N338Y (AAC→TAC) secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base A (46/48);  total (48/50)
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.58e-01

ATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAA  >  NC_000913/3442989‑3443181
                                                                                                 |                                                                                               
aTTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGAGGAATTCCTGTTACAAATGCACCGGACTTCTTCAGTTTATCTGCTGTTTCACTCGGGTAGa                                                                                               >  2:386330/1‑100 (MQ=255)
 ttACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGaa                                                                                              >  2:381765/1‑100 (MQ=255)
   aCTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGaaaa                                                                                            <  2:2800000/100‑1 (MQ=255)
    ctgtATCGATATACTTCGCCGTGTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAAc                                                                                           <  2:1934380/97‑1 (MQ=255)
    cTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAAc                                                                                           >  2:2444428/1‑100 (MQ=255)
    cTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAAc                                                                                           <  2:1863974/100‑1 (MQ=255)
    cTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAAc                                                                                           <  2:1866449/100‑1 (MQ=255)
         tCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACg                                                                                      <  1:17092/100‑1 (MQ=255)
              ataCTTCGCCGTTTGTTTTCTCGGATGATTTTCTGTTTCAAATGCACCGGACTTCTTCAGGTTATCTGCTGGTGCACGCGGGTAGAAAACCAACGCCgtg                                                                                 <  1:633415/100‑1 (MQ=255)
              ataCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtg                                                                                 <  2:243031/100‑1 (MQ=255)
              ataCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTACCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtg                                                                                 >  2:1552945/1‑100 (MQ=255)
               taCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtgt                                                                                >  1:2702706/1‑100 (MQ=255)
                aCTTCGCTGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTTAGGTTATCTGCTGTTTTACGCGGGTAGAAAACCAACGCCGTGTa                                                                               <  2:2868314/100‑1 (MQ=255)
                aCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTa                                                                               >  1:2020289/1‑100 (MQ=255)
                aCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTa                                                                               <  2:1270616/100‑1 (MQ=255)
                 cTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTACTGTTTCACGCGGGTAGAAAACCAACGCCGTGTag                                                                              >  1:2793588/1‑100 (MQ=255)
                  ttCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTaga                                                                             <  2:1600459/100‑1 (MQ=255)
                   tCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTagaa                                                                            >  1:2741146/1‑100 (MQ=255)
                    cgtcGTTTGCTCTACAGGACGAAAACCTGGCACAAATGCACCGGACTTCTTCAGGTTATGTGCTGTTTCAAGCGGGTAGAAAACCAATACCGTGTagaag                                                                           >  2:1585015/4‑100 (MQ=255)
                    cGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCATGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTagaan                                                                           <  1:1479155/100‑2 (MQ=255)
                    cGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCTGGTAGAAAACCAACGCCGTGTagaag                                                                           <  1:1850750/100‑1 (MQ=255)
                     gCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTagaa                                                                            >  2:820437/1‑98 (MQ=255)
                     gCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGAAGAAAACCAACGCCGTGTagaaaa                                                                          >  1:2876013/1‑98 (MQ=255)
                     gCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGATGTTTCACGCGGGTAGAAAACCAACGCCGTGTagaaga                                                                          <  2:1866818/100‑1 (MQ=255)
                        gtgtGCTATCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAAc                                                                       <  2:2759595/97‑1 (MQ=255)
                        gTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCAGCGGACTTCTTTAGGTTATCTGCTGTTTCACGCTGGTATAAATCCAACGCGGTTTAGAAGAAAc                                                                       <  2:476455/100‑1 (MQ=255)
                        gTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAAc                                                                       >  2:105659/1‑100 (MQ=255)
                         tttGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACa                                                                      <  2:1096704/100‑1 (MQ=255)
                             ctctCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaag                                                                  >  2:1432911/1‑100 (MQ=255)
                              tctcCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaaga                                                                 >  1:1694736/1‑100 (MQ=255)
                                tcCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaagaag                                                               <  1:1346522/100‑1 (MQ=255)
                                tcCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAAAAAACAAACGCCGTGTAGAAGAAACagaagaag                                                               >  1:516409/1‑100 (MQ=255)
                                  ccGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAgat                                                             <  1:123488/100‑1 (MQ=255)
                                   cGGACGAATTCCTTGTATAAATTCATCGGATTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAgat                                                             <  1:1265931/99‑1 (MQ=255)
                                   cGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAgatg                                                            <  1:505045/100‑1 (MQ=255)
                                    ggACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAgatg                                                            >  2:357955/1‑99 (MQ=255)
                                    ggACGAATTCCTGGTACAAATGCACCGACCTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaagaag                                                               >  2:2360584/1‑96 (MQ=255)
                                     gACGAATTACGAGTACAAATGCAGCGGACTTCTTCAGGTTATCTGCTATTTCACGCGGGTAGAAAACCAAATCCGTATAGAAGAAACAGAAGAAgatgag                                                          <  2:1336413/100‑2 (MQ=255)
                                      aCGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGAtt                                                         >  2:1462025/1‑100 (MQ=255)
                                       cGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTg                                                        <  2:549924/100‑1 (MQ=255)
                                        gAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCCACTCCGTGTAGAAGAAACAGAAGAAGATGATTGc                                                       <  2:200641/100‑1 (MQ=255)
                                        gAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGc                                                       <  2:2893097/100‑1 (MQ=255)
                                        gAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGc                                                       <  2:1554745/100‑1 (MQ=255)
                                        gAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGc                                                       <  2:1659176/100‑1 (MQ=255)
                                          aTTCCTGGTACAAATTCACCGGACTTATTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCTGTGTAGCAGAAACAGAAGAAGATGAtt                                                         <  2:1858521/96‑1 (MQ=255)
                                          aTTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAg                                                     >  1:1430926/1‑100 (MQ=255)
                                          aTTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGAGTAAAAGAAACAGAAGAAGATGATTGCAg                                                     >  1:1050200/1‑100 (MQ=255)
                                             ccTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACg                                                  >  2:407809/1‑100 (MQ=255)
                                             ccTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACg                                                  <  2:2470957/100‑1 (MQ=255)
                                             ccTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACg                                                  >  2:406743/1‑100 (MQ=255)
                                              cTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTACAGACg                                                  >  1:2999838/1‑99 (MQ=255)
                                               tGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCa                                                >  2:860039/1‑100 (MQ=255)
                                                 gTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGCGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATa                                              <  2:737650/100‑1 (MQ=255)
                                                  tACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATag                                             <  2:259203/100‑1 (MQ=255)
                                                   aCAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATaga                                            <  1:716746/100‑1 (MQ=255)
                                                   aCAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATaga                                            <  2:662668/100‑1 (MQ=255)
                                                    cAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATagag                                           <  2:57897/100‑1 (MQ=255)
                                                     aaaTGCACCGGACTTCTTCAGGTTATCTGTTGTTTCACGCGGGTAGAAAACCAACGCCGTGTATAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGt                                          <  1:2134380/100‑1 (MQ=255)
                                                     aaaTGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCTGACGCATAGAGt                                          <  1:1999199/100‑1 (MQ=255)
                                                     aaaTGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGGAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGt                                          <  1:1239720/100‑1 (MQ=255)
                                                      aaTTCACCGGACTTCTTCAGGTGATCTGCTGTTTCACGCGGGTAGAAATCCAACGCCGTGCAGAAGAAACAGAAGAAGATGATTGCAGACGCATagag                                           <  2:2928961/98‑1 (MQ=255)
                                                      aaTGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGAAGACGCATAGAGTa                                         >  1:2275728/1‑100 (MQ=255)
                                                      aaTGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAAAGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTa                                         >  2:2018413/1‑100 (MQ=255)
                                                        acgaCCGGGTTTCTTCAGGTGATCTGCTGTATCACGCGGGTAGACAGCCAACGCTGTGTAGAAGAAACATAAGATGATGTTTGCAGACGCATAGAGt                                          <  2:181973/94‑1 (MQ=255)
                                                          cACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATAATTGCAGACGCATAGAGTAacac                                     >  1:1229435/1‑100 (MQ=255)
                                                           aCCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAAcaca                                    <  2:1177604/100‑1 (MQ=255)
                                                            ccGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACAt                                   >  1:1156784/1‑100 (MQ=255)
                                                            ccGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGAAGCATAGAGTAACACAt                                   >  2:1775246/1‑100 (MQ=255)
                                                              ggACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATaa                                 <  1:868692/100‑1 (MQ=255)
                                                               gACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATaaa                                >  1:1541991/1‑100 (MQ=255)
                                                                 cttcttCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGc                              >  2:2378359/1‑100 (MQ=255)
                                                                    cttcAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCgg                            >  2:93242/1‑99 (MQ=255)
                                                                    cttcAGGTTATCTGCTGTTTAACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGAAACACATAAAGc                              >  2:1602255/1‑97 (MQ=255)
                                                                     ttcAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGtt                          <  1:344669/100‑1 (MQ=255)
                                                                      tcAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTg                         <  1:1221377/100‑1 (MQ=255)
                                                                       cAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGc                        <  1:2436133/100‑1 (MQ=255)
                                                                        agttgaTCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATTGAGTATCACATAAAGCGGTTGcc                       <  2:2689409/95‑1 (MQ=255)
                                                                        aGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACAAATAAAGCGGTTGcc                       >  1:2340281/1‑100 (MQ=255)
                                                                        aGGTTATCTGCTGTATCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGcc                       <  2:1722532/100‑1 (MQ=255)
                                                                          gTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCa                     <  1:1988569/100‑1 (MQ=255)
                                                                            tATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAgg                   <  1:881094/100‑1 (MQ=255)
                                                                             aTCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGc                  <  1:1717637/100‑1 (MQ=255)
                                                                               ctgctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATAGAAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTg                >  1:1811939/1‑100 (MQ=255)
                                                                                tgctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTg                <  1:1479485/99‑1 (MQ=255)
                                                                                tgctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGc               >  1:2987246/1‑100 (MQ=255)
                                                                                tgctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGc               >  1:979155/1‑100 (MQ=255)
                                                                                tgctgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGc               >  2:2425851/1‑100 (MQ=255)
                                                                                   tgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCaaa            >  2:2501502/1‑100 (MQ=255)
                                                                                   tgTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCaaa            >  1:1935629/1‑100 (MQ=255)
                                                                                    gTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAAt           >  2:932736/1‑100 (MQ=255)
                                                                                     tttCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATa          >  2:2903681/1‑100 (MQ=255)
                                                                                     tttCACGCGGGTAGAAAACCAACGCCGAGTAGAAGAAACATAAGAAGATGATTGCAGACGCAAAGAGTAACACATAAAGCGGATGCCAAGGATGAAAATa          >  2:159309/1‑100 (MQ=255)
                                                                                       tCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACa        >  2:2609175/1‑100 (MQ=255)
                                                                                       tCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGATGCCCAGGCTGCAAATACa        >  1:1406768/1‑100 (MQ=255)
                                                                                        cACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAg       >  1:1265620/1‑100 (MQ=255)
                                                                                        cACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAg       >  2:2601331/1‑100 (MQ=255)
                                                                                            cgGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGaa   >  1:2850173/1‑100 (MQ=255)
                                                                                             gggTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGaaa  <  2:2102203/100‑1 (MQ=255)
                                                                                                 |                                                                                               
ATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAA  >  NC_000913/3442989‑3443181

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