Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,321 C→T 19.6% A22V (GCT→GTT)  yffL → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3210CT19.6% 370.5 / 74.8 168A22V (GCT→GTT) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (58/77);  new base T (18/15);  total (76/92)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

Only 100 of 186 total aligned reads displayed.

TGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCA  >  NC_000913/2560222‑2560420
                                                                                                   |                                                                                                   
tGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgc                                                                                                     <  1:1543234/100‑1 (MQ=255)
tGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgc                                                                                                     <  1:19579/100‑1 (MQ=255)
 ggAATAAATTATTCAATATCAAGAACTTCAATCCATATTTGCAACGAAAGATCCAGAGTTCAAAAATAGGATCAACACAAATAACAGCACCCGGAAtgc                                                                                                     >  2:2887301/1‑99 (MQ=255)
 ggAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgct                                                                                                    >  1:1721369/1‑100 (MQ=255)
  gAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctg                                                                                                   <  1:2473736/100‑1 (MQ=255)
  gAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctg                                                                                                   >  1:3858875/1‑100 (MQ=255)
  gAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctg                                                                                                   >  2:1334503/1‑100 (MQ=255)
   aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                                  <  2:283318/100‑1 (MQ=255)
    aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                 <  2:3916092/100‑1 (MQ=255)
    aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                 <  1:2309838/100‑1 (MQ=255)
    aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                 >  2:3804101/1‑100 (MQ=255)
     tAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                                                                                                <  1:3372374/100‑1 (MQ=255)
       aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                              >  2:2890591/1‑100 (MQ=255)
       aaTTATTAAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                              <  2:2447851/100‑1 (MQ=255)
           attCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGcc                                                                                          >  1:2666679/1‑100 (MQ=255)
           attCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGcc                                                                                          >  2:1776478/1‑100 (MQ=255)
           attCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGcc                                                                                          >  1:134048/1‑100 (MQ=255)
           attCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGc                                                                                           >  2:3850128/1‑99 (MQ=255)
            ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                         >  2:918946/1‑100 (MQ=255)
            ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                         >  1:1172550/1‑100 (MQ=255)
              cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                       <  1:3322778/100‑1 (MQ=38)
              cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                       <  2:3213649/100‑1 (MQ=255)
              cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                       >  1:3131084/1‑100 (MQ=255)
              cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                       >  1:2251819/1‑100 (MQ=255)
              cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                       >  2:2003530/1‑100 (MQ=255)
              cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgag                                                                                       >  2:2540932/1‑98 (MQ=255)
              cAATAACAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCTGCATGTTGCAGCATTCCggg                                                                                       <  2:1443723/100‑1 (MQ=255)
               aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                      <  1:297711/100‑1 (MQ=38)
                atatCAAGAACTTAAATCCATGTTTTCAACGAAAGATCCAGAGTCCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                     <  2:268330/100‑1 (MQ=37)
                atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                     <  2:1561246/100‑1 (MQ=38)
                atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                     <  2:3106467/100‑1 (MQ=38)
                atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                     <  2:1966091/100‑1 (MQ=38)
                atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggggg                                                                                     <  2:4170251/100‑1 (MQ=255)
                 tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggggn                                                                                    <  1:576600/100‑2 (MQ=38)
                 tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGc                                                                                    <  2:3402653/100‑1 (MQ=38)
                 tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                    <  1:2723288/100‑1 (MQ=255)
                 tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                    >  1:701514/1‑100 (MQ=255)
                   tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                  <  1:2013486/100‑1 (MQ=255)
                   tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                  <  2:4068292/100‑1 (MQ=255)
                    cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGcttt                                                                                 <  2:1680476/100‑4 (MQ=37)
                    cAAGAACTTAAATCCATGTTTGCAACGAAAGAGCCAGAGTGCGAAAATAGGACCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                  <  1:3346412/99‑1 (MQ=255)
                     aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                                <  2:1589072/100‑1 (MQ=255)
                     aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                                <  2:3159288/100‑1 (MQ=255)
                     aaGAACTTAAATCCATGAATGCAACGAAAGATCCAGAATTCGAAAAAAGGAACAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGCCGGTa                                                                                <  2:1345058/100‑1 (MQ=255)
                      aGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTAt                                                                               <  1:820006/100‑1 (MQ=255)
                       gAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                              <  1:277904/100‑1 (MQ=255)
                         aCTTAAATCCATGTTTGCAACGAAAGATCCGGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                            <  1:3978984/100‑1 (MQ=255)
                          ctgaaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                           <  1:1896304/97‑1 (MQ=255)
                          cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                           <  2:3015210/100‑1 (MQ=255)
                           ttAAATCCATGTTTGCAACGAAAGAACCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGcttttatttgt                                                                          <  2:569061/100‑11 (MQ=21)
                                  cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCt                                                                   <  1:1314/100‑1 (MQ=255)
                                     gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                                <  2:550868/100‑1 (MQ=255)
                                       ttGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAgg                                                              >  2:2252963/1‑100 (MQ=255)
                                          cAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgaga                                                           >  2:3056617/1‑100 (MQ=255)
                                          cAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgaga                                                           >  2:1949391/1‑100 (MQ=255)
                                            aCGAAAGAGCCAGAGTAAGAAAAGAGGATCAACACAAGTAAAAGCCCCCAGAATGCTGCAACATGCCGGGGGCAGTATGAGCAACCAGCTAAAGGAGAAt                                                         <  2:1211081/100‑1 (MQ=255)
                                                aaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                     <  2:3668552/100‑1 (MQ=255)
                                                 aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                    <  2:3523741/100‑1 (MQ=255)
                                                  gATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                   <  1:2861797/100‑1 (MQ=255)
                                                   aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                                  >  2:775454/1‑100 (MQ=255)
                                                        gagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATGTCTCATGTCTg                                             <  2:2192175/100‑1 (MQ=255)
                                                           ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGata                                          <  1:3820613/100‑1 (MQ=255)
                                                            tCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGatat                                         <  1:961578/100‑1 (MQ=255)
                                                             cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCGACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATg                                        >  2:152724/1‑100 (MQ=255)
                                                               aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                      >  1:922835/1‑100 (MQ=255)
                                                                  aTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGc                                   <  2:1353024/100‑1 (MQ=255)
                                                                  aTAGGATCAACACAAATAAAAGCACCCGGAATGCTGCAACATGCCGGGGGCAGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGc                                   <  2:3962505/100‑1 (MQ=255)
                                                                   tAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCa                                  <  1:3060385/100‑1 (MQ=255)
                                                                    aGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAg                                 <  1:1410203/100‑1 (MQ=255)
                                                                     ggATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGt                                <  1:857440/100‑1 (MQ=255)
                                                                     agctCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCa                                  <  2:1640643/95‑1 (MQ=255)
                                                                       aTCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAg                              <  1:3646770/100‑1 (MQ=255)
                                                                       aTCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAg                              >  2:598078/1‑100 (MQ=255)
                                                                        tCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGa                             <  1:152482/100‑1 (MQ=255)
                                                                           acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                          >  1:2183998/1‑100 (MQ=255)
                                                                           acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGATATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                          >  2:277019/1‑100 (MQ=255)
                                                                             ttctcataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTATGATATGCTAGCAGTAGAACAAg                        >  2:2419915/11‑100 (MQ=21)
                                                                             actgaTAAAAGCCCCCGGAAGGCTGCAACATTCCTGGGGCGGTATGTGAAACAAGCGACAGTCGAATTTCGCATGTCTGCTATGCTAGCAGGAGAACaa                         >  2:1873485/5‑99 (MQ=255)
                                                                             acaAATAAAAGCCCCCGGAATGTTGCAACATGCCGGGGGCGGTATGAGTAACAAGCTAACGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                        >  2:1871148/1‑100 (MQ=255)
                                                                             acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                        <  2:2973752/100‑1 (MQ=255)
                                                                             acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                        >  1:2216984/1‑100 (MQ=255)
                                                                              caAATAAAAGCCCCGGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                       >  2:1809488/1‑100 (MQ=255)
                                                                              caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                       >  2:3236685/1‑100 (MQ=255)
                                                                              caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGAATGAGAAACAAGCTAAAGGAGAATTTCACATGTCTGATATGCTAGCAGTAGAACAAGa                       >  2:603038/1‑100 (MQ=255)
                                                                                aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                     >  1:380101/1‑100 (MQ=255)
                                                                                 aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                    <  1:1776482/100‑1 (MQ=255)
                                                                                  tAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaaca                   >  1:1650406/1‑100 (MQ=255)
                                                                                  tAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaaca                   <  1:1672787/100‑1 (MQ=255)
                                                                                      aGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac                <  1:4151171/99‑1 (MQ=255)
                                                                                       gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa              >  2:503141/1‑100 (MQ=38)
                                                                                        cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat             >  2:2698318/1‑100 (MQ=38)
                                                                                        cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat             >  1:3941369/1‑100 (MQ=38)
                                                                                         ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGACTTTCTCATGTCCGATATGCTAACAGTAGAACAAGCAACAAACAatg            >  2:881362/1‑100 (MQ=25)
                                                                                         ccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg            >  2:3459536/1‑100 (MQ=255)
                                                                                          cccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga           >  2:4115009/1‑100 (MQ=38)
                                                                                           ccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgat          >  1:1093381/1‑100 (MQ=255)
                                                                                            cGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGGCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgatg         <  1:773332/100‑1 (MQ=255)
                                                                                               aaTGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTc      <  1:2628755/100‑1 (MQ=255)
                                                                                                  gctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTc   >  1:3961108/1‑100 (MQ=255)
                                                                                                   ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCa  >  2:3563219/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCA  >  NC_000913/2560222‑2560420

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: