Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,803 A→C 15.4% intergenic (+430/+64) 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,8030AC15.4% 538.1 / 93.0 215intergenic (+430/+64)icd/ymfDisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (90/92);  new base C (18/15);  total (108/107)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

Only 100 of 219 total aligned reads displayed.

TATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGA  >  NC_000913/1196704‑1196902
                                                                                                   |                                                                                                   
tATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGa                                                                                                     >  1:2903134/1‑100 (MQ=255)
   cAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAgcg                                                                                                  >  1:3287088/1‑100 (MQ=255)
    aaCTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAGAAGGCCGCGAgcgc                                                                                                 <  1:2220025/100‑1 (MQ=255)
      ctactCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCgg                                                                                               <  2:260392/100‑1 (MQ=255)
         ctCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGccc                                                                                            <  2:2673313/100‑1 (MQ=255)
         ctCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGccc                                                                                            >  1:2751362/1‑100 (MQ=255)
          tCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCACTCACCTTAGCGTTCAGGTACATACAAAATGTAAAAAAGGCCGCGCTCGCGGcctt                                                                                           <  2:630241/100‑3 (MQ=255)
           cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGccttt                                                                                          <  2:1664793/100‑4 (MQ=255)
           cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGccttt                                                                                          <  1:2237810/100‑4 (MQ=255)
           cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCt                                                                                          >  2:847806/1‑100 (MQ=255)
             cccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTc                                                                                        >  2:1425040/1‑100 (MQ=255)
              ccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATTCGGAATGTAAAAAAGGCCGTGCTCGCGGcctttttt                                                                                       <  1:1396037/100‑7 (MQ=255)
              ccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcctttttt                                                                                       <  1:1076688/100‑7 (MQ=255)
              ccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcctttttt                                                                                       <  2:2844346/100‑7 (MQ=255)
              ccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTca                                                                                       <  1:2091927/100‑1 (MQ=255)
              ccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTca                                                                                       <  2:3558587/100‑1 (MQ=255)
               cGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcctttttta                                                                                      <  2:2311982/100‑8 (MQ=255)
               cGGGACTCCGCAGGGGTCAGACAACAGGCATTGGCTGCAAGTATCTTAGCCTTCAGGTACATGCGGAATGCCAAAAAGGCCGCGCTCgccgccttttt                                                                                        <  1:1883017/98‑6 (MQ=255)
                    cTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCacataca                                                                                 >  1:4335517/1‑100 (MQ=255)
                    cTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCacataca                                                                                 <  2:1875211/100‑1 (MQ=255)
                     tCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGTGGCCCCTTCacataca                                                                                 <  1:1317708/99‑1 (MQ=255)
                      cGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATc                                                                               >  2:3381439/1‑100 (MQ=255)
                       gCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATc                                                                               <  1:3598779/99‑1 (MQ=255)
                       gCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCt                                                                              >  1:1763501/1‑100 (MQ=255)
                         cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                            >  1:3921978/1‑100 (MQ=255)
                         cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                            <  2:2725090/100‑1 (MQ=255)
                         cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                            <  2:1710125/100‑1 (MQ=255)
                         cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                            <  2:1993972/100‑1 (MQ=255)
                         cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                            >  2:799585/1‑100 (MQ=255)
                           ggggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAg                                                                          <  2:3030223/100‑1 (MQ=255)
                            gggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGt                                                                         >  2:3193611/1‑100 (MQ=255)
                            gggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGt                                                                         <  2:413044/100‑1 (MQ=255)
                                cAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTg                                                                     <  2:1142461/100‑1 (MQ=255)
                                cAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTg                                                                     <  1:60062/100‑1 (MQ=255)
                                cAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTg                                                                     <  1:4224982/100‑1 (MQ=255)
                                cAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTg                                                                     <  1:3518882/100‑1 (MQ=255)
                                 aGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTga                                                                    >  2:3243849/1‑100 (MQ=255)
                                  gCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAACAAAGGCCGCGAGCGCGGCCCCTTCACATAAAACTTTAGTACt                                                                      >  1:1278192/1‑97 (MQ=255)
                                  gCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTgag                                                                   <  2:1578588/100‑1 (MQ=255)
                                  gCCAACAGGCATTGGGTGCAATCACATTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGACGCGAGCGCGGCCACTTCACATACATCTTTAGTACTgag                                                                   <  1:889593/100‑1 (MQ=255)
                                   ccAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTgaga                                                                  >  2:2505433/1‑100 (MQ=255)
                                     aaCAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGAc                                                                 >  1:411819/1‑99 (MQ=255)
                                          gCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTa                                                           <  1:2007941/100‑1 (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)
                                                ggTGCAATCACCTTAGCATTCAAGTAAATAATGAATTTAAAAAATACCACGAACGGAACCTTTTTACATACTTCTTTAATTCTCATACTTTTTAACTTa                                                      >  1:198109/1‑99 (MQ=255)
                                                ggTGCAAGCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAg                                                     <  1:3553393/100‑1 (MQ=255)
                                                 gTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAgg                                                    >  2:1158409/1‑100 (MQ=255)
                                                 gTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAgg                                                    <  1:855362/100‑1 (MQ=255)
                                                 gTGCAATCACCGTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAgg                                                    <  1:949180/100‑1 (MQ=255)
                                                 gTGAAATCAACTTAGCGTTCAGGTACAAAAAAAAAAAATAAAAGGCCGCGAGCGCGGCTCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAgn                                                    <  1:2396913/100‑2 (MQ=255)
                                                   gCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGa                                                  <  2:1066519/100‑1 (MQ=255)
                                                    cAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTGTAACCTAGGGaa                                                 >  1:254007/1‑100 (MQ=255)
                                                      aTCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAatt                                               >  2:4203744/1‑100 (MQ=255)
                                                       tCACCTTAGCGTTCAGGTACATCCACAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAatta                                              <  1:3965982/100‑1 (MQ=255)
                                                          ccTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAAttatta                                           <  2:719043/100‑1 (MQ=255)
                                                          ccTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTAATGAGACTGTTTAACATAGGGAAttatta                                           >  2:3960412/1‑100 (MQ=255)
                                                             tAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                        <  2:520571/100‑1 (MQ=255)
                                                             tAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                        >  2:2096820/1‑100 (MQ=255)
                                                             tAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                        <  1:460728/100‑1 (MQ=255)
                                                             tAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                        <  1:784/100‑1 (MQ=255)
                                                              aGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGt                                       >  2:3517883/1‑100 (MQ=255)
                                                                cGTTCAGGTACATGCCGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGtat                                     <  2:4340133/100‑1 (MQ=255)
                                                                     aGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTatat                                >  2:1101821/1‑100 (MQ=255)
                                                                     aGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTatat                                >  2:2948505/1‑100 (MQ=255)
                                                                      ggTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATAtt                               >  1:3144660/1‑100 (MQ=255)
                                                                         aCATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAATACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCa                            <  1:2436415/100‑1 (MQ=255)
                                                                          cATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCAt                           >  2:4099990/1‑100 (MQ=255)
                                                                          cATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCAt                           <  1:3300073/100‑1 (MQ=255)
                                                                            tGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTa                         >  1:1332146/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)
                                                                                   tGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACACACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTcc                  >  1:2633237/1‑100 (MQ=255)
                                                                                    gTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTccc                 <  2:2751718/100‑1 (MQ=255)
                                                                                     tAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCa                >  2:3046753/1‑100 (MQ=255)
                                                                                       gaaggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAAc              >  2:2268552/6‑100 (MQ=255)
                                                                                       ccgatctccGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAAc              <  1:1905532/93‑1 (MQ=255)
                                                                                       aaaaaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAAc              <  2:2374299/100‑1 (MQ=255)
                                                                                         aggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAg            >  2:3125267/4‑100 (MQ=255)
                                                                                         aggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAg            >  2:3023195/4‑100 (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)
                                                                                               cgcgAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACat      <  2:1112093/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)
                                                                                                   |                                                                                                   
TATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGA  >  NC_000913/1196704‑1196902

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: