Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,321 C→T 19.0% 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.0% 418.1 / 81.6 188A22V (GCT→GTT) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (65/87);  new base T (14/22);  total (79/109)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-02

Only 100 of 199 total aligned reads displayed.

TGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCA  >  NC_000913/2560222‑2560420
                                                                                                   |                                                                                                   
tGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgc                                                                                                     >  2:2112794/1‑100 (MQ=255)
  gAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctg                                                                                                   <  1:2019984/100‑1 (MQ=255)
   aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                                  >  1:573495/1‑100 (MQ=255)
    aTAAATTATTCAATATCAAGAACTTAAATCCATTTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                 <  2:894430/100‑1 (MQ=255)
    aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                 >  2:692337/1‑100 (MQ=255)
    aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                 <  1:3602104/100‑1 (MQ=255)
    aTAAATTAATCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGAACCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                 <  1:1097640/100‑1 (MQ=255)
      aaaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAAc                                                                                               >  2:2030439/1‑100 (MQ=255)
       aaTTATTCAATATCAAGAACTTGAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                              <  2:2645846/100‑1 (MQ=255)
       aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                              <  2:3829360/100‑1 (MQ=255)
       aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                              >  1:2877187/1‑100 (MQ=255)
       aaTTATTCAATATAAAGAACTTAAATACATGTTTGCAACGAAAGATACAGAGTTCAAAAATAGGATAAAAACAAATAAAAGCGAGCGGAAtgcaaaaaca                                                                                              >  1:4026738/1‑93 (MQ=255)
        attattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                              <  1:893461/99‑1 (MQ=255)
            ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                         >  2:2566827/1‑100 (MQ=255)
               atcttCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                      <  2:4110383/96‑1 (MQ=38)
               aaTATCAAGAACTTAAATTTATTTTTTCTACTAATGATCCTGCGTTGGCAAATAAGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                       <  1:3900360/99‑1 (MQ=12)
                atatCAAGAACTTAAATCCATGTTTTCAACGAAAGATCCAGAGTTCGAAAATAGGATCACCACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                     <  1:2331877/100‑1 (MQ=37)
                atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                     <  2:3242502/100‑1 (MQ=38)
                atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                     <  2:1450202/100‑1 (MQ=38)
                atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                     <  2:1046285/100‑1 (MQ=38)
                atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                     <  2:2738834/100‑1 (MQ=38)
                atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggggg                                                                                     >  1:3995884/1‑100 (MQ=255)
                 tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGc                                                                                    <  2:4305471/100‑1 (MQ=38)
                 tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGc                                                                                    <  1:1314270/100‑1 (MQ=38)
                  ctCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                                                                                   <  2:444231/99‑1 (MQ=255)
                  atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGct                                                                                   <  2:750113/100‑2 (MQ=38)
                  atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGc                                                                                    <  1:4265533/99‑1 (MQ=38)
                  atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCGGCATGTTGCAGCATTCCGGGGGctt                                                                                  <  2:1630198/100‑3 (MQ=37)
                   tCAAGAACTTAAATCCATGTTTGCAACGCAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                  >  1:2844282/1‑100 (MQ=255)
                   tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                  <  2:4158322/100‑1 (MQ=255)
                   tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                  >  1:2490302/1‑100 (MQ=255)
                   tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGCTAGAAAATAGGATCAACACAAGTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                  <  1:1052170/100‑1 (MQ=255)
                    cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCGCCGGAATGCTGCAACATGCCgggg                                                                                      >  1:3734700/1‑95 (MQ=255)
                     aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                                >  2:2934018/1‑100 (MQ=255)
                       gAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTAtn                                                                              <  1:164624/100‑2 (MQ=255)
                       gAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                              <  2:830915/100‑1 (MQ=255)
                        aaCTTAAATCCCTGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAATAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                              <  1:3717810/99‑1 (MQ=255)
                        aaCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                             <  1:933251/100‑1 (MQ=255)
                         aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                            <  2:2825725/100‑1 (MQ=255)
                         aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                             <  1:1842651/99‑1 (MQ=255)
                               aTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                      <  2:2391224/100‑1 (MQ=255)
                               aTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                      <  2:545150/100‑1 (MQ=255)
                               aTCCAGGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACan                                                                      <  2:1669962/100‑2 (MQ=255)
                                tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAg                                                                     <  1:4054315/100‑1 (MQ=255)
                                 ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                    <  1:4333270/100‑1 (MQ=255)
                                 ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                    <  1:2699958/100‑1 (MQ=255)
                                  cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCt                                                                   >  1:3347186/1‑100 (MQ=255)
                                         gCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATTCCGGGGGCGGTATGAGAAACAAGCTAAAGgag                                                            <  2:4179326/100‑1 (MQ=255)
                                         gCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgag                                                            >  2:3014572/1‑100 (MQ=255)
                                         gCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgag                                                            <  1:752098/100‑1 (MQ=255)
                                          cAACGAACGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgag                                                            >  1:2302400/1‑99 (MQ=255)
                                          cAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgaga                                                           >  1:1431994/1‑100 (MQ=255)
                                                aaaaTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                     <  1:3068107/97‑1 (MQ=255)
                                                aaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAGGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                     <  1:873371/100‑1 (MQ=255)
                                                 aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                    <  2:3868355/100‑1 (MQ=255)
                                                 aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                    <  2:2525665/100‑1 (MQ=255)
                                                 aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                    <  2:2186050/100‑1 (MQ=255)
                                                      cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCTGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                               <  1:2074859/100‑1 (MQ=255)
                                                      cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                               >  2:374358/1‑100 (MQ=255)
                                                       agagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCt                                              <  2:2959034/100‑1 (MQ=255)
                                                        gagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTg                                             <  2:1113688/100‑1 (MQ=255)
                                                        gagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTg                                             <  1:3014555/100‑1 (MQ=255)
                                                          gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                           >  1:3490651/1‑100 (MQ=255)
                                                          gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                           >  1:1243830/1‑100 (MQ=255)
                                                             cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACCAGCTAAAGGAGAATTTCTCATGTCTGATATg                                        <  1:4102585/100‑1 (MQ=255)
                                                             cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGACCAACAAGCTAAAGGAGACTTTCTCATGTCTGCTATg                                        >  1:4137623/1‑100 (MQ=255)
                                                               aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                      <  1:1579893/100‑1 (MQ=255)
                                                               aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                      >  2:1350272/1‑100 (MQ=255)
                                                                 aaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                                    <  1:2944729/100‑1 (MQ=255)
                                                                  aTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGc                                   <  2:3697123/100‑1 (MQ=255)
                                                                  aTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGc                                   <  2:1829042/100‑1 (MQ=255)
                                                                    aGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAg                                 <  1:1400481/100‑1 (MQ=255)
                                                                     ggATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGt                                >  2:2693160/1‑100 (MQ=255)
                                                                      gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                               <  2:755703/100‑1 (MQ=255)
                                                                      gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                               <  2:957005/100‑1 (MQ=255)
                                                                      gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                               >  1:2152669/1‑100 (MQ=255)
                                                                       aTCAACACAAATAAAAGCCCCCTGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACCAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGt                                <  1:2511650/98‑1 (MQ=255)
                                                                       aTCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAg                              <  2:3376366/100‑1 (MQ=255)
                                                                         cAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAg                              <  1:448555/98‑1 (MQ=255)
                                                                            cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                         >  1:49950/1‑100 (MQ=255)
                                                                              caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                       >  1:3862471/1‑100 (MQ=255)
                                                                              caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                       >  2:2154617/1‑100 (MQ=255)
                                                                                aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                     <  2:1793440/100‑1 (MQ=255)
                                                                                aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                     >  1:4192684/1‑100 (MQ=255)
                                                                                   aaaaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaa                  >  2:3291488/1‑100 (MQ=255)
                                                                                     ccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGATAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaa                 >  2:3285716/3‑99 (MQ=37)
                                                                                      cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAGTTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca               >  1:2816525/2‑100 (MQ=37)
                                                                                      cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca               >  2:1882011/2‑100 (MQ=38)
                                                                                        cccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat             <  1:2417083/100‑1 (MQ=255)
                                                                                         ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg            >  2:3446544/1‑100 (MQ=38)
                                                                                         ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg            >  2:2439237/1‑100 (MQ=38)
                                                                                         ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATCTCTGACATGCTAGCAGTAGAACAAGAAACAAACAatg            >  1:2960068/1‑100 (MQ=38)
                                                                                         acccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGAAGAATTTCTAATGTATGATATGCTAACAGTAGAACAAGAAacaaaaaa              >  1:2659441/2‑98 (MQ=25)
                                                                                          cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGCAGAACAAGAAACAAACAat             >  2:3770125/1‑98 (MQ=255)
                                                                                           ccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgat          >  2:2863400/1‑100 (MQ=255)
                                                                                                aTGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTTTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTc      >  1:1868362/1‑99 (MQ=255)
                                                                                                aTGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCg     <  2:2714223/100‑1 (MQ=255)
                                                                                                aTGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCg     <  1:1706267/100‑1 (MQ=255)
                                                                                                   ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCa  >  2:3556824/1‑100 (MQ=255)
                                                                                                   ctgcAACATGCCGGGGGCGGAATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAAGGATGTTCGTCa  >  1:98606/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: