Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 926,174 A→T 18.8% intergenic (‑203/+51) 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,1740AT18.8% 461.3 / 98.1 203intergenic (‑203/+51)serW/infAtRNA‑Ser/translation initiation factor IF‑1
Reads supporting (aligned to +/- strand):  ref base A (77/87);  new base T (24/14);  total (101/102)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

Only 100 of 212 total aligned reads displayed.

ATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTC  >  NC_000913/926075‑926273
                                                                                                   |                                                                                                   
 tACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaa                                                                                                    >  2:380035/1‑100 (MQ=255)
  aCTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTaaa                                                                                                   >  1:997365/1‑100 (MQ=255)
   cTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAAc                                                                                                  >  2:353727/1‑100 (MQ=255)
   cTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAAc                                                                                                  >  1:1263301/1‑100 (MQ=255)
    tGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAAcc                                                                                                 <  2:2125094/100‑1 (MQ=255)
     gTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCg                                                                                                >  2:1064309/1‑100 (MQ=255)
      tAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCgg                                                                                               <  1:962519/100‑2 (MQ=255)
       aaaCAAGTTGTTCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGAc                                                                                              <  1:2924292/100‑1 (MQ=255)
        aaCAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATAGATATCACAAAATAAAAAGGCCGGTTAAACCGAcc                                                                                             >  1:772485/1‑100 (MQ=255)
        aaCAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGAcc                                                                                             >  1:4330428/1‑100 (MQ=255)
           aaGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCttt                                                                                          <  1:4120776/100‑1 (MQ=255)
           aaGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCttt                                                                                          <  1:3238409/100‑1 (MQ=255)
           aaGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCttt                                                                                          >  1:3192435/1‑100 (MQ=255)
            aGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttt                                                                                         <  2:166522/100‑1 (MQ=255)
            aGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCtttt                                                                                         >  2:2250656/1‑100 (MQ=255)
             gTTGTGCAGTATATCTACATCGAGACAAGTTACGTACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCtttt                                                                                         <  1:2881827/99‑1 (MQ=255)
              ttGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTAc                                                                                       <  2:1606332/100‑1 (MQ=255)
              ttGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTAc                                                                                       <  2:2545231/100‑1 (MQ=255)
                gtgCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttatt                                                                                     <  2:3135524/100‑5 (MQ=38)
                 tgCAGTATATGTACATCGAGACAAGTTACCGACATAAACTTTCAAAGTACTTCATACATATCACATAATAAAAAGGTCCGTTTTACCGGCCttttgattt                                                                                    <  1:4300511/100‑6 (MQ=255)
                 tgCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCg                                                                                    <  1:3896001/100‑1 (MQ=255)
                   cAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttattttg                                                                                  <  2:968186/100‑8 (MQ=37)
                     gTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttattttgtg                                                                                <  1:3384698/100‑10 (MQ=37)
                     gTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGTCGGTTTAACCGGCCtttttattttgtg                                                                                <  2:1442434/100‑10 (MQ=37)
                     gTATATCTACAACGAGACAAGTTACGCACCTATACTTCCAAAGTAATTCATACATATCACAATATAAAAAGGTCGGTTTAACCGGCCtttttattt                                                                                    <  1:2218054/96‑6 (MQ=25)
                      tataTCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctt                                                                               <  2:3669591/100‑1 (MQ=255)
                       atatCTACATCGAGACAAGTTACGGACTCATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGCCCTTTTAATCGttctat                                                                              >  1:2801012/1‑98 (MQ=255)
                        tatCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttc                                                                             <  1:3030778/100‑1 (MQ=255)
                        tatCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttc                                                                             <  2:859137/100‑1 (MQ=255)
                        tatCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttc                                                                             <  1:1737374/100‑1 (MQ=255)
                         atCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttct                                                                            >  1:365933/1‑100 (MQ=255)
                         atCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttctttct                                                                            <  2:3987178/100‑1 (MQ=255)
                          tCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAACAAGGCCGGTTAAACCGATCTTTTACTCGttctttct                                                                            >  1:612933/1‑99 (MQ=255)
                          tCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTtctc                                                                           >  2:2239963/1‑100 (MQ=255)
                          tCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTtctc                                                                           >  2:2871597/1‑100 (MQ=255)
                            tACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCtt                                                                         >  2:2146866/1‑100 (MQ=255)
                            tACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCtt                                                                         >  2:3129339/1‑100 (MQ=255)
                             aCATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTAATACATATCACAAAATACAAAGGCCGGTTAAACCGACCTTTTACTCGt                                                                                   >  1:3415573/1‑89 (MQ=255)
                              cATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCCTTTACTCGTTCTTTCTCTTCg                                                                       >  1:2508221/1‑100 (MQ=255)
                                 cGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGccc                                                                    >  2:229489/1‑100 (MQ=255)
                                  gagaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCa                                                                   >  1:3364539/1‑100 (MQ=255)
                                     aCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCa                                                                >  2:319337/1‑100 (MQ=255)
                                        aGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGc                                                             >  1:1250677/1‑100 (MQ=255)
                                         gTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCg                                                            >  1:1402898/1‑100 (MQ=255)
                                         gTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCg                                                            <  2:3731707/100‑1 (MQ=255)
                                           tACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGt                                                          >  2:423471/1‑100 (MQ=255)
                                            aCGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTa                                                         <  2:2744329/100‑1 (MQ=255)
                                            aCGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTa                                                         <  2:1537140/100‑1 (MQ=255)
                                             cGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTaa                                                        >  2:2847199/1‑100 (MQ=255)
                                                aCTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAAc                                                     <  1:380034/100‑1 (MQ=255)
                                                 cTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACa                                                    >  1:3139359/1‑100 (MQ=255)
                                                      aCTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAg                                               <  2:1920400/100‑1 (MQ=255)
                                                       cTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGc                                              <  1:2555770/100‑1 (MQ=255)
                                                        ttCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCg                                             <  1:260078/100‑1 (MQ=255)
                                                            aaaGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTa                                         <  2:1704863/100‑1 (MQ=255)
                                                               gTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACgg                                      <  1:4008789/100‑1 (MQ=255)
                                                                 aCTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGaa                                    >  2:587172/1‑100 (MQ=255)
                                                                    tCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGAc                                 >  2:3313375/1‑100 (MQ=255)
                                                                    tCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGAc                                 <  1:1646781/100‑1 (MQ=255)
                                                                       tacataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAAt                              >  2:3147584/1‑100 (MQ=255)
                                                                         cataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGc                            >  2:1467291/1‑100 (MQ=255)
                                                                         cataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGc                            >  1:489833/1‑100 (MQ=255)
                                                                            atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                         <  1:4253580/100‑1 (MQ=255)
                                                                            atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                         >  2:3619045/1‑100 (MQ=255)
                                                                             aagaacgagtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGcc                        >  2:4116886/11‑100 (MQ=37)
                                                                              cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCt                       >  1:2333038/1‑100 (MQ=255)
                                                                              cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCt                       <  2:1488433/100‑1 (MQ=255)
                                                                              cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCt                       >  2:277361/1‑100 (MQ=255)
                                                                               gaacgaataaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCtt                      >  1:1745077/9‑100 (MQ=38)
                                                                               acaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCtt                      >  1:2634582/1‑100 (MQ=255)
                                                                                caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCttt                     <  2:2642719/100‑1 (MQ=255)
                                                                                caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCttt                     <  1:1580483/100‑1 (MQ=255)
                                                                                  cgagtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                   >  1:1482117/6‑100 (MQ=38)
                                                                                  cgagtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                   >  1:1700461/6‑100 (MQ=38)
                                                                                  cgagtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                   >  1:2441513/6‑100 (MQ=38)
                                                                                  aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                   >  1:1726692/1‑100 (MQ=255)
                                                                                  aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                   <  2:3850287/100‑1 (MQ=255)
                                                                                    agtaaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTc                 >  2:684315/4‑100 (MQ=38)
                                                                                    agtaaaaGGTCGGTTTAACCGGCCCTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTc                 >  2:1187149/4‑100 (MQ=255)
                                                                                    aTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCCGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTc                 >  2:241504/1‑100 (MQ=255)
                                                                                     tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCa                <  2:3062705/100‑1 (MQ=255)
                                                                                     tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATAAGGCGATAAAACAACCCGCGACTACGGAAGACAATGCGGCCTTTGCTc                 >  2:4072769/1‑99 (MQ=255)
                                                                                       aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAgg              >  2:2111925/1‑100 (MQ=255)
                                                                                       aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAgg              >  1:1144183/1‑100 (MQ=255)
                                                                                       aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAgg              >  2:1818209/1‑100 (MQ=255)
                                                                                       aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAgg              >  2:2464541/1‑100 (MQ=255)
                                                                                       aaaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAgg              <  2:1726699/100‑1 (MQ=255)
                                                                                       aaaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAAACAATGCGGCCTTTGCTCAgg              >  2:2570343/1‑100 (MQ=255)
                                                                                        aaaGGTCGGTTTAACCGGCCTTGTACGCGTTCTGTCGCTGCGCCCATCAGGCGGGAAAACAATCAGCGACTACGGAAGACAATGCGGCCCTTGCTCAgg              >  1:3679488/1‑99 (MQ=255)
                                                                                        aaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGt             <  2:3492514/100‑1 (MQ=255)
                                                                                         aaGGCCGGTTAAACCGACTTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                   <  2:3675976/93‑1 (MQ=255)
                                                                                         aaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTc            >  2:3971902/1‑100 (MQ=255)
                                                                                                gTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACggg     <  2:2333048/100‑1 (MQ=255)
                                                                                                gTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACggg     >  1:1883362/1‑100 (MQ=255)
                                                                                                gTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgg      <  1:1033495/99‑1 (MQ=255)
                                                                                                 ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgggg    <  1:4259203/100‑1 (MQ=255)
                                                                                                 ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgggg    <  1:540459/100‑1 (MQ=255)
                                                                                                 ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgggg    <  1:1006700/100‑1 (MQ=255)
                                                                                                 ntAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACagg     >  2:1374786/2‑99 (MQ=255)
                                                                                                   aaaCCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTc  >  2:457758/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTC  >  NC_000913/926075‑926273

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: