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,5160AG97.7% 559.7 / 0.3 168K551K (AAA→AAG) yqiKPHB family membrane protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (97/67);  total (99/69)
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.65e-01

Only 100 of 172 total aligned reads displayed.

CGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGC  >  NC_000913/3194419‑3194615
                                                                                                 |                                                                                                   
 gCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCa                                                                                                  >  1:1555027/1‑100 (MQ=255)
  cTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCag                                                                                                 <  1:2261766/100‑1 (MQ=255)
    gATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCagag                                                                                               >  2:473534/1‑100 (MQ=255)
    gATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCagag                                                                                               <  1:2674646/100‑1 (MQ=255)
     aTTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGt                                                                                              <  2:1069535/100‑1 (MQ=255)
      ttGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGCAAAGGCAGAGTa                                                                                             >  2:3121013/1‑100 (MQ=255)
        gACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAGCTCCCGTCGAAGAAAAGGCAGAGTa                                                                                             <  2:3489342/98‑1 (MQ=255)
        gACTCCTGGCGGAATGATATGGGCTTTTGAGGCGGCGCACTGGCGGCATTGATTTCACCCTTAGCATCAACAACTCCCGTCGAAGAAAAGGCAGAGt                                                                                              <  2:3544637/97‑1 (MQ=255)
         aCTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAc                                                                                          <  2:3258684/100‑1 (MQ=255)
          cTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGTCATTGACTTCACCCTTAACCTCAACTACTCCCGTCGAAGAAAAGGCAGAGTAAAc                                                                                          <  1:1315581/99‑1 (MQ=255)
          cTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAcc                                                                                         <  1:3307263/100‑1 (MQ=255)
            ccTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                       >  1:2095168/1‑100 (MQ=255)
             ctcgCTCAATGAAATTGGCTTTTCGGGCGAATCACTGGCGGCATTGAGTTCACCCTTAACCTCCACAACTCACGTCGAAGAAAAGGCAGAGTAAGccc                                                                                        <  2:3175077/95‑1 (MQ=255)
                 cTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTcc                                                                                  <  1:1764944/100‑1 (MQ=255)
                   cAATGAAATTGGCGTTTCAGGCGGCTCACTTGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                                <  1:319239/100‑1 (MQ=255)
                   cAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                                <  1:2008602/100‑1 (MQ=255)
                   cAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                                >  2:2012057/1‑100 (MQ=255)
                    aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                               >  2:736743/1‑100 (MQ=255)
                    aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                               >  1:1036805/1‑100 (MQ=255)
                    aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                               <  2:1121012/100‑1 (MQ=255)
                    aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                               <  2:809959/100‑1 (MQ=255)
                    aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                               <  2:1100542/100‑1 (MQ=255)
                      tGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaa                                                                             <  1:1943485/100‑1 (MQ=255)
                      tGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaa                                                                             >  1:3374419/1‑100 (MQ=255)
                      tGAAATTGGCGTGTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTACCCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                               <  2:3709129/98‑1 (MQ=255)
                       gAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaat                                                                            <  2:1275412/100‑1 (MQ=255)
                       gAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaat                                                                            <  1:1495389/100‑1 (MQ=255)
                       gAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaat                                                                            >  2:2983821/1‑100 (MQ=255)
                        aaaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaata                                                                           >  2:3408081/1‑100 (MQ=255)
                        aaaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaata                                                                           <  1:176923/100‑1 (MQ=255)
                        aaaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaata                                                                           >  1:1853726/1‑100 (MQ=255)
                         aaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaataa                                                                          >  1:3579177/1‑100 (MQ=255)
                          aTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaataat                                                                         >  1:3294655/1‑100 (MQ=255)
                             ggCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATccc                                                                      >  1:2373782/1‑100 (MQ=255)
                             ggCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATccc                                                                      >  1:1986236/1‑100 (MQ=255)
                              gCGTTTCAGGCGGCTCACTGGCGGCATTGCCTTCACCCTTAACCTCAACAACTCCCGTCGAAGAACAGGCAGAGTAAACCCTCTTCCGGGAATAATCCCt                                                                     >  2:60958/1‑100 (MQ=255)
                              gCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCt                                                                     <  1:1634228/100‑1 (MQ=255)
                               cGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCt                                                                     >  1:764760/1‑99 (MQ=255)
                                gTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAATTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTgg                                                                   >  2:2730218/1‑100 (MQ=255)
                                gTTTCAGGCGGCTCACTGGCGGCATTGACTTCAACCTCAACCTCAACAACTCCCGTCGAAGAAAAGGCAGTGTAAACCCTCTTCCAGGAATAATCCCTgg                                                                   >  2:2301291/1‑100 (MQ=255)
                                 tttCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACATCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGc                                                                  <  1:3593541/100‑1 (MQ=255)
                                 tttCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGc                                                                  <  1:638282/100‑1 (MQ=255)
                                  ttCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGcc                                                                 >  1:2971943/1‑100 (MQ=255)
                                   tCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGccc                                                                <  1:1403650/100‑1 (MQ=255)
                                    cAGGCGGCTCACTGGCGGGATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACGCTCTTCCAGGAATAATCCCTgg                                                                   >  2:2193042/1‑96 (MQ=255)
                                       gcggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaaa                                                            <  2:2095155/100‑1 (MQ=255)
                                        cggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaaaa                                                           >  1:2467427/1‑100 (MQ=255)
                                         ggcTCTCTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAg                                                          >  2:1227030/1‑100 (MQ=255)
                                          gcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGc                                                         >  2:70410/1‑100 (MQ=255)
                                           cTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCTAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAg                                                          <  2:729559/98‑1 (MQ=255)
                                           cTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGcc                                                        >  1:1233507/1‑100 (MQ=255)
                                             cACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCgg                                                      <  2:3563958/100‑1 (MQ=255)
                                             cACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCgg                                                      >  1:1035014/1‑100 (MQ=255)
                                              aCTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCggg                                                     >  2:1956921/1‑100 (MQ=255)
                                                 ngcggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAt                                                  >  1:3012795/2‑100 (MQ=255)
                                                 ggcggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAt                                                  >  2:3096862/1‑100 (MQ=255)
                                                    ggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttt                                               >  1:432321/1‑100 (MQ=255)
                                                    ggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttt                                               >  1:481041/1‑100 (MQ=255)
                                                       aTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttttt                                             <  1:728053/99‑1 (MQ=255)
                                                       aTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                            >  2:15178/1‑100 (MQ=255)
                                                       aTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                            >  2:1712442/1‑100 (MQ=255)
                                                        ttGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                            <  2:1139246/99‑1 (MQ=255)
                                                          gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtttggga                                            >  2:2222163/1‑93 (MQ=255)
                                                             ttCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTg                                      <  1:2703554/100‑1 (MQ=255)
                                                              tCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGc                                     <  1:1075716/100‑1 (MQ=255)
                                                              tCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGc                                     >  2:2016951/1‑100 (MQ=255)
                                                               cACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCg                                    >  1:2940368/1‑100 (MQ=255)
                                                                 cccTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGtt                                  >  1:542795/1‑100 (MQ=255)
                                                                      aaCCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGc                              <  2:2896711/99‑1 (MQ=255)
                                                                      aaCCTCAACAACTCCCGCCGAAGAAAAGGCAGATTAAACCCGCTTCCAGGAATAAAAAAAAAAAAAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGc                              <  2:1865444/99‑1 (MQ=12)
                                                                       aCCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTa                            >  1:2551383/1‑100 (MQ=255)
                                                                        ccTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaa                           >  1:739921/1‑100 (MQ=255)
                                                                        ccTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaa                           >  1:2873893/1‑100 (MQ=255)
                                                                        ccTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCC‑AAAAGCGGGGGACTTTTTATATCTGCGTTCCGCTaaa                          >  2:1409848/1‑100 (MQ=37)
                                                                         cTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaa                          >  1:370624/1‑100 (MQ=255)
                                                                           caacaaCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaa                        >  1:700380/1‑100 (MQ=255)
                                                                           caacaaCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaa                        <  2:101543/100‑1 (MQ=255)
                                                                            aacaacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAg                       >  1:936785/1‑100 (MQ=255)
                                                                             acaacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaggt                      >  2:2608925/1‑97 (MQ=255)
                                                                             acaacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                      >  1:2542861/1‑100 (MQ=255)
                                                                             acaacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                      >  2:3746157/1‑100 (MQ=255)
                                                                             acaacTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                      >  2:539496/1‑100 (MQ=255)
                                                                                 cTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCaa                  >  1:3434545/1‑100 (MQ=255)
                                                                                   cccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt                >  2:2424492/1‑100 (MQ=255)
                                                                                    ccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCTAGGAAAAAAAAAAAAAAAAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATg               <  2:922574/100‑1 (MQ=9)
                                                                                     cGTCAAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGa              <  1:3264444/100‑1 (MQ=255)
                                                                                         gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGCTTTGGAAGAGCTGCGTGCCGGGAAAAAGTGCAGATg               >  2:3799801/1‑95 (MQ=21)
                                                                                         gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg          <  1:3526101/100‑1 (MQ=255)
                                                                                         gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg          >  1:3158192/1‑100 (MQ=255)
                                                                                         gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg          >  1:2765107/1‑100 (MQ=255)
                                                                                         gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg          <  2:3031150/100‑1 (MQ=255)
                                                                                          aagaaAAGGCAGAGTAAACCCTCTTCCAGTAATAATCCCTAACCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg          <  2:3368807/99‑1 (MQ=21)
                                                                                           agaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgt        <  2:361614/100‑1 (MQ=255)
                                                                                           agaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTATATATCTGCGGTCCGCTAAAACGTGCAAATGATCAgg         >  2:1574122/1‑99 (MQ=21)
                                                                                            gaaAAGGCAGAGTAAACCCTCTTCCCGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt       >  2:2856365/1‑100 (MQ=255)
                                                                                            gaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt       >  1:3397982/1‑100 (MQ=255)
                                                                                            gaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt       >  1:1143168/1‑100 (MQ=255)
                                                                                             aaaagGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAgtt        >  2:3140599/1‑96 (MQ=255)
                                                                                             aaaagGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttg      <  1:2318479/100‑1 (MQ=255)
                                                                                                 aGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgc  <  2:542789/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
CGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGC  >  NC_000913/3194419‑3194615

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