Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,765,014 T→C D47G (GAC→GGC)  yfjL ← CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,765,0140TC98.2% 398.9 / ‑2.9 114D47G (GAC→GGC) yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (60/52);  total (60/54)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 114 total aligned reads displayed.

AGCACCTTTTGAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGTCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTCCTGAAATAAGA  >  NC_000913/2764917‑2765110
                                                                                                 |                                                                                                
aGCACCTTTTGAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCg                                                                                                >  2:1108732/1‑100 (MQ=255)
 gCACCTTTTGAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGa                                                                                               >  2:705188/1‑100 (MQ=255)
  cACCTTTTGAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGat                                                                                              <  1:1059234/100‑1 (MQ=255)
  cACCTTTTGAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGat                                                                                              <  2:2165192/100‑1 (MQ=255)
   aCCTTTTGAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGata                                                                                             <  1:285469/100‑1 (MQ=255)
    ccTTTTGAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGata                                                                                             <  1:631358/99‑1 (MQ=255)
     cTTTTGAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATc                                                                                           >  1:1110875/1‑100 (MQ=255)
     cTTTTGAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATc                                                                                           >  2:325398/1‑100 (MQ=255)
     cTTTTGAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATc                                                                                           <  1:1387218/100‑1 (MQ=255)
          gAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCtg                                                                                      <  1:2815435/100‑1 (MQ=255)
          gAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCGGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCtg                                                                                      <  2:1001800/100‑1 (MQ=255)
           aaTTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCtgt                                                                                     >  2:2518013/1‑100 (MQ=255)
            aTTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCtgtg                                                                                    <  1:2006841/100‑1 (MQ=255)
            aTTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCtgtg                                                                                    <  1:565766/100‑1 (MQ=255)
             ttCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGTCGATATCGTCTGTGa                                                                                   <  1:1415023/100‑1 (MQ=255)
             ttCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCtgtg                                                                                    >  1:2609308/1‑99 (MQ=255)
             ttCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGa                                                                                   <  2:1015088/100‑1 (MQ=255)
              tCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGAc                                                                                  <  1:287340/100‑1 (MQ=255)
              tCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGAc                                                                                  >  1:194319/1‑100 (MQ=255)
              tCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGAc                                                                                  <  2:2076754/100‑1 (MQ=255)
               cTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGAcc                                                                                 <  2:847768/100‑1 (MQ=255)
                tGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCt                                                                                <  1:1128223/100‑1 (MQ=255)
                 gAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGATCGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTc                                                                               >  2:2892543/1‑100 (MQ=255)
                 gAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTc                                                                               >  1:618224/1‑100 (MQ=255)
                       tGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAgcg                                                                         <  2:2512427/100‑1 (MQ=255)
                         tCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGca                                                                       >  1:528252/1‑100 (MQ=255)
                          ccGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGcac                                                                      >  1:1309602/1‑100 (MQ=255)
                            gAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCacac                                                                    <  1:1377394/100‑1 (MQ=255)
                            gAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCacac                                                                    <  2:802711/100‑1 (MQ=255)
                               aTGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGc                                                                 >  1:2742723/1‑100 (MQ=255)
                                tGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCg                                                                <  2:2077909/100‑1 (MQ=255)
                                 gTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGcag                                                               >  1:1283994/1‑98 (MQ=255)
                                 gTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCgg                                                               >  1:1102864/1‑100 (MQ=255)
                                 gTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCgg                                                               >  2:2298973/1‑100 (MQ=255)
                                   tCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGAt                                                             >  1:830751/1‑100 (MQ=255)
                                     cAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAg                                                           >  1:632449/1‑100 (MQ=255)
                                        cGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCt                                                        <  2:805243/100‑1 (MQ=255)
                                        cGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCt                                                        <  2:81215/100‑1 (MQ=255)
                                         gtgtTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCtt                                                       >  1:934977/1‑100 (MQ=255)
                                          tgtTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTg                                                      >  2:1616749/1‑100 (MQ=255)
                                          tgtTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTg                                                      >  1:1635070/1‑100 (MQ=255)
                                           gtTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGt                                                     >  1:2203931/1‑100 (MQ=255)
                                           gtTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGt                                                     >  2:339862/1‑100 (MQ=255)
                                             tttGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCa                                                   <  1:2517981/100‑1 (MQ=255)
                                               tGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACg                                                 >  2:2878755/1‑100 (MQ=255)
                                               tGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACg                                                 >  2:2878473/1‑100 (MQ=255)
                                               tGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACg                                                 <  1:581935/100‑1 (MQ=255)
                                               tGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACg                                                 <  2:238794/100‑1 (MQ=255)
                                                  cGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAg                                              <  2:1649430/100‑1 (MQ=255)
                                                   ggTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGc                                             >  2:2372576/1‑100 (MQ=255)
                                                   ggTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGc                                             >  2:2527358/1‑100 (MQ=255)
                                                   ggTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGc                                             >  2:473110/1‑100 (MQ=255)
                                                   ggTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGc                                             >  1:456816/1‑100 (MQ=255)
                                                    gTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCa                                            <  2:1994376/100‑1 (MQ=255)
                                                     tGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCAt                                           <  2:2814705/100‑1 (MQ=255)
                                                      gCTTTTTACCTGAACGAAGTCGGTGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCAtt                                          >  2:706569/1‑100 (MQ=255)
                                                        tttttACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCg                                        <  1:881734/100‑1 (MQ=255)
                                                        tttttACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCg                                        <  2:618230/100‑1 (MQ=255)
                                                        tttttACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCg                                        >  2:1823213/1‑100 (MQ=255)
                                                         ttttACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCgg                                       >  2:1073700/1‑100 (MQ=255)
                                                         ttttACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCgg                                       <  2:354522/100‑1 (MQ=255)
                                                          tttACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCggg                                      <  1:1227351/100‑1 (MQ=255)
                                                            tACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTg                                    >  2:333149/1‑100 (MQ=255)
                                                            tACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTg                                    <  1:1684311/100‑1 (MQ=255)
                                                            tACCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTg                                    <  1:1967226/100‑1 (MQ=255)
                                                             nccTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGa                                   >  1:937899/2‑100 (MQ=255)
                                                             aCCTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGa                                   <  2:718582/100‑1 (MQ=255)
                                                              ccTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGAc                                  >  1:2440861/1‑100 (MQ=255)
                                                              ccTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGAc                                  <  1:1607568/100‑1 (MQ=255)
                                                               cTGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACa                                 >  1:2821545/1‑100 (MQ=255)
                                                                tGAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACAt                                >  2:1156793/1‑100 (MQ=255)
                                                                 gAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATg                               >  2:1336931/1‑100 (MQ=255)
                                                                 gAACGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATg                               >  2:2533736/1‑100 (MQ=255)
                                                                  aaCGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATgg                              >  1:2567020/1‑100 (MQ=255)
                                                                   aCGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGt                             >  1:1533751/1‑100 (MQ=255)
                                                                    cGAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTa                            <  1:1435901/100‑1 (MQ=255)
                                                                     gAAGTCGATGTATCCATCAGATTCAATGCCGGTATCGACTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTaa                           <  1:1951938/100‑1 (MQ=255)
                                                                     gAAGTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTaa                           <  2:1974390/100‑1 (MQ=255)
                                                                        gTCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCg                        <  2:236050/100‑1 (MQ=255)
                                                                         tCGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCgg                       <  2:332437/100‑1 (MQ=255)
                                                                          cGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGc                      >  1:1252951/1‑100 (MQ=255)
                                                                          cGATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGc                      <  2:1468775/100‑1 (MQ=255)
                                                                           gATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCt                     >  1:2347734/1‑100 (MQ=255)
                                                                           gATGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCt                     <  2:2241861/100‑1 (MQ=255)
                                                                             tGTATCCATCAGATACAATGTCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTAc                   <  2:293424/100‑1 (MQ=255)
                                                                             tGTATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTAc                   <  1:705182/100‑1 (MQ=255)
                                                                               tATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGc                 <  2:2821583/100‑1 (MQ=255)
                                                                               tATCCATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGc                 >  1:2086590/1‑100 (MQ=255)
                                                                                  ccATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGt              >  2:1538814/1‑100 (MQ=255)
                                                                                   cATCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTc             >  2:1789118/1‑100 (MQ=255)
                                                                                    aTCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTcc            <  2:2705394/100‑1 (MQ=255)
                                                                                    aTCAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCAATCAAGTGACATAGTAAGCGACTACGCAGTcc            >  1:2587320/1‑100 (MQ=255)
                                                                                      cAGATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTCCTg          <  1:2527326/100‑1 (MQ=255)
                                                                                        gATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTCCTGaa        <  1:414533/100‑1 (MQ=255)
                                                                                        gATACAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCAGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTCCTGaa        <  2:1284012/100‑1 (MQ=255)
                                                                                           aCAATGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTCCTGAAATa     <  1:483932/100‑1 (MQ=255)
                                                                                             aaTGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTCCTGAAATAAg   <  1:2833223/100‑1 (MQ=255)
                                                                                             aaTGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTCCTGAAATAAg   >  1:2301063/1‑100 (MQ=255)
                                                                                              aTGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTCCTGAAATAAGa  >  2:1390072/1‑100 (MQ=255)
                                                                                              aTGCCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTCCTGAAATAAGa  >  2:1265854/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AGCACCTTTTGAATTCTGAACGATGTCCGAAATGTTCCAACGTGTTTTGCCGGTGCTTTTTACCTGAACGAAGTCGATGTATCCATCAGATACAATGTCGATATCGTCTGTGACCTCGCAGCGCACACTGCGGATAGTCTTGTCACGAAGCATTCGGGTGACATGGTAAGCGGCTACGCAGTCCTGAAATAAGA  >  NC_000913/2764917‑2765110

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