Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,194,516 A→G 100% K551K (AAA→AAG yqiK → PHB family membrane protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,194,5160AG99.5% 595.0 / ‑5.2 176K551K (AAA→AAGyqiKPHB family membrane protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (100/75);  total (101/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

Only 100 of 181 total aligned reads displayed.

CGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGT  >  NC_000913/3194419‑3194612
                                                                                                 |                                                                                                
cGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCGCAACAACTCCCGTCGAAGAAAAGGc                                                                                                <  1:71106/100‑1 (MQ=255)
 acTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGAGACATTGACTTCACCCTTAAACTCAACAACTACCGTCGAAGAAAAGGCa                                                                                               >  2:1566119/2‑100 (MQ=255)
  cTGATTGACTCCTTGCTCAATGAAATTGGCGTTTTAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGTAGACAAGGCag                                                                                              <  2:2055676/100‑1 (MQ=255)
  cTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCag                                                                                              >  1:1438573/1‑100 (MQ=255)
  cTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCag                                                                                              <  2:3222302/100‑1 (MQ=255)
        gACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTaaa                                                                                        >  1:3397405/1‑100 (MQ=255)
        gACTCCTTGCTCAATGAAATTGGCGTTTCAGGCCGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTaaa                                                                                        >  2:2576405/1‑100 (MQ=255)
          cTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCAGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAcc                                                                                      <  2:2754255/100‑1 (MQ=255)
           tCCTTGCTCAATGGAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAccc                                                                                     >  2:2411991/1‑100 (MQ=255)
            ccTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAcc                                                                                      >  2:1239228/1‑98 (MQ=255)
            ccTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                    <  2:1518308/100‑1 (MQ=255)
            ccTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                    >  2:1242296/1‑100 (MQ=255)
                gCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTc                                                                                <  1:3404593/100‑1 (MQ=255)
                 cTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTcc                                                                               >  1:329307/1‑100 (MQ=255)
                 cTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTcc                                                                               >  1:819527/1‑100 (MQ=255)
                 cTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTcc                                                                               >  2:2442726/1‑100 (MQ=255)
                  tCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCAGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCa                                                                              >  1:428353/1‑100 (MQ=255)
                   cAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                             >  2:839954/1‑100 (MQ=255)
                   cAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                             >  1:2727222/1‑100 (MQ=255)
                    aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACATCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                            <  1:3860595/100‑1 (MQ=255)
                    aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACACCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                            >  1:1152904/1‑100 (MQ=255)
                    aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACACCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                            <  1:3160780/100‑1 (MQ=255)
                     aTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGa                                                                           >  1:2698589/1‑100 (MQ=255)
                      tGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaa                                                                          <  2:2175688/100‑1 (MQ=255)
                       gAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaat                                                                         <  1:2592588/100‑1 (MQ=255)
                       gAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaan                                                                         <  1:3819279/100‑2 (MQ=255)
                        aaaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGTCTTCACTCTTAACCTCAACAACGCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaata                                                                        <  1:1001944/100‑1 (MQ=255)
                        aaaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaata                                                                        <  2:519654/100‑1 (MQ=255)
                           ttGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaataatn                                                                     <  2:2019157/100‑2 (MQ=255)
                              gCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATccc                                                                   <  1:1330682/99‑1 (MQ=255)
                              gCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCt                                                                  >  1:2869309/1‑100 (MQ=255)
                               cGTTTCAGGCGGCTCACTGGCTGCATTGACTTCACCCTTAACCTCAACAACTCCCGACGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTg                                                                 <  1:3556288/100‑1 (MQ=255)
                               cGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTg                                                                 <  1:2736446/100‑1 (MQ=255)
                               cGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTg                                                                 >  2:915925/1‑100 (MQ=255)
                               cGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTg                                                                 >  2:3150472/1‑100 (MQ=255)
                                gTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTgg                                                                >  2:749509/1‑100 (MQ=255)
                                gTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTgg                                                                >  1:2096923/1‑100 (MQ=255)
                                gTTTCAGGCGGCTCACGGGCGGCATTGACTTAACCCTTAACCTCATCAACTCCCGTCTAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTgg                                                                <  2:391678/100‑1 (MQ=255)
                                 tttCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGc                                                               <  2:47298/100‑1 (MQ=255)
                                   tCAGGCGGCTCACTGGCGGCATTGACTTCATCCTTATCCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGccc                                                             <  1:774985/100‑1 (MQ=255)
                                    cAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCa                                                            >  2:3897144/1‑100 (MQ=255)
                                    cAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCa                                                            >  1:1550177/1‑100 (MQ=255)
                                    cAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCa                                                            >  1:40377/1‑100 (MQ=255)
                                    cAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCa                                                            >  1:597009/1‑100 (MQ=255)
                                    cAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCa                                                            >  1:636289/1‑100 (MQ=255)
                                     aGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaa                                                           >  1:951683/1‑100 (MQ=255)
                                      ggcggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaa                                                          >  2:144388/1‑100 (MQ=255)
                                         ggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAg                                                       >  1:2259041/1‑100 (MQ=255)
                                           cTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGcc                                                     >  1:393492/1‑100 (MQ=255)
                                             cACTGGCGGCATTGACTTCACCCTTAACCTCAACAGCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGTAATAATCCCTGGCCCAAAAAGCCgg                                                   <  2:368540/100‑1 (MQ=255)
                                              aCTGGCGGCATTGACTTCACCCTTAACCTCATCATCTCCCTTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCggg                                                  <  1:2145093/100‑1 (MQ=255)
                                              aCTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCggg                                                  <  1:1937974/100‑1 (MQ=255)
                                              aCTGGCGGCATTGACTTCACCCTTAACCTAAACAACTCGCGTCGAAGAAAAGGCAGATTAAAACCTCTTCCAGAAATAATACCTAGCACAAAAAgac                                                     >  2:3848916/1‑95 (MQ=255)
                                                 ggcggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGa                                                >  1:1921887/1‑99 (MQ=255)
                                                 ggcggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAt                                               >  1:3283135/1‑100 (MQ=255)
                                                  gcggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGCAGAAAAGGCAGAGTAAACCCTCTTCAAGTAATAATCCCTGGCCCaaaaataagaaga                                                >  1:1674257/1‑90 (MQ=255)
                                                  gcggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtt                                              >  1:3910945/1‑100 (MQ=255)
                                                      tccgatcTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttt                                          <  2:647825/94‑1 (MQ=255)
                                                      cATTGACTTCACCCTTAACCTCAACAACTCCCGTCTAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttt                                          <  2:816343/100‑1 (MQ=255)
                                                      cATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttt                                          <  2:2709730/100‑1 (MQ=255)
                                                       aTTGACTTCACCCTTAACCTCAACATCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttt                                          <  1:2889748/99‑1 (MQ=255)
                                                       aTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttt                                            >  1:1784803/1‑97 (MQ=255)
                                                       aTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                         >  2:635002/1‑100 (MQ=255)
                                                         tGACTTCACCCTGAACCTCAACAACTCCCGTCGAAGAAAAGGCATAGTAAACCGTCTTCCAGGAATAATCGCTTGCCCATAAAGCCGGGGATTTTTtata                                       <  2:367226/100‑1 (MQ=255)
                                                          gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTatat                                      >  1:164272/1‑100 (MQ=255)
                                                          gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTatat                                      >  2:3495849/1‑100 (MQ=255)
                                                          gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTatat                                      >  2:2393617/1‑100 (MQ=255)
                                                              tCACCCTTAACCTCAACAACTCCAGTAGAAGAAAAGGCAGGGTAAACCCTCTTCCAGGAATAATACCTCGCCCAATAAGCCGGGGTTTTTTTATATCTga                                  <  2:1297670/100‑2 (MQ=255)
                                                               cACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCg                                 <  1:1003951/100‑1 (MQ=255)
                                                               cACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCg                                 <  2:2458491/100‑1 (MQ=255)
                                                                aCCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCg                                 <  1:3352168/99‑1 (MQ=255)
                                                                aCCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGt                                >  1:2448569/1‑100 (MQ=255)
                                                                 cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGtt                               >  1:3686845/1‑100 (MQ=255)
                                                                 cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTAATATCTGCGtt                               >  1:1243227/1‑100 (MQ=255)
                                                                   cTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTcc                             >  1:3233335/1‑100 (MQ=255)
                                                                     tAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGc                           <  2:2850186/100‑1 (MQ=255)
                                                                     tAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGc                           <  1:3632888/100‑1 (MQ=255)
                                                                        ccTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaa                        >  1:3717393/1‑100 (MQ=255)
                                                                        ccTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTAGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGATaa                        >  1:3327159/1‑100 (MQ=255)
                                                                        ccTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCACTGGCACAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaa                        >  1:1040943/1‑100 (MQ=255)
                                                                         cTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCATGGCCCAAAAAGCCGGGGGATTTTTATATCTGCGTTCCGCTaaa                       >  2:1398938/1‑100 (MQ=37)
                                                                          tCAACAAGTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaa                      >  1:2046467/1‑100 (MQ=255)
                                                                           caacaaCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaa                     >  1:3772659/1‑100 (MQ=255)
                                                                             acaacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAACAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                   >  2:2120664/1‑100 (MQ=255)
                                                                              caacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTg                  >  2:1120485/1‑100 (MQ=255)
                                                                              caacCCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGTAATAATCCCTGGCCTAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaa                     <  2:1022552/97‑1 (MQ=25)
                                                                               aacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGc                 <  2:2316046/100‑1 (MQ=255)
                                                                               aacTCCCATCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGAGTTCAAATAAAAAGTGc                 <  2:3359586/100‑1 (MQ=25)
                                                                                acTCCCGTCGAAGAAAAGGCAGAGTAAACCCGCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATTTCTGCGTTCCGCTAAAAAGTGCa                <  2:399069/100‑1 (MQ=255)
                                                                                 cTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCaa               <  2:1935529/100‑1 (MQ=255)
                                                                                   cccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt             >  2:2871323/1‑100 (MQ=255)
                                                                                        cGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCa        >  1:1798688/1‑100 (MQ=255)
                                                                                         gaagaaAAGGCCGAGTAAACCCTCTTCCAGGAATAATCCCTAGCCCAAAAAGCCGGTGATTTTTTGGATTTTCGATCCGCTAGTAGGTGCAGATGATCa        >  2:320277/1‑99 (MQ=16)
                                                                                         gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGAt          >  2:3578420/1‑97 (MQ=255)
                                                                                         gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg       <  2:126945/100‑1 (MQ=255)
                                                                                         gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg       >  2:3070980/1‑100 (MQ=255)
                                                                                          aagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAgg      >  2:3289550/1‑100 (MQ=255)
                                                                                          aagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCACTAAAAAGTGCAAATGATCAgg      <  2:3560000/100‑1 (MQ=255)
                                                                                             aaaagGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAACAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttg   <  2:3327548/100‑1 (MQ=255)
                                                                                             aaaagGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttg   <  2:3139898/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGT  >  NC_000913/3194419‑3194612

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