Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,194,516 A→G 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.3% 510.3 / NA 143K551K (AAA→AAG) yqiKPHB family membrane protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (80/61);  total (82/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

CGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTG  >  NC_000913/3194419‑3194611
                                                                                                 |                                                                                               
cGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCCAAGAAAAGGc                                                                                               <  2:3891677/100‑1 (MQ=255)
 gCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCa                                                                                              >  1:3882266/1‑100 (MQ=255)
  cTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCag                                                                                             <  2:1587626/100‑1 (MQ=255)
  cTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCag                                                                                             <  2:405301/100‑1 (MQ=255)
         aCTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAAc                                                                                      <  1:329462/100‑1 (MQ=255)
            ccTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                   >  2:976077/1‑100 (MQ=255)
            ccTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                   >  2:2131562/1‑100 (MQ=255)
            ccTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCct                                                                                   >  1:1799759/1‑100 (MQ=255)
               tGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCtt                                                                                <  1:1022304/100‑1 (MQ=255)
                gCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTc                                                                               >  1:4081027/1‑100 (MQ=255)
                  tCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCa                                                                             <  1:1562426/100‑1 (MQ=255)
                   naaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTcccg                                                                            >  2:1754432/2‑98 (MQ=255)
                   cAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                            >  1:114007/1‑100 (MQ=255)
                   cAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                            >  1:3933807/1‑100 (MQ=255)
                    aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                           <  1:3957246/100‑1 (MQ=255)
                    aaTGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCGTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAgg                                                                           <  2:765857/100‑1 (MQ=255)
                    aaTAAAATTGGCGTTTCGGGCGGCTCTCTGGCGGCATTGAGTTCACCCTTAACCGCACCATCTCCCGTGGAAGAAAAGGCAGAGTAAACCCTCTTCCAg                                                                            <  2:1301081/99‑1 (MQ=255)
                      tGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaa                                                                         <  2:3996670/100‑1 (MQ=255)
                      tGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaa                                                                         <  1:3552366/100‑1 (MQ=255)
                      tGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaa                                                                         <  2:2353703/100‑1 (MQ=255)
                       gAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaat                                                                        >  1:1054520/1‑100 (MQ=255)
                        aaaTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaata                                                                       >  1:2930175/1‑100 (MQ=255)
                          aTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaataat                                                                     >  2:462417/1‑100 (MQ=255)
                          aTTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGaataat                                                                     >  2:2908687/1‑100 (MQ=255)
                           ttGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATc                                                                    <  2:1275879/100‑1 (MQ=255)
                           ttGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATc                                                                    <  2:1277208/100‑1 (MQ=255)
                             ggCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATccc                                                                  <  1:1005914/100‑1 (MQ=255)
                               cGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTg                                                                >  2:2408846/1‑100 (MQ=255)
                                gTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTgg                                                               >  1:721258/1‑100 (MQ=255)
                                 tttCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGc                                                              >  2:653154/1‑100 (MQ=255)
                                  ttCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGcc                                                             <  2:1452761/100‑1 (MQ=255)
                                   tCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGccc                                                            >  1:78199/1‑100 (MQ=255)
                                   tCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGccc                                                            >  1:2689162/1‑100 (MQ=255)
                                   tCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTATTCCAGGAATAATCCCTGGccc                                                            >  2:2813205/1‑100 (MQ=255)
                                    cAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCa                                                           >  1:951818/1‑100 (MQ=255)
                                    cAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCa                                                           >  1:545696/1‑100 (MQ=255)
                                    cAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGACAGAGGAAACCCTCATCCAGGAATAATCCCTGGCCCa                                                           >  2:2814869/1‑100 (MQ=255)
                                      ggcggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCaaa                                                         >  2:3446060/1‑100 (MQ=255)
                                         ggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAg                                                      <  2:1799745/100‑1 (MQ=255)
                                         ggcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCCGGctttttggg                                                      >  2:926376/1‑92 (MQ=255)
                                           gtCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGcc                                                    >  1:2690985/2‑100 (MQ=255)
                                           cTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAgac                                                    >  2:3943822/1‑98 (MQ=255)
                                           cTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGcc                                                    <  2:664197/100‑1 (MQ=255)
                                            tCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCg                                                   <  1:2996812/100‑1 (MQ=255)
                                              aCTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCggg                                                 <  1:1304931/100‑1 (MQ=255)
                                                 ggcggcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAt                                              >  1:1014411/1‑100 (MQ=255)
                                                 ggcggcAATGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTGCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAt                                              >  2:2200289/1‑100 (MQ=255)
                                                  gcggcATTGACTTCACCCTTAACCTCAACACCTCCCGTCGCAGAAAACGCAGAGTCAACCCTCTTCCAGGACTAATCCCTGGCCCacaacccc                                                    >  2:2746593/1‑89 (MQ=255)
                                                  gcggcATTGACTTCACCCTTAACCTCAACAACGCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAtt                                             <  1:2746625/100‑1 (MQ=255)
                                                     gcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGAttttt                                          <  1:929513/100‑1 (MQ=255)
                                                     gcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCCTCCAGGAATAATACCTGCCCCAAAAAGCCGGGGAttttt                                          >  2:3535806/1‑100 (MQ=255)
                                                       aTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTta                                        >  1:1237498/1‑100 (MQ=255)
                                                         tGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTACACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTtata                                      >  2:1884212/1‑100 (MQ=255)
                                                         tGACCTCACCCTTAACCTCAACAGCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTtatn                                      <  1:1754445/100‑2 (MQ=255)
                                                          gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTatat                                     >  1:906354/1‑100 (MQ=255)
                                                          gACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTGTTTatat                                     <  1:2200309/100‑1 (MQ=255)
                                                           aCTTCACCCTTAACCTCAACAACTCCCGTCAAAGAAAAGGCAGAGTAAACCCTCTTCCAAAAATAAACACAAAAACAAAAAGCCGGGGATTTTTTATATc                                    <  2:863773/100‑1 (MQ=255)
                                                            cTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCt                                   >  1:2873664/1‑100 (MQ=255)
                                                            cTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCt                                   >  2:1369476/1‑100 (MQ=255)
                                                            cTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCt                                   >  1:3420368/1‑100 (MQ=255)
                                                             ttCAGCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCt                                   <  2:3423620/99‑1 (MQ=255)
                                                              tCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGc                                 <  1:3496404/100‑1 (MQ=255)
                                                                 cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGtt                              >  1:1913201/1‑100 (MQ=255)
                                                                 cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGtt                              >  1:1583912/1‑100 (MQ=255)
                                                                 cccTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGtt                              >  1:4098829/1‑100 (MQ=255)
                                                                 cccTTAACCTCAACAACTCCCGTCGAAAAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCACTGGCCCAAAAAACCGGGGATTTTTTATATCTGc                                 >  2:1863718/1‑97 (MQ=38)
                                                                   cTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTcc                            >  1:242231/1‑100 (MQ=255)
                                                                   cTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTcc                            <  2:1945610/100‑1 (MQ=255)
                                                                   cTTAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTcc                            >  1:542982/1‑100 (MQ=255)
                                                                     tctccTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGc                          <  2:2031341/97‑1 (MQ=255)
                                                                     tAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGc                          <  2:3166735/100‑1 (MQ=255)
                                                                     tAACCTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGc                          >  1:1103271/1‑100 (MQ=255)
                                                                        ccTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGcaga                       >  1:2031356/1‑97 (MQ=255)
                                                                         cTCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaa                      <  1:2579891/100‑1 (MQ=255)
                                                                          tCAACAACTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaa                     <  1:1869287/100‑1 (MQ=255)
                                                                           caacaaCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaa                    >  1:2412508/1‑100 (MQ=255)
                                                                           caacaaCTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTaaaaa                    >  2:3584584/1‑100 (MQ=255)
                                                                            aacaacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAg                   <  1:4111111/100‑1 (MQ=255)
                                                                            aacaacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAg                   <  2:1054515/100‑1 (MQ=255)
                                                                            aacaacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAg                   <  2:1237488/100‑1 (MQ=255)
                                                                             acaacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                  >  1:2595422/1‑100 (MQ=255)
                                                                             acaacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                  >  1:2911138/1‑100 (MQ=255)
                                                                             acaacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                  >  1:3921594/1‑100 (MQ=255)
                                                                               aacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGc                <  2:109484/100‑1 (MQ=255)
                                                                               aacTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGc                <  2:1249476/100‑1 (MQ=255)
                                                                                acTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCa               <  1:763651/100‑1 (MQ=255)
                                                                                acTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCa               <  1:2163623/100‑1 (MQ=255)
                                                                                 cTCCCGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCaa              >  1:2776114/1‑100 (MQ=255)
                                                                                   cccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCaaa             <  2:3742085/99‑1 (MQ=255)
                                                                                   cccGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt            <  1:28259/100‑1 (MQ=255)
                                                                                      gTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGAt         >  2:1707633/1‑100 (MQ=255)
                                                                                      gTCGAAGAAAAGACAGAGTAAACCCTCTTCCAGGAATAATCCCCGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGAt         >  1:3659548/1‑100 (MQ=21)
                                                                                        cGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCa       >  1:801880/1‑100 (MQ=255)
                                                                                        cGAAGAAAAGGCAGAGTAAACCCGCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCa       <  2:2978157/100‑1 (MQ=255)
                                                                                         gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg      >  2:3263934/1‑100 (MQ=255)
                                                                                         gaagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTCTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg      >  2:3466108/1‑100 (MQ=255)
                                                                                          aagaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAgg     >  2:1385800/1‑100 (MQ=255)
                                                                                           agaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgt    >  1:1486649/1‑100 (MQ=255)
                                                                                            gaaAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt   >  2:3775480/1‑100 (MQ=255)
                                                                                             aaaagGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttg  <  1:115534/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTG  >  NC_000913/3194419‑3194611

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