Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,197,340 T→C 63.0% T65A (ACT→GCT)  ymfD ← e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,3400TC63.0% 102.6 / 116.4 119T65A (ACT→GCT) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (19/25);  new base C (36/39);  total (55/64)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

Only 100 of 120 total aligned reads displayed.

TTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACT  >  NC_000913/1197241‑1197438
                                                                                                   |                                                                                                  
ttGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGt                                                                                                    >  2:2848505/1‑100 (MQ=255)
 tGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTa                                                                                                   >  2:516156/1‑100 (MQ=255)
   cATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAAc                                                                                                 <  1:2390978/100‑1 (MQ=255)
   cATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAAc                                                                                                 >  2:2023586/1‑100 (MQ=255)
    aTTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACt                                                                                                <  2:474477/100‑1 (MQ=255)
       tttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtca                                                                                             <  1:3261852/100‑1 (MQ=255)
        ttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtcat                                                                                            <  1:36475/100‑1 (MQ=255)
         tataATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtcat                                                                                            <  1:2190014/99‑1 (MQ=255)
          ataataTCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTtcatca                                                                                          >  1:92078/1‑100 (MQ=255)
           taataTCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGTCTTTTCGAGTCTAGAAAAGCAACTTCATCaa                                                                                         >  1:2501845/1‑100 (MQ=255)
           taataTCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAATTTCATCaa                                                                                         >  1:990145/1‑100 (MQ=255)
           taataTCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGACTTTTCGAGTCTAGAAAAGCAACTTCATCaa                                                                                         >  1:2992348/1‑100 (MQ=255)
             atatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAAt                                                                                       <  1:1259705/100‑1 (MQ=255)
             atatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCaact                                                                                       >  1:1640587/1‑98 (MQ=255)
              tatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAAtt                                                                                      <  1:1368171/100‑1 (MQ=255)
              tatCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAtt                                                                                      <  1:951156/100‑1 (MQ=255)
               atCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAAttt                                                                                     >  2:3266315/1‑100 (MQ=255)
                 cGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTAc                                                                                   <  2:1432055/100‑1 (MQ=255)
                  tccgatcttAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACt                                                                                  <  1:2838272/92‑1 (MQ=255)
                  gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACt                                                                                  >  2:58159/1‑100 (MQ=255)
                  gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACt                                                                                  <  2:1256352/100‑1 (MQ=255)
                    gggACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGa                                                                                >  1:3161553/1‑100 (MQ=255)
                     ggACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGAt                                                                               >  1:3377354/1‑100 (MQ=255)
                      gACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGAtt                                                                              <  1:3051324/100‑1 (MQ=255)
                       aCATATTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGAt                                                                               <  1:2882519/98‑1 (MQ=255)
                       aCATAGTCAATAATTTTAGTTTTAATACCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTa                                                                             <  2:2750133/100‑1 (MQ=255)
                        cATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTaa                                                                            <  2:49350/100‑1 (MQ=255)
                          tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTTTAGAAAAGCAACTTCATCAAATTTACTGATTAAtt                                                                          >  2:1974538/1‑100 (MQ=255)
                          tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAAtt                                                                          >  2:2896670/1‑100 (MQ=255)
                          tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACtagatcgga                                                                          >  2:2838295/1‑92 (MQ=255)
                          tAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAAtt                                                                          >  1:1356605/1‑100 (MQ=255)
                           aGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATtc                                                                         >  2:320470/1‑100 (MQ=255)
                           aGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATtc                                                                         >  2:3373397/1‑100 (MQ=255)
                           aGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtc                                                                         >  2:2546230/1‑100 (MQ=255)
                             tCAATAATTTTAGTTTTNATTNCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATNAAATTTACTGATTAATtcat                                                                       <  1:2130764/100‑1 (MQ=255)
                             tCAATAATTTTAGTTTTAATTTCTCTAATAATTAGCGATATTTCAAGTTACCTTATCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAAt                                                                           <  1:244318/96‑1 (MQ=255)
                              cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                      >  2:1733718/1‑100 (MQ=255)
                              cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATtcatc                                                                      <  1:119802/100‑1 (MQ=255)
                               aataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATtcatca                                                                     <  1:58159/100‑1 (MQ=255)
                                ataatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAg                                                                    <  2:2101222/100‑1 (MQ=255)
                                   atTTTACTTTTGATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCTAGTCTAGAAAAGCAACTTCATCATATTTACTGATTAATtcatca                                                                     <  1:3091594/96‑1 (MQ=255)
                                    ttttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAata                                                                <  1:436712/100‑1 (MQ=255)
                                     tttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTACCTTCATCAAATTTACTGATTAATTCATAAGAatat                                                               >  1:1131805/1‑100 (MQ=255)
                                     tttAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAatat                                                               >  2:1638685/1‑100 (MQ=255)
                                         gTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTaa                                                           >  2:1960475/1‑100 (MQ=255)
                                           tttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGt                                                         >  1:793372/1‑100 (MQ=255)
                                             tAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGttt                                                       <  1:272195/100‑1 (MQ=255)
                                             tAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGttt                                                       >  2:319980/1‑100 (MQ=255)
                                             ccgatcttctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGttt                                                       <  1:104559/93‑1 (MQ=255)
                                              aaTTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGtttt                                                      >  1:2901359/1‑100 (MQ=255)
                                               aTTTCTCTAATAATTTGATCTCTTTCAAGTTGCCTTATAGAGTCTAGAAAAGTAACTTCATAAAATTTACTGATTAATTCATCAGAATATCTAAGtttgc                                                     <  1:31472/100‑3 (MQ=255)
                                                 tCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTcc                                                    >  1:55580/1‑99 (MQ=255)
                                                 tCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCg                                                   <  2:2375499/100‑1 (MQ=255)
                                                  ccTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGc                                                  <  1:694067/100‑1 (MQ=255)
                                                   ctctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGCGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCa                                                 >  2:2389363/1‑100 (MQ=255)
                                                   ctctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCa                                                 >  1:2252350/1‑100 (MQ=255)
                                                   ctctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCa                                                 <  2:136056/100‑1 (MQ=255)
                                                    tctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGtttagatcgg                                                >  2:104559/1‑93 (MQ=255)
                                                     ctactaATTTGCTCTCTTTCTAGTTGCCTTTTCGAGCCTACAAAAGCAACTTCATCACATTTACTAGCTAAGTCATGAGAATATATAAGTTTTCCGc                                                  <  1:973274/97‑1 (MQ=255)
                                                     ctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGATTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGc                                               <  1:1432261/100‑1 (MQ=255)
                                                     ctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTTATCAGAATATCTAAGTTTTCCGCAGc                                               <  2:2061057/100‑1 (MQ=255)
                                                     ctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGc                                               <  2:875199/100‑1 (MQ=255)
                                                        ataatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCaa                                            >  1:2831229/1‑100 (MQ=255)
                                                           atgtGCTCTCTTTAAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCATATGTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAAt                                         <  1:1090329/97‑1 (MQ=255)
                                                           atTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAAt                                         <  2:3096689/100‑1 (MQ=255)
                                                           atTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAAt                                         <  2:1510047/100‑1 (MQ=255)
                                                            tttGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATc                                        <  2:2555175/100‑1 (MQ=255)
                                                              tGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCaa                                      <  2:3235671/100‑1 (MQ=255)
                                                              tGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCaa                                      <  2:2689242/100‑1 (MQ=255)
                                                              tGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCaa                                      <  2:2349999/100‑1 (MQ=255)
                                                               gCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAg                                     >  2:1247019/1‑100 (MQ=255)
                                                                ctctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAgg                                    <  2:1792393/100‑1 (MQ=255)
                                                                  ctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGc                                  >  2:299167/1‑100 (MQ=255)
                                                                   tctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCg                                 <  1:779279/100‑1 (MQ=255)
                                                                    ctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGt                                >  2:2259109/1‑100 (MQ=255)
                                                                    ctTTCAAGTTGCATTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATTAGAATATATAAGTTTTCCGCAGCAAAAATAAAgg                                    >  1:708105/1‑96 (MQ=255)
                                                                         aaGTTGCTTTTCTAGACTCGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAt                          >  1:44721/1‑100 (MQ=255)
                                                                             tGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAtttt                       <  2:2605725/100‑1 (MQ=255)
                                                                             tGCCTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAtttt                       >  2:655942/1‑100 (MQ=255)
                                                                                cTTTTCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtct                    >  2:2327203/1‑100 (MQ=255)
                                                                                 ttttCGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTtctc                   <  2:488210/100‑1 (MQ=255)
                                                                                     cGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATg               <  1:2327186/100‑1 (MQ=255)
                                                                                     cGAGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTCGCCATTCTTCTCAAt                >  2:630967/1‑99 (MQ=255)
                                                                                       aGTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCt             >  2:1292438/1‑100 (MQ=255)
                                                                                        gTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAttttttttttttttt            >  1:1464438/1‑90 (MQ=255)
                                                                                        gTCTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCt             >  1:632291/1‑99 (MQ=255)
                                                                                         tCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTc           <  1:1353564/100‑1 (MQ=255)
                                                                                          cTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt          <  2:793375/100‑1 (MQ=255)
                                                                                          cTAGAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt          >  2:1527068/1‑100 (MQ=255)
                                                                                           tAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTa         >  1:204300/1‑100 (MQ=255)
                                                                                           tAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTa         >  1:73720/1‑100 (MQ=255)
                                                                                            aGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCGGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCt          <  1:2222367/98‑1 (MQ=255)
                                                                                             gAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat       >  1:1143133/1‑100 (MQ=255)
                                                                                             gAAAAGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAat       >  1:985415/1‑100 (MQ=255)
                                                                                              aaaaGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCTCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAata      <  1:1384455/100‑1 (MQ=255)
                                                                                              aaaaGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAata      <  1:285385/100‑1 (MQ=255)
                                                                                              aaaaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAata      <  1:1726618/100‑1 (MQ=255)
                                                                                               aaaGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAata      <  2:2863194/99‑1 (MQ=255)
                                                                                               aaaGCAACTTCATGAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCTGCAGCCAAAATCTAGGGTGTAGCCATTTTTCTCAATGCTTCTAatat     <  1:840413/100‑1 (MQ=255)
                                                                                               aaaGCAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAatat     >  1:2211291/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
TTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACT  >  NC_000913/1197241‑1197438

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: