Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,804 G→T 16.0% intergenic (+431/+63) icd → / ← ymfD isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,8040GT16.0% 575.1 / 87.1 213intergenic (+431/+63)icd/ymfDisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (85/93);  new base T (19/15);  total (105/108)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 219 total aligned reads displayed.

ATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAA  >  NC_000913/1196705‑1196903
                                                                                                   |                                                                                                   
  cAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAgcg                                                                                                   >  1:3287088/1‑100 (MQ=255)
     ctactCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCgg                                                                                                <  2:260392/100‑1 (MQ=255)
        ctCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGccc                                                                                             <  2:2673313/100‑1 (MQ=255)
        ctCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGccc                                                                                             >  2:1104887/1‑100 (MQ=255)
         tCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCACTCACCTTAGCGTTCAGGTACATACAAAATGTAAAAAAGGCCGCGCTCGCGGcctt                                                                                            <  2:630241/100‑3 (MQ=255)
          cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGccttt                                                                                           <  2:1664793/100‑4 (MQ=255)
          cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGccttt                                                                                           <  1:2233140/100‑4 (MQ=255)
          cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCt                                                                                           >  2:847806/1‑100 (MQ=255)
            cccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTc                                                                                         >  2:1425040/1‑100 (MQ=255)
            cccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTc                                                                                         >  1:1993964/1‑100 (MQ=255)
             ccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcctttttt                                                                                        <  1:1076688/100‑7 (MQ=255)
             ccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcctttttt                                                                                        <  2:2844346/100‑7 (MQ=255)
             ccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTca                                                                                        <  2:3558587/100‑1 (MQ=255)
              cGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcctttttta                                                                                       <  2:2311982/100‑8 (MQ=255)
               gggACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTcaca                                                                                      <  1:1322616/100‑1 (MQ=255)
               gggAATCGTCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACAGGCGGAATGTAAAAAAGGCCGCGCTCGCGGccttttttac                                                                                      <  2:1087764/100‑9 (MQ=255)
                   cTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCacataca                                                                                  >  1:4335517/1‑100 (MQ=255)
                   cTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCacataca                                                                                  <  2:1875211/100‑1 (MQ=255)
                     cGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATc                                                                                >  2:3381439/1‑100 (MQ=255)
                     cGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATc                                                                                >  1:1944637/1‑100 (MQ=255)
                      gCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATc                                                                                <  1:3598779/99‑1 (MQ=255)
                      gCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGAAAAAAAAAAAAAAGGCCGCGAGCGCGGCCCCTTCacataca                                                                                  <  1:776422/97‑1 (MQ=255)
                        cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                             <  1:2842626/100‑1 (MQ=255)
                        cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                             >  1:3921978/1‑100 (MQ=255)
                        cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                             >  2:799585/1‑100 (MQ=255)
                        cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                             <  2:1710125/100‑1 (MQ=255)
                        cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                             <  2:1993972/100‑1 (MQ=255)
                        cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                             <  2:2725090/100‑1 (MQ=255)
                          ggggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAg                                                                           <  2:3030223/100‑1 (MQ=255)
                          gacgACAGCCAACAGGCATTGGGTGCAATCGCCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAg                                                                           <  1:2731300/97‑1 (MQ=255)
                           gggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGt                                                                          >  2:3193611/1‑100 (MQ=255)
                           gggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGt                                                                          <  2:413044/100‑1 (MQ=255)
                           gggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGt                                                                          >  1:92497/1‑100 (MQ=255)
                           gcgACAGCCACTAGTCATTGGGTGCTATCACCTTTGCGTTCAGGTACAGCCCTTATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAg                                                                           <  1:882783/97‑1 (MQ=255)
                               cAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTg                                                                      <  1:3518882/100‑1 (MQ=255)
                               cAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTg                                                                      <  1:60062/100‑1 (MQ=255)
                               cAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTg                                                                      <  1:4224982/100‑1 (MQ=255)
                                aGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTga                                                                     >  2:3243849/1‑100 (MQ=255)
                                aGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTga                                                                     >  2:1127930/1‑100 (MQ=255)
                                 gCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAACAAAGGCCGCGAGCGCGGCCCCTTCACATAAAACTTTAGTACt                                                                       >  1:1278192/1‑97 (MQ=255)
                                 gCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTgag                                                                    >  1:2254013/1‑100 (MQ=255)
                                 gCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTgag                                                                    <  2:1578588/100‑1 (MQ=255)
                                  ccAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTgaga                                                                   >  2:2505433/1‑100 (MQ=255)
                                    aaCAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGAc                                                                  >  1:411819/1‑99 (MQ=255)
                                         gCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTa                                                            >  1:3750503/1‑100 (MQ=255)
                                          cATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTaa                                                           >  2:3890476/1‑100 (MQ=255)
                                           aTTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAAc                                                          >  1:3379582/1‑100 (MQ=255)
                                              gggTGCAATCACCTTAGCGTTCAGGTACATGCGGCATGTAAAAAAGCCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAAcc                                                         >  1:1108556/1‑98 (MQ=255)
                                              gggTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTa                                                       >  2:4053544/1‑100 (MQ=255)
                                               ggTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAg                                                      >  2:3783781/1‑100 (MQ=255)
                                               ggTGCAAGCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAg                                                      <  1:3553393/100‑1 (MQ=255)
                                                gTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAgg                                                     >  2:1158409/1‑100 (MQ=255)
                                                 tGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACAGCTTTAGTACTGAGACTGTTTAACCTAggg                                                    <  1:1794954/100‑1 (MQ=255)
                                                     aTCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAatt                                                >  2:4203744/1‑100 (MQ=255)
                                                      tCACCTTAGCGTTCAGGTACATCCACAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAatta                                               <  1:3965982/100‑1 (MQ=255)
                                                         ccTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAAttatta                                            <  2:719043/100‑1 (MQ=255)
                                                         ccTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAAttatta                                            >  1:1374903/1‑100 (MQ=255)
                                                         ccTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTAATGAGACTGTTTAACATAGGGAAttatta                                            >  2:3960412/1‑100 (MQ=255)
                                                          cTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATtattat                                           >  1:2415862/1‑100 (MQ=255)
                                                            tAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                         >  2:2096820/1‑100 (MQ=255)
                                                            tAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                         <  1:460728/100‑1 (MQ=255)
                                                            tAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                         <  2:520571/100‑1 (MQ=255)
                                                             aGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGt                                        >  2:3517883/1‑100 (MQ=255)
                                                               cGTTCAGGTACATGCCGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGtat                                      <  2:4340133/100‑1 (MQ=255)
                                                                    aGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTatat                                 >  2:2948505/1‑100 (MQ=255)
                                                                     ggTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATAtt                                >  1:3144660/1‑100 (MQ=255)
                                                                       tACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGc                              <  1:2077721/100‑1 (MQ=255)
                                                                         cATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCAt                            >  2:4099990/1‑100 (MQ=255)
                                                                         cATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCAt                            <  1:3300073/100‑1 (MQ=255)
                                                                         cATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCAt                            >  2:1059323/1‑100 (MQ=255)
                                                                            gCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTAc                         >  1:1312205/1‑100 (MQ=255)
                                                                            gCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTAc                         <  2:777095/100‑1 (MQ=255)
                                                                              ggAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGACTTATTATCGTATTATATTGCATTACaa                       >  1:3136594/1‑100 (MQ=255)
                                                                              ggAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACaa                       <  2:3705569/100‑1 (MQ=255)
                                                                                aaTGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACt                     >  1:39432/1‑100 (MQ=255)
                                                                                  tGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTcc                   <  1:3764167/100‑1 (MQ=255)
                                                                                   gTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTccc                  <  2:2751718/100‑1 (MQ=255)
                                                                                   gTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTcc                   >  1:834937/1‑99 (MQ=255)
                                                                                    tAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCa                 >  2:3046753/1‑100 (MQ=255)
                                                                                    tAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCa                 >  1:2496182/1‑100 (MQ=255)
                                                                                      gaaggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAAc               >  2:2268552/6‑100 (MQ=255)
                                                                                      aaaaaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAAc               <  2:2374299/100‑1 (MQ=255)
                                                                                      aaaaaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAAc               <  1:1578275/100‑1 (MQ=255)
                                                                                        aggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAg             >  2:3125267/4‑100 (MQ=255)
                                                                                        aggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAg             >  1:21142/4‑100 (MQ=255)
                                                                                        aggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAg             >  2:3023195/4‑100 (MQ=255)
                                                                                         ngggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGCCTGTTTAACCTAGGGAATTCTTATCGTATTATATTGCATTACGACTCCCAACAg             >  1:2592869/3‑99 (MQ=255)
                                                                                         ngggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGt            >  2:1222567/3‑100 (MQ=255)
                                                                                         ggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGt            >  2:3316143/3‑100 (MQ=255)
                                                                                          aGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTg           >  1:400247/1‑100 (MQ=255)
                                                                                          aGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTg           <  1:4348112/100‑1 (MQ=255)
                                                                                           ggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGa          >  2:437735/1‑100 (MQ=255)
                                                                                           ggCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGa          <  1:664681/100‑1 (MQ=255)
                                                                                           ggCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGa          <  2:3287104/100‑1 (MQ=255)
                                                                                             ccGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACa        <  2:4022823/100‑1 (MQ=255)
                                                                                               gcgAGCGCGGCTCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACa        <  1:3338252/98‑1 (MQ=255)
                                                                                                cgAGCGCGGCTCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACatat     <  1:3219747/100‑1 (MQ=255)
                                                                                                cgAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACatat     >  2:1912320/1‑100 (MQ=255)
                                                                                                 gAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATCGCATTACAACTCCCAACAGTGACATATg    <  2:3223318/100‑1 (MQ=255)
                                                                                                  aGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGa   <  1:1887095/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAA  >  NC_000913/1196705‑1196903

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: