Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 926,168 C→T 18.0% intergenic (‑197/+57) serW ← / ← infA tRNA‑Ser/translation initiation factor IF‑1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913926,1680CT18.0% 489.0 / 81.3 211intergenic (‑197/+57)serW/infAtRNA‑Ser/translation initiation factor IF‑1
Reads supporting (aligned to +/- strand):  ref base C (80/92);  new base T (24/14);  total (105/106)
Fisher's exact test for biased strand distribution p-value = 7.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

TTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTAC  >  NC_000913/926069‑926267
                                                                                                   |                                                                                                   
 tGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGcc                                                                                                    <  1:4064830/100‑1 (MQ=255)
   cccATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCgg                                                                                                  <  1:2676366/98‑1 (MQ=255)
   cTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCgg                                                                                                  <  1:13787/100‑1 (MQ=255)
   cTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCgg                                                                                                  >  2:2385575/1‑100 (MQ=255)
    tCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGt                                                                                                 <  2:1294150/100‑1 (MQ=255)
        aCTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaaa                                                                                             >  1:997365/1‑100 (MQ=255)
         cTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAAc                                                                                            >  1:1263301/1‑100 (MQ=255)
          tGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAAcc                                                                                           <  2:2125094/100‑1 (MQ=255)
          tGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAAc                                                                                            >  1:3383328/1‑99 (MQ=255)
           gTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCg                                                                                          >  2:1064309/1‑100 (MQ=255)
           gTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCg                                                                                          <  1:327167/100‑1 (MQ=255)
            tAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCgg                                                                                         <  1:962519/100‑2 (MQ=255)
            tAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGa                                                                                         <  1:1659588/100‑1 (MQ=255)
             aaaCAAGTTGTTCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGAc                                                                                        <  1:2924292/100‑1 (MQ=255)
              aaCAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATAGATATCACAAAATAAAAAGGCCGGTTAAACCGAcc                                                                                       >  1:772485/1‑100 (MQ=255)
                  aGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttt                                                                                   <  2:166522/100‑1 (MQ=255)
                  aGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCtttt                                                                                   >  2:3804523/1‑100 (MQ=255)
                  aGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCttt                                                                                    >  1:245406/1‑99 (MQ=255)
                   gTTGTGCAGTATATCTACATCGAGACAAGTTACGTACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCtttt                                                                                   <  1:2881827/99‑1 (MQ=255)
                   gTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCttttt                                                                                  <  2:2320713/100‑2 (MQ=39)
                   gTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCttttt                                                                                  <  2:2564679/100‑2 (MQ=39)
                   gTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCttttt                                                                                  <  2:259975/100‑2 (MQ=39)
                    ttGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTan                                                                                 <  2:2773763/100‑2 (MQ=255)
                    ttGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTAc                                                                                 <  2:1606332/100‑1 (MQ=255)
                      gtgCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttatt                                                                               <  2:3135524/100‑5 (MQ=38)
                       tgCAGTATATCTACATGGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttattt                                                                              <  1:463699/100‑6 (MQ=38)
                       tgCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCg                                                                              <  1:3896001/100‑1 (MQ=255)
                        gCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCga                                                                             >  2:3896030/1‑99 (MQ=255)
                         cAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttattttg                                                                            <  2:968186/100‑8 (MQ=37)
                          aGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttc                                                                           <  1:1426062/100‑1 (MQ=255)
                           gTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttattttgtg                                                                          <  2:1442434/100‑10 (MQ=37)
                           gTATATCTACAACGAGACAAGTTACGCACCTATACTTCCAAAGTAATTCATACATATCACAATATAAAAAGGTCGGTTTAACCGGCCtttttattt                                                                              <  1:2218054/96‑6 (MQ=25)
                             atatCTACATCGAGACAAGTTACGGACTCATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGCCCTTTTAATCGttctat                                                                        >  1:2801012/1‑98 (MQ=255)
                              tatCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttc                                                                       <  2:859137/100‑1 (MQ=255)
                              tatCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttc                                                                       <  1:3030778/100‑1 (MQ=255)
                              tatCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttc                                                                       <  1:1737374/100‑1 (MQ=255)
                                tCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTtctc                                                                     >  2:3985173/1‑100 (MQ=255)
                                 cTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTctct                                                                    <  2:2647353/100‑1 (MQ=255)
                                  tACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCtt                                                                   >  2:3129339/1‑100 (MQ=255)
                                  tACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCtt                                                                   >  2:2146866/1‑100 (MQ=255)
                                    cATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCCTTTACTCGTTCTTTCTCTTCg                                                                 >  1:2508221/1‑100 (MQ=255)
                                      tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGcc                                                               >  2:2765726/1‑100 (MQ=255)
                                        gagaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCa                                                             >  1:3385485/1‑100 (MQ=255)
                                          gaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAGGGGGGGTTAAACCGACCTTTTACTCGTTCTTGCTCTTCGCCCGTc                                                           >  2:3512477/1‑100 (MQ=255)
                                           aCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCa                                                          >  2:319337/1‑100 (MQ=255)
                                              aGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGc                                                       >  1:1250677/1‑100 (MQ=255)
                                              aGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGc                                                       <  2:2854234/100‑1 (MQ=255)
                                                 tACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGt                                                    >  2:423471/1‑100 (MQ=255)
                                                  aCGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTa                                                   <  2:1537140/100‑1 (MQ=255)
                                                     gACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTaaaa                                                <  2:365934/100‑1 (MQ=255)
                                                      aCTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAAc                                               <  1:380034/100‑1 (MQ=255)
                                                       cTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACa                                              >  1:1491744/1‑100 (MQ=255)
                                                         tataCTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAAt                                            >  2:2258343/1‑100 (MQ=255)
                                                            aCTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAg                                         <  2:1920400/100‑1 (MQ=255)
                                                             cTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGc                                        <  2:36716/100‑1 (MQ=255)
                                                             cTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGc                                        <  1:2378823/100‑1 (MQ=255)
                                                             cTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGc                                        <  1:2555770/100‑1 (MQ=255)
                                                              ttCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCg                                       <  1:3582927/100‑1 (MQ=255)
                                                               tCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGa                                      >  1:554251/1‑100 (MQ=255)
                                                                  aaaGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTa                                   <  2:1704863/100‑1 (MQ=255)
                                                                  aaaGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTa                                   >  1:4282079/1‑100 (MQ=255)
                                                                     gTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACgg                                <  1:4008789/100‑1 (MQ=255)
                                                                      tACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGa                               <  2:4034572/100‑1 (MQ=255)
                                                                       aCTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGaa                              >  2:587172/1‑100 (MQ=255)
                                                                          tCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGAc                           >  2:3313375/1‑100 (MQ=255)
                                                                          tCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGAc                           <  1:2623257/100‑1 (MQ=255)
                                                                          tCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGAc                           <  1:4330299/100‑1 (MQ=255)
                                                                          tCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGAc                           <  1:1646781/100‑1 (MQ=255)
                                                                             tacataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAAt                        >  2:3147584/1‑100 (MQ=255)
                                                                              acataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATg                       <  1:319336/100‑1 (MQ=255)
                                                                               cataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGc                      >  2:1467291/1‑100 (MQ=255)
                                                                               cataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGc                      <  2:3713533/100‑1 (MQ=255)
                                                                                 tatCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCgg                    <  2:2947727/100‑1 (MQ=255)
                                                                                  atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                   >  2:3619045/1‑100 (MQ=255)
                                                                                  atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                   <  2:497096/100‑1 (MQ=255)
                                                                                  atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                   <  1:4253580/100‑1 (MQ=255)
                                                                                  atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                   >  1:1018709/1‑100 (MQ=255)
                                                                                    cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCt                 >  1:2333038/1‑100 (MQ=255)
                                                                                    cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCt                 <  1:3121563/100‑1 (MQ=255)
                                                                                    cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCt                 <  2:1488433/100‑1 (MQ=255)
                                                                                     gaacgaataaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCtt                >  1:1745077/9‑100 (MQ=38)
                                                                                        aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc             >  1:1726692/1‑100 (MQ=255)
                                                                                         gagtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCt            >  1:493980/5‑100 (MQ=38)
                                                                                          agtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTc           >  2:684315/4‑100 (MQ=38)
                                                                                          agtaaaaGGTCGGTTTAACCGGCCCTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTc           >  2:1187149/4‑100 (MQ=255)
                                                                                          aTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCCGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTc           >  2:241504/1‑100 (MQ=255)
                                                                                          aTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTc           <  2:3546148/100‑1 (MQ=255)
                                                                                          aTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTc           <  2:1782595/100‑1 (MQ=255)
                                                                                           tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCa          <  2:3062705/100‑1 (MQ=255)
                                                                                           tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATAAGGCGATAAAACAACCCGCGACTACGGAAGACAATGCGGCCTTTGCTc           >  2:4072769/1‑99 (MQ=255)
                                                                                           gtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCt            >  1:3859120/3‑98 (MQ=38)
                                                                                           gtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCa          >  2:1830120/3‑100 (MQ=39)
                                                                                           gtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCa          >  1:3141533/3‑100 (MQ=38)
                                                                                           ataaaaGGTCGGTTTAACCGGCCTCTTACTCGTTCTTTCTCTTCGCCCATCAAGCGGTAAAACAATCATCGACTACGGAAGACAATGCGGCCTTTGCTCa          >  1:1702857/3‑100 (MQ=255)
                                                                                            aaaaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAg         >  2:2246086/1‑100 (MQ=255)
                                                                                             aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAgg        >  1:1144183/1‑100 (MQ=255)
                                                                                             aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAgg        >  2:4093566/1‑100 (MQ=255)
                                                                                             aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAgg        >  2:2464541/1‑100 (MQ=255)
                                                                                             aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAgg        >  2:2111925/1‑100 (MQ=255)
                                                                                              aaaGGTCGGTTTAACCGGCCTTGTACGCGTTCTGTCGCTGCGCCCATCAGGCGGGAAAACAATCAGCGACTACGGAAGACAATGCGGCCCTTGCTCAgg        >  1:3679488/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
TTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTAC  >  NC_000913/926069‑926267

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