Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,443,086 T→A 18.8% 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,0860TA18.8% 322.1 / 65.1 138N338Y (AAC→TAC) secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (46/66);  new base A (14/12);  total (60/78)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 144 total aligned reads displayed.

CATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGT  >  NC_000913/3442988‑3443185
                                                                                                  |                                                                                                   
cATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTg                                                                                                    >  2:2625717/1‑100 (MQ=255)
cATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTg                                                                                                    <  1:1336879/100‑1 (MQ=255)
cATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTg                                                                                                    >  1:622426/1‑100 (MQ=255)
 aTTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGa                                                                                                   <  2:2806101/100‑1 (MQ=255)
 aTTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATATACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGa                                                                                                   <  1:1851981/100‑1 (MQ=255)
     cTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGaaaa                                                                                                >  2:445551/1‑99 (MQ=255)
     cTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAAc                                                                                               <  1:1803810/100‑1 (MQ=255)
     cTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAAc                                                                                               >  2:227144/1‑100 (MQ=255)
     cTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAAc                                                                                               <  1:2771698/100‑1 (MQ=255)
      tttctggATATACGTAGGCGGGTGCGCTCCTGTACGAATTCCTGGTACAAATGCACCTGACTTCGTCAGGTTATCTGCTTTTTCACGCGTGTTGAAAAc                                                                                               <  1:2610518/93‑1 (MQ=255)
       ttATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCa                                                                                             <  2:394998/100‑1 (MQ=255)
            gATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGcc                                                                                        >  2:2930874/1‑100 (MQ=255)
              tataCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCgt                                                                                      <  2:2205988/100‑1 (MQ=255)
                taCTTCGTCGGTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCgtgt                                                                                    <  2:1971729/100‑1 (MQ=255)
                taCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCgtgt                                                                                    <  2:2590873/100‑1 (MQ=255)
                taCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtgt                                                                                    >  2:2511854/1‑100 (MQ=255)
                taCTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtgt                                                                                    <  1:400610/100‑1 (MQ=255)
                  cTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTag                                                                                  >  1:1443438/1‑100 (MQ=255)
                  cTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTag                                                                                  >  1:2660835/1‑100 (MQ=255)
                  cTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCgtgt                                                                                    >  2:1325421/1‑98 (MQ=255)
                  cTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTag                                                                                  >  2:1670333/1‑100 (MQ=255)
                    tCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTagaa                                                                                <  1:2662395/100‑1 (MQ=255)
                     cGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTagaag                                                                               <  1:2449045/100‑1 (MQ=255)
                     cGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTagaag                                                                               >  1:1136890/1‑100 (MQ=255)
                      gCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTagaaga                                                                              <  1:292822/100‑1 (MQ=255)
                       ccGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAgaagaa                                                                             <  1:772934/100‑1 (MQ=255)
                         gTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAAc                                                                           >  2:1715407/1‑100 (MQ=255)
                         gTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAAc                                                                           >  2:193473/1‑100 (MQ=255)
                           ttGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACag                                                                         <  1:3059050/100‑1 (MQ=255)
                            tGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACaga                                                                        <  2:585405/100‑1 (MQ=255)
                            tGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACaga                                                                        <  2:1185187/100‑1 (MQ=255)
                             gCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACagaa                                                                       >  1:1065892/1‑100 (MQ=255)
                                ctcCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACagaagaa                                                                    >  2:1645879/1‑100 (MQ=255)
                                ctcCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACagaagaa                                                                    <  2:2034669/100‑1 (MQ=255)
                                ctcCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACagaagaa                                                                    >  2:1645589/1‑100 (MQ=255)
                                ctcCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACagaagaa                                                                    >  1:1489917/1‑100 (MQ=255)
                                  cccGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACagaagaaga                                                                  >  2:502475/1‑100 (MQ=255)
                                  cccGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACagaagaaga                                                                  <  1:448128/100‑1 (MQ=255)
                                     ggACGAATTCCTGGTACAAATGCACCGGACTTCTTCGGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAgatga                                                               >  2:2385462/1‑100 (MQ=255)
                                     ggACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAgatga                                                               <  2:1043354/100‑1 (MQ=255)
                                      gACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAgatgat                                                              <  2:1692333/100‑1 (MQ=255)
                                      gACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAgatga                                                               <  1:1018250/99‑1 (MQ=255)
                                        cGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTg                                                            <  1:3286132/100‑1 (MQ=255)
                                         gAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGc                                                           <  2:732182/100‑1 (MQ=255)
                                          aaTTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCCCGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCa                                                          >  2:2109298/1‑100 (MQ=255)
                                          aaTTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCa                                                          <  1:188891/100‑1 (MQ=255)
                                           aTTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAg                                                         <  2:2419829/100‑1 (MQ=255)
                                            ttCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGa                                                        <  2:123407/100‑1 (MQ=255)
                                             tCCTGGTACAAATGCACCGTACGTCGTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGa                                                        <  2:1464011/99‑1 (MQ=255)
                                               cTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGc                                                     >  1:344573/1‑100 (MQ=255)
                                                tGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCa                                                    >  2:1451403/1‑100 (MQ=255)
                                                tGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCa                                                    >  2:1042516/1‑100 (MQ=255)
                                                 ggTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCAt                                                   >  2:2040099/1‑100 (MQ=255)
                                                 ggTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCAt                                                   >  1:1639600/1‑100 (MQ=255)
                                                  gTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATa                                                  <  2:802403/100‑1 (MQ=255)
                                                   tACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATag                                                 <  1:2583627/100‑1 (MQ=255)
                                                   tACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATag                                                 <  1:1160956/100‑1 (MQ=255)
                                                         tGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGCAGAAACAGAAGAAGATGATTGCAGTCGCATAGAGTAac                                           <  1:2090448/100‑1 (MQ=255)
                                                         tGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGATGAATACCAACGCCGTGTAGAAGAAACGGAAGAAGATGATTGCAGACGCATAGAGTAac                                           >  1:1272064/1‑100 (MQ=255)
                                                            aCCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAacac                                         <  1:3265242/99‑1 (MQ=255)
                                                             ccGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACAt                                       <  1:2378253/100‑1 (MQ=255)
                                                               ggACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATaa                                     >  2:601713/1‑100 (MQ=255)
                                                                  cttcttCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGc                                  >  2:2422810/1‑100 (MQ=255)
                                                                  cttcttCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGc                                  >  2:300281/1‑100 (MQ=255)
                                                                  cttcttCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGc                                  <  2:1588365/100‑1 (MQ=255)
                                                                  cttcttCAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGc                                  <  2:2902503/100‑1 (MQ=255)
                                                                     cttcAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGt                               >  1:1712289/1‑100 (MQ=255)
                                                                     cttcAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGt                               >  1:873737/1‑100 (MQ=255)
                                                                     cttcAGGTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGt                               >  1:3289190/1‑100 (MQ=255)
                                                                      ttcAGGTTATCTGCTGTTTCACGCGCGCTGACAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGtt                              >  1:716240/1‑100 (MQ=255)
                                                                       tcAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTg                             <  2:2407116/100‑1 (MQ=255)
                                                                         aGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGcc                           <  1:1555888/100‑1 (MQ=255)
                                                                          ggTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGccc                          >  1:2023433/1‑100 (MQ=255)
                                                                           gTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGccga                         <  2:2289666/100‑3 (MQ=255)
                                                                           gTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCa                         <  2:1425333/100‑1 (MQ=255)
                                                                           gTTATCTGCTGTTTCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCa                         <  1:220199/100‑1 (MQ=255)
                                                                            ttATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAg                        <  1:2960232/100‑1 (MQ=255)
                                                                              aTCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGc                      <  1:3119607/100‑1 (MQ=255)
                                                                              aTCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGc                      >  1:1712468/1‑100 (MQ=255)
                                                                                ctgctgTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTg                    >  2:1509306/1‑100 (MQ=255)
                                                                                 tgctgTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGc                   >  2:767068/1‑100 (MQ=255)
                                                                                  gctgTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCa                  <  2:523468/100‑1 (MQ=255)
                                                                                    tgTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCaaa                <  1:2764483/100‑1 (MQ=255)
                                                                                       ttCACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATAc             >  2:990326/1‑100 (MQ=255)
                                                                                         cACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAg           >  1:941197/1‑100 (MQ=255)
                                                                                         cACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAg           >  1:589266/1‑100 (MQ=255)
                                                                                         cACGCGGGTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAg           <  2:2138/100‑1 (MQ=255)
                                                                                         cACGCGGGCTGAAAACCAACGCCGTGGAGAAGAAACAGAAGAAGATGATTGCACACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCCAATACAg           >  1:811727/1‑100 (MQ=255)
                                                                                         cACGCGGGATGAAACCCAACGCCGTGCAGAATAAACAGGAGTAGCTGATTGCAGACGCATAGTGTAACACATAAAAAGGTTGCCCTCGATTCAAATACa            >  1:1494763/1‑99 (MQ=255)
                                                                                          aCGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGc          <  2:1569317/100‑1 (MQ=255)
                                                                                          aCGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGc          >  2:1547442/1‑100 (MQ=255)
                                                                                          aCGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGc          <  2:3236728/100‑1 (MQ=255)
                                                                                           cgcgGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCg         >  1:1910728/1‑100 (MQ=255)
                                                                                             cgGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCGGCAAATACAGCGaa       >  1:1170791/1‑100 (MQ=255)
                                                                                              gggTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATGGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAg           <  1:3262972/95‑1 (MQ=255)
                                                                                               ggTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAAt     >  1:2217444/1‑100 (MQ=255)
                                                                                               ggTAGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAAt     >  1:208961/1‑100 (MQ=255)
                                                                                                 ttGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAAttg   <  2:1317350/100‑1 (MQ=255)
                                                                                                 ttGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAAttg   >  1:2813767/1‑100 (MQ=255)
                                                                                                  tGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAAttgt  >  1:2551526/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CATTACTTTATCGATATACTTCGCCGTTTGCTCTCCCGGACGAATTCCTGGTACAAATGCACCGGACTTCTTCAGGTTATCTGCTGTTTCACGCGGGTTGAAAACCAACGCCGTGTAGAAGAAACAGAAGAAGATGATTGCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGT  >  NC_000913/3442988‑3443185

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: