Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,707,321 A→T 6.3% intergenic (‑197/+1) lepA ← / ← rseC elongation factor 4/protein RseC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,707,3210AT 6.3% 159.7 / 2.1 62intergenic (‑197/+1)lepA/rseCelongation factor 4/protein RseC
Reads supporting (aligned to +/- strand):  ref base A (33/24);  new base T (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

GTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATAC  >  NC_000913/2707179‑2707461
                                                                                                                                              |                                                                                                                                            
gTTCTGAGCTTGACGTGTGTAGTAATAGAAGCGCCGCACTATACACTTCAACATTGAGGCTAGAAAATGTTCAGACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGaaa                                                                                                                                              <  2:924157/143‑1 (MQ=255)
 ttCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAAt                                                                                                                                             >  1:189408/1‑143 (MQ=255)
 ttCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGGGGTATCCGGAAACGACCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAAt                                                                                                                                             >  2:336379/1‑143 (MQ=255)
 ttCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTAATGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAAt                                                                                                                                             >  2:663670/1‑143 (MQ=255)
         ttAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCt                                                                                                                                     >  1:390754/1‑143 (MQ=255)
          tAACGTCTGTTTTAATAGAAACGCCGCATTCTACACTACAACGGTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtt                                                                                                                                    <  1:69100/143‑2 (MQ=255)
          tAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtt                                                                                                                                    <  2:370419/143‑2 (MQ=255)
             cggccgTATCAATAGAAACGCCGCATTCTACACTACAACGTTGAGGCGAGGAAATGGTCATACCGTATGGACAGGGGGAGGTGGTCACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTcgcg                                                                                                                                 <  2:493499/138‑1 (MQ=255)
                 tGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCtt                                                                                                                             >  1:702306/1‑143 (MQ=255)
                   tATTAATAGAAACGCCGCATTCTACACTACAACCTACAGCCGAGGAAAAGTTCATACCGTATGGATCGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAACCCAGTGAAATCACTGCCTCACGGCTTcc                                                                                                                           >  2:1081480/1‑143 (MQ=255)
                   tATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAAGGCAACAAAGCCAGTGAAATAACTGGCTCGAGTCTTcc                                                                                                                           >  1:343608/1‑143 (MQ=255)
                     ttAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGa                                                                                                                         >  2:697107/1‑143 (MQ=255)
                     ttAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGACCCGCGTCTTCCGa                                                                                                                         >  2:432866/1‑143 (MQ=255)
                     ctAATAGAAACGCCGCATTCTAGACGACAACAGGGAGGCGACGAAATGTTCATACCGACAGGATCGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGa                                                                                                                         <  1:1099907/142‑1 (MQ=255)
                          aGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATg                                                                                                                    >  1:403494/1‑143 (MQ=255)
                           gAAACGCCGCATTCTACATTACAACATTGAGGCGTGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCTAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGt                                                                                                                   <  2:791469/143‑1 (MQ=255)
                            aaaCACCGCATTCTACACTAAAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCCAAATGCAACAAAGCGAGTGACATCACTGGCTCGCGTAGTCCGAAGCTGtt                                                                                                                  >  2:25256/1‑143 (MQ=255)
                             caCGCCGCTTTCTTCACTCCCAGAGTGAGGCGTGGAAATGTTCCTACCGTATGGATTGTGGTATCTGGACACGTCCGCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGttt                                                                                                                 <  2:271828/142‑1 (MQ=255)
                             aaCGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTAGCATATGTTATGAAAAATGCAACAAAGCCAGTGAAATCACGGGCTCGCGTCGTACGAAAAGGttg                                                                                                                 >  2:1036614/1‑142 (MQ=255)
                              aCGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTc                                                                                                                <  2:590916/143‑1 (MQ=255)
                              aCGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGGTTc                                                                                                                >  1:870537/1‑143 (MQ=255)
                                     aTTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGc                                                                                                         <  2:395286/143‑1 (MQ=255)
                                      ttCTACACTACAACATTGGGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGACGAGGATTCAAATCGcg                                                                                                        >  2:574994/1‑142 (MQ=255)
                                      ttCTACACTACAACAGTGAGGCGAGGAAATGTTCATACCGTATGGATTGGGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGca                                                                                                        <  1:825496/143‑1 (MQ=255)
                                          accccacaacaTTGCGGCGAATAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTTTTATCTAAAATGAAACATATCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAg                                                                                                    <  2:1078038/138‑1 (MQ=255)
                                          acactacaACATTGAGGCGCGACAATGACCATACCTTATGGATTGTGGTATCTGGAAACGTCATCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAg                                                                                                    <  1:559499/143‑1 (MQ=255)
                                            actacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGc                                                                                                  <  2:501332/143‑1 (MQ=255)
                                            actacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTCCCGAAGATGCCTCAAATCGCACAACGc                                                                                                  >  2:954241/1‑143 (MQ=255)
                                              tacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGGCAGTGAATTCACTGTCTCGCGTCTTCCGAAGATGTTTCAAATCGCACCAGGccg                                                                                                >  1:465619/1‑142 (MQ=255)
                                              tacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCCCGCGTCTTCCGAAGATGCTTCAAAACACACAAGGccc                                                                                                >  2:291322/1‑142 (MQ=255)
                                               acaacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAg                                                                                               <  2:270349/143‑1 (MQ=255)
                                                caacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAgg                                                                                              >  1:927403/1‑143 (MQ=255)
                                                 aacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTGTGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCGGAAGATGTTTCAAATCTCACAAGGCCAggc                                                                                             >  1:676168/1‑143 (MQ=255)
                                                   cagggAGGCGAGGAAGTGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAggcgg                                                                                           <  1:679548/139‑1 (MQ=255)
                                                                    gTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCCGTGAAATCACTGGCTCGCGTATTCCCAAGAGGTTTCAAATCCCACAAGGCCCGGCGCCCCTCCCACGCCTCAga                                                                          >  2:845004/1‑143 (MQ=255)
                                                                      tCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgatg                                                                        >  1:699628/1‑143 (MQ=255)
                                                                        aTACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgatgat                                                                      <  1:1032705/143‑1 (MQ=255)
                                                                         tACCATATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGGTc                                                                     >  2:186303/1‑143 (MQ=255)
                                                                                    ttGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTc                                                          <  1:585244/143‑1 (MQ=255)
                                                                                     tgtgGTATCTGGAAACGTCCTCGTATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCt                                                         <  2:309511/143‑1 (MQ=255)
                                                                                     tgtgGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGTCACGCTTAAGATGATCGGTTGCCATTCt                                                         <  2:189408/143‑1 (MQ=255)
                                                                                     tgtgGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCt                                                         <  1:663670/143‑1 (MQ=255)
                                                                                     tcgggTATCCGGACACTTCCCCGCATTTGTTATGCCAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCATGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCt                                                         <  1:217984/140‑1 (MQ=255)
                                                                                           aTCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCCATGCCACGCTTAAGCTGATCGGTTGCCCTTCTGCCCgg                                                   >  1:539273/1‑143 (MQ=255)
                                                                                             cTGGAAACGTCCTCGCCTTGGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgc                                                 <  2:1067053/143‑1 (MQ=255)
                                                                                               ggAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgcag                                               <  2:600491/143‑1 (MQ=255)
                                                                                                   aCGCCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCaaa                                           <  1:244420/143‑1 (MQ=255)
                                                                                                      tCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGATTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACtt                                        >  2:248991/1‑143 (MQ=255)
                                                                                                         tCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACACGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTAcg                                     >  2:891718/1‑143 (MQ=255)
                                                                                                           gCATTTGTTATGCAAAATGCAACAAAACCAGTGAAATCACTGGCTCGAGTCTCCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCGACGCTTAAGATGATCAGTTGGCATTCCGCCCGGGCACCAAAccaacacg                                   >  2:209083/1‑136 (MQ=255)
                                                                                                            cATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGa                                  >  2:891310/1‑143 (MQ=255)
                                                                                                            cATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGa                                  <  1:328756/143‑1 (MQ=255)
                                                                                                                tGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCCGGCGGCAACGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAACCCGACGCGAGTAg                              >  2:654334/1‑143 (MQ=255)
                                                                                                                    aTGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGACCGGGCAGCAAACTTACGCGAGTAGCcgc                          >  1:1058240/1‑143 (MQ=255)
                                                                                                                     tGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCACACCTACGCGACTAGCcgcg                         >  2:876234/1‑143 (MQ=255)
                                                                                                                     tGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCcgcg                         >  1:652717/1‑143 (MQ=255)
                                                                                                                     tGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAACAAACTTACGCGAGTAGCcgcg                         >  1:489894/1‑143 (MQ=255)
                                                                                                                                acaaCGCCAGTGAACACTCTGGCTCGCGTCTGCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGa              <  2:1009567/143‑1 (MQ=255)
                                                                                                                                 caaAGCCAGTGAAATCGCCGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGaa             <  2:511305/143‑1 (MQ=255)
                                                                                                                                  aaaGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGCTGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCCCAATCAGGAAc            >  1:709632/1‑143 (MQ=255)
                                                                                                                                   gaGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCACTGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTCACGCGAGGAGCCGCGCGCAACCAGGAAcc           >  2:303857/2‑143 (MQ=255)
                                                                                                                                        aGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGCACCCCGCCa      >  2:689788/1‑143 (MQ=255)
                                                                                                                                           gAAATCACTGTCTCGCGTCTTCCGAAGATGTTTCAAATCGCATAAGGCCAGGCGGCAATGACACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATa   <  1:631570/143‑1 (MQ=255)
                                                                                                                                            acaTCCCTGGCCCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATac  <  1:291322/141‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATAC  >  NC_000913/2707179‑2707461

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