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.1% 457.6 / 3.0 137pseudogene (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/3);  new base G (62/70);  total (63/74)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

Only 100 of 137 total aligned reads displayed.

TGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTTCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGC  >  NC_000913/576041‑576231
                                                                                                |                                                                                              
tGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTGTGCCGGAAGCAAATACTGCAg                                                                                             <  1:3480118/100‑1 (MQ=255)
 gggTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTTCAgg                                                                                            >  2:1831288/1‑100 (MQ=255)
  ggTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGga                                                                                           <  1:713197/100‑1 (MQ=255)
  ggTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGga                                                                                           <  2:1238178/100‑1 (MQ=255)
    tttCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGgaag                                                                                         <  2:3365637/100‑1 (MQ=255)
    tttCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGgaag                                                                                         <  1:1942988/100‑1 (MQ=255)
         gAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACtt                                                                                    >  1:1054778/1‑100 (MQ=255)
         gAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACtt                                                                                    >  2:572599/1‑100 (MQ=255)
         gAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACtt                                                                                    >  1:1956316/1‑100 (MQ=255)
             aGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTccc                                                                                <  1:551093/100‑1 (MQ=255)
              ggTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCt                                                                               >  1:1110218/1‑100 (MQ=255)
              ggTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCt                                                                               >  2:207746/1‑100 (MQ=255)
               gTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTg                                                                              >  1:126710/1‑100 (MQ=255)
               gTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTg                                                                              >  1:1802499/1‑100 (MQ=255)
                 aGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCa                                                                            <  1:366512/100‑1 (MQ=255)
                  gTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCAt                                                                           <  1:2380450/100‑1 (MQ=255)
                    ggCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTa                                                                         <  1:2295053/100‑1 (MQ=255)
                      ccAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACTGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAAc                                                                       >  2:3191956/1‑100 (MQ=255)
                      ccAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAAc                                                                       <  2:1043205/100‑1 (MQ=255)
                       cAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACt                                                                      <  1:1355346/100‑1 (MQ=255)
                       cAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACt                                                                      <  1:1399306/100‑1 (MQ=255)
                       cAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGACCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACt                                                                      <  2:1181352/100‑1 (MQ=255)
                         ggTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTg                                                                    >  1:3612424/1‑100 (MQ=255)
                           tAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAg                                                                  <  2:2892499/100‑1 (MQ=255)
                           tAAATGTTGTTAGCGTCATATTGCAGACCCGGGGCCCAAACCTCTGAGTTTTTACCGGAACCAAATACTGCAGGAAGAACAGTCCCTGCATTAACTTGAg                                                                  <  2:1619391/100‑1 (MQ=255)
                            aaaTGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCGTTAACTTGAgt                                                                 >  2:3189764/1‑100 (MQ=255)
                             aaTGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTTCAGGAAGAACTTTCCCTGCTTTAACTTGAgtg                                                                <  1:3514123/100‑1 (MQ=255)
                             aaTGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTTCAGGAAGAACTTTCCCTGCTTTAACTTGAgtg                                                                <  1:3514775/100‑1 (MQ=255)
                             aaTGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAgtg                                                                <  2:1176639/100‑1 (MQ=255)
                             aaTGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAgtg                                                                >  1:515760/1‑100 (MQ=255)
                               tgttgtTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGACGAACTTTCCCTGCATTAACTTGAGTGTc                                                              >  2:3030086/1‑100 (MQ=255)
                                gttgttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCg                                                             <  1:797147/100‑1 (MQ=255)
                                 ttgttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCgg                                                            <  2:3797616/100‑1 (MQ=255)
                                 ttgttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCgg                                                            <  1:2600974/100‑1 (MQ=255)
                                   gttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTa                                                          >  1:2572868/1‑100 (MQ=255)
                                    ttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTAc                                                         <  2:3185622/100‑1 (MQ=255)
                                     tAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACg                                                        >  1:3172174/1‑100 (MQ=255)
                                     tAGCGTCATATTGCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACg                                                        <  2:3473497/100‑1 (MQ=255)
                                       gCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGAt                                                      <  1:1405275/100‑1 (MQ=255)
                                       gCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGAt                                                      <  2:3548388/100‑1 (MQ=255)
                                         gTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCa                                                    >  1:1403360/1‑100 (MQ=255)
                                           cATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGa                                                  >  2:823835/1‑100 (MQ=255)
                                           cATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGa                                                  >  1:3996779/1‑100 (MQ=255)
                                           cATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGGCGGTACGCTCAGa                                                  >  1:3671812/1‑100 (MQ=255)
                                            atatTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGAt                                                 <  2:1349560/100‑1 (MQ=255)
                                            atatTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGAt                                                 <  1:2171928/100‑1 (MQ=255)
                                            atatTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGAt                                                 <  2:2371227/100‑1 (MQ=255)
                                            atatTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGAt                                                 <  2:1329812/100‑1 (MQ=255)
                                            atatTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTAAGATCAGAt                                                 <  2:3837009/100‑1 (MQ=255)
                                              atTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGGTCAGAttt                                               >  1:3890499/1‑100 (MQ=255)
                                              atTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGAttt                                               >  1:3090021/1‑100 (MQ=255)
                                              atTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGAttt                                               <  1:94741/100‑1 (MQ=255)
                                              atTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGAttt                                               <  1:96525/100‑1 (MQ=255)
                                               tttCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTc                                              <  1:3192580/100‑1 (MQ=255)
                                                ttCTGACCTGCGGCCCAAGGTGCGGGGCGTGTACCGGAAGCAAATAGTGCAGGAAGATCTGTCCCTGCATTAACTTGGGTGTCGGTACGATCAGATTGCg                                             <  2:1797334/100‑1 (MQ=255)
                                                ttCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCg                                             <  1:1571693/100‑1 (MQ=255)
                                                ttCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTCCAGGAAGAACCTTCCCTGCATTAACTTGCGTGTCGGTACGATCAGATTTCg                                             <  2:1622480/100‑1 (MQ=255)
                                                 tCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTAGGTACGATCAGATTTCGc                                            >  2:1597476/1‑100 (MQ=255)
                                                  cAGACCTGCGGCCCAAACTTCTGCAGTTTTACCGGAAGCAAATACTGCAGTAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCa                                           <  1:599588/100‑1 (MQ=255)
                                                    gACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATa                                         >  1:3838025/1‑100 (MQ=255)
                                                     aCCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATCTCGCATaa                                        >  2:168520/1‑100 (MQ=255)
                                                      ccTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAg                                       >  2:1085727/1‑100 (MQ=255)
                                                      ccTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAg                                       >  1:3353299/1‑100 (MQ=255)
                                                         gCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTg                                    <  1:3087364/100‑1 (MQ=255)
                                                         gCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTg                                    <  1:1683421/100‑1 (MQ=255)
                                                          cGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGc                                   <  2:1084750/100‑1 (MQ=255)
                                                           ggCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCa                                  <  1:3805796/100‑1 (MQ=255)
                                                            gCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCAc                                 <  1:2605869/100‑1 (MQ=255)
                                                            gCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCAc                                 <  1:1334830/100‑1 (MQ=255)
                                                             cccAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCAcc                                <  1:3448305/100‑1 (MQ=255)
                                                             cccAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCAcc                                <  1:3055227/100‑1 (MQ=255)
                                                               cAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGa                              >  1:479125/1‑100 (MQ=255)
                                                                   tgtCTGCAGTTTTACCGGAAGCAAATACTGCAGTAAGAACTTCCCCTGCATTAACTTGAGTGTCGGTACGATCACAGTTCGCATAACTTGCACCGATAcc                          <  2:813846/98‑1 (MQ=255)
                                                                   cTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATAcc                          >  1:1750370/1‑100 (MQ=255)
                                                                    ttCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCg                         <  2:2507176/100‑1 (MQ=255)
                                                                      cTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGaa                       >  1:766044/1‑100 (MQ=255)
                                                                      cTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGaa                       <  2:3908036/100‑1 (MQ=255)
                                                                       tGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAAt                      >  1:3830928/1‑100 (MQ=255)
                                                                       tGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAAt                      >  2:2908303/1‑100 (MQ=255)
                                                                       tGCATTTTTAACGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAAt                      >  2:3680692/1‑100 (MQ=255)
                                                                         cATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATcc                    >  2:2579761/1‑100 (MQ=255)
                                                                         cATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATcc                    >  1:1981967/1‑100 (MQ=255)
                                                                         cATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATcc                    <  2:619976/100‑1 (MQ=255)
                                                                          aTTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCt                   >  1:2350925/1‑100 (MQ=255)
                                                                            ttttACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTc                 >  1:3190093/1‑100 (MQ=255)
                                                                             tttACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCg                <  2:2073180/100‑1 (MQ=255)
                                                                             tttACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCg                <  2:2287186/100‑1 (MQ=255)
                                                                             tttACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCg                <  2:2221404/100‑1 (MQ=255)
                                                                             tttACCAGAAGCAAATACTGCAGGAAGAACTGTACCTGCATTAACTTGAGTGTCGGTACGATCTGATTTCGCATAGGTTGCGCCGATACCGAATCCTTCg                <  2:2426302/100‑1 (MQ=255)
                                                                               tACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTa              >  1:770418/1‑100 (MQ=255)
                                                                                aCCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTAt             >  1:604021/1‑100 (MQ=255)
                                                                                 ccGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTAtt            >  1:1786028/1‑100 (MQ=255)
                                                                                 ccGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTAtt            >  2:3879825/1‑100 (MQ=255)
                                                                                   ggAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTAGTCa          >  1:3013376/1‑100 (MQ=255)
                                                                                    gAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTGGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCAt         >  2:3081550/1‑100 (MQ=255)
                                                                                      aGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAg       >  2:2920019/1‑100 (MQ=255)
                                                                                      aGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAg       <  1:1147875/100‑1 (MQ=255)
                                                                                       gCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAgg      <  1:2203838/100‑1 (MQ=255)
                                                                                       gCAAATACTGCAGGAAAAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAgt      >  2:3872247/1‑99 (MQ=255)
                                                                                        cAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGt     <  1:2721439/100‑1 (MQ=255)
                                                                                                |                                                                                              
TGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTTCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGC  >  NC_000913/576041‑576231

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: