Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 576,137 T→G pseudogene (689/1040 nt) nmpC ← DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913576,1370TG97.5% 692.8 / ‑0.8 196pseudogene (689/1040 nt)nmpCDLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base G (94/97);  total (95/101)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

CTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTTCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAG  >  NC_000913/576040‑576235
                                                                                                 |                                                                                                  
ntGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCa                                                                                                  >  1:559649/2‑100 (MQ=255)
cttggTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCa                                                                                                  <  1:104463/97‑1 (MQ=255)
cTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCa                                                                                                  <  2:4618065/100‑1 (MQ=255)
cTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCa                                                                                                  >  2:4481818/1‑100 (MQ=255)
 tGGGTTTCAGAATAGTTATTGTCCATGTAAATGTTGTTAGCGTCATATTTCAGACCTTCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAg                                                                                                 <  2:539194/100‑1 (MQ=255)
 tGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAg                                                                                                 >  1:3868138/1‑100 (MQ=255)
 tGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAg                                                                                                 <  1:2890376/100‑1 (MQ=255)
   ggTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAggg                                                                                               >  2:3682212/1‑99 (MQ=255)
   ggTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGga                                                                                               >  1:4212658/1‑100 (MQ=255)
   ggTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGga                                                                                               <  2:4596014/100‑1 (MQ=255)
      ttCAGAATAGGTGGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGTCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGgaaga                                                                                            <  2:4103172/100‑1 (MQ=255)
      ttCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGgaaga                                                                                            >  1:4097788/1‑100 (MQ=255)
      ttCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGgaaga                                                                                            <  1:2618793/100‑1 (MQ=255)
         aGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACt                                                                                         >  2:3496941/1‑100 (MQ=255)
             tAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTcc                                                                                     <  1:3577587/100‑1 (MQ=255)
             tAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTcc                                                                                     <  1:3725215/100‑1 (MQ=255)
             tAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTcc                                                                                     <  1:1220057/100‑1 (MQ=255)
              aGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTccc                                                                                    <  2:2634236/100‑1 (MQ=255)
              aGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTccc                                                                                    <  1:3103544/100‑1 (MQ=255)
               ggTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCt                                                                                   <  1:3767015/100‑1 (MQ=255)
                gTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTg                                                                                  >  1:4873518/1‑100 (MQ=255)
                gTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTg                                                                                  <  1:772352/100‑1 (MQ=255)
                gTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTg                                                                                  >  2:513584/1‑100 (MQ=255)
                 tAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGc                                                                                 <  1:4393787/100‑1 (MQ=255)
                 tAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGc                                                                                 <  2:1191984/100‑1 (MQ=255)
                  aGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCa                                                                                <  2:1819269/100‑1 (MQ=255)
                   gTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCAt                                                                               >  1:3998596/1‑100 (MQ=255)
                     ggCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCTTTa                                                                             >  2:634956/1‑100 (MQ=255)
                     ggCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTa                                                                             >  1:2578254/1‑100 (MQ=255)
                     ggCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCAGTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTa                                                                             >  1:5106698/1‑100 (MQ=255)
                      gCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTaa                                                                            >  2:3283357/1‑100 (MQ=255)
                       ccAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAAc                                                                           >  2:3168757/1‑100 (MQ=255)
                        cAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACt                                                                          <  2:4426556/100‑1 (MQ=255)
                         aGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGGAAATACTTCAGGAAGAACTTTCCCTGCATTAACtt                                                                         <  1:5386151/100‑1 (MQ=255)
                         aGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACtt                                                                         <  2:3222898/100‑1 (MQ=255)
                          ggTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACttn                                                                        <  2:2236959/100‑2 (MQ=255)
                          ggTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTg                                                                        <  2:2547722/100‑1 (MQ=255)
                          ggTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTg                                                                        >  2:3961487/1‑100 (MQ=255)
                          ggTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTg                                                                        >  1:4804981/1‑100 (MQ=255)
                          ggTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTg                                                                        >  1:4338398/1‑100 (MQ=255)
                          ggTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTg                                                                        >  2:4321540/1‑100 (MQ=255)
                          ggTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTg                                                                        <  1:3429622/100‑1 (MQ=255)
                          ggTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTg                                                                        >  2:4827689/1‑100 (MQ=255)
                            tAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAg                                                                      <  2:3099252/100‑1 (MQ=255)
                             aaaTGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAgt                                                                     <  2:692074/100‑1 (MQ=255)
                              aaTGTTGTTAGCGTCATATCTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTTCAGGAAGAACTTTCCCTGCATTAACTTGAgtg                                                                    <  2:4914490/100‑1 (MQ=255)
                                tgttgtTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTc                                                                  <  2:1996356/100‑1 (MQ=255)
                                tgttgtTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTc                                                                  >  2:1324284/1‑100 (MQ=255)
                                  ttgttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCgg                                                                <  2:4162442/100‑1 (MQ=255)
                                  ttgttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCgg                                                                <  2:1914583/100‑1 (MQ=255)
                                  ttgttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCgg                                                                <  1:1409117/100‑1 (MQ=255)
                                    gttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAATTTTCCCTGCATTAACTTGAGTGTCGGTa                                                              >  2:1527685/1‑100 (MQ=255)
                                      tAGCGTCATATTTCAGACCTGCGGCCCACACTTCTGCATTTGTGCCGGATGCAAATACTGCAGGAAGCACTTTCTCTGCATTAGCTTGAGTCTCGGTACg                                                            <  2:837995/100‑1 (MQ=255)
                                      tAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACg                                                            <  2:3623040/100‑1 (MQ=255)
                                          gTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCa                                                        >  1:1964270/1‑100 (MQ=255)
                                          gTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCa                                                        >  1:4146821/1‑100 (MQ=255)
                                          gTCATATTTCAGACCTGCGGCACAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCa                                                        >  2:278166/1‑100 (MQ=255)
                                           tCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAg                                                       <  1:1135209/100‑1 (MQ=255)
                                            cATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGa                                                      <  1:2122629/100‑1 (MQ=255)
                                              tatTTCAGACCTGCGGCCCAAAGTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGAtt                                                    >  1:4652849/1‑100 (MQ=255)
                                               atTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGAttt                                                   >  1:3754067/1‑100 (MQ=255)
                                               atTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGAttt                                                   >  1:5145123/1‑100 (MQ=255)
                                                  tCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGc                                                <  2:553078/100‑1 (MQ=255)
                                                  tCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGc                                                >  2:3043869/1‑100 (MQ=255)
                                                  tCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGc                                                <  1:3198329/100‑1 (MQ=255)
                                                    aGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCan                                              <  1:224758/100‑2 (MQ=255)
                                                    aGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCAt                                              >  2:4927983/1‑100 (MQ=255)
                                                    aGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCAt                                              >  2:159184/1‑100 (MQ=255)
                                                    aGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCAt                                              <  1:4158813/100‑1 (MQ=255)
                                                     gACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATa                                             >  1:1166755/1‑100 (MQ=255)
                                                     gACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATa                                             <  1:2632716/100‑1 (MQ=255)
                                                       ccTGCGGCCCAAACTTCTGTATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAATTTGATTGTCGGTACGATCAGATTTCGCACAAg                                           >  2:5035373/1‑100 (MQ=255)
                                                        cTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGt                                          >  2:338683/1‑100 (MQ=255)
                                                         tGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGtt                                         >  1:1580281/1‑100 (MQ=255)
                                                            ggCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCa                                      <  1:1833994/100‑1 (MQ=255)
                                                              cccAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTCGTGTGTCGGTACGATCAGATTTCGCATAAGTTGCAcc                                    <  2:5370901/100‑1 (MQ=255)
                                                              cccAAACTTCTGCATTATTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCAcc                                    >  2:5318236/1‑100 (MQ=255)
                                                               ccAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGTACCg                                   >  1:2654052/1‑100 (MQ=255)
                                                                 aaaCTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGAt                                 <  1:4896584/100‑1 (MQ=255)
                                                                    cTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATAcc                              <  2:2067601/100‑1 (MQ=255)
                                                                       cTGCATTTTTACCGGAAGCAAATCCTGCAGGAAGAACTTTCACTGCATTAACTTGAGTGTAGGTACGATCAGATTTCGCATAAGTTGCACCGATACTGaa                           >  2:2087145/1‑100 (MQ=255)
                                                                        tGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAAt                          <  2:3832741/100‑1 (MQ=255)
                                                                          cATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATcc                        >  2:1687105/1‑100 (MQ=255)
                                                                          cATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATcc                        <  1:4257241/100‑1 (MQ=255)
                                                                             ttttACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTc                     >  1:5007425/1‑100 (MQ=255)
                                                                              tttACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCg                    <  1:3970172/100‑1 (MQ=255)
                                                                                tACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTa                  <  1:3952010/100‑1 (MQ=255)
                                                                                  ccGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTAtt                >  2:1743287/1‑100 (MQ=255)
                                                                                  ccGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTAtt                >  1:5035466/1‑100 (MQ=255)
                                                                                  ccGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTAtt                >  1:950156/1‑100 (MQ=255)
                                                                                   cGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTc               <  2:119763/100‑1 (MQ=255)
                                                                                   cGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTc               <  2:3446636/100‑1 (MQ=255)
                                                                                    ggAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCa              >  2:2615766/1‑100 (MQ=255)
                                                                                    ggAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCa              >  1:617360/1‑100 (MQ=255)
                                                                                      aaGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATa            <  2:2388860/100‑1 (MQ=255)
                                                                                         cAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGt         <  1:470251/100‑1 (MQ=255)
                                                                                         cAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGt         >  2:654740/1‑100 (MQ=255)
                                                                                           aaTACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAg       <  1:3722054/100‑1 (MQ=255)
                                                                                             tACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCa     >  1:991225/1‑100 (MQ=255)
                                                                                                ttCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCagag  >  1:1359172/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
CTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTTCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAG  >  NC_000913/576040‑576235

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: