Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,707,321 A→T 19.8% intergenic (‑197/+1) lepA ← / ← rseC back‑translocating elongation factor EF4, GTPase/SoxR iron‑sulfur cluster reduction factor component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,707,3210AT19.8% 370.9 / 91.0 162intergenic (‑197/+1)lepA/rseCback‑translocating elongation factor EF4, GTPase/SoxR iron‑sulfur cluster reduction factor component
Reads supporting (aligned to +/- strand):  ref base A (74/56);  new base T (9/23);  total (83/79)
Fisher's exact test for biased strand distribution p-value = 5.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

Only 100 of 173 total aligned reads displayed.

CACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAA  >  NC_000913/2707222‑2707420
                                                                                                   |                                                                                                   
cactacaACATTGAGGCGAGGAAATGTTCATAGCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGTAAAATGCAACAAAGCCAGTGaaa                                                                                                     >  1:3399652/1‑100 (MQ=255)
 actacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGaaa                                                                                                     <  2:1445056/99‑1 (MQ=255)
  ctacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGaaagc                                                                                                   >  1:2358893/1‑98 (MQ=255)
  ctacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATc                                                                                                   >  2:1428065/1‑100 (MQ=255)
    acaacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCAc                                                                                                 >  1:1414292/1‑100 (MQ=255)
    acaacaTTGAGGCGAGGAAATGTTCATACCGTATGGAATGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCAc                                                                                                 >  1:1866015/1‑100 (MQ=255)
      aacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTg                                                                                               <  2:720012/100‑1 (MQ=255)
      aacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTg                                                                                               >  1:2668183/1‑100 (MQ=255)
      aacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTg                                                                                               <  1:2190751/100‑1 (MQ=255)
        caTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGc                                                                                             <  1:3751893/100‑1 (MQ=255)
         aTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCt                                                                                            >  1:2170227/1‑100 (MQ=255)
          ttGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCtt                                                                                           <  2:1871918/100‑2 (MQ=255)
           tGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCttt                                                                                          <  1:1363999/100‑3 (MQ=255)
           tGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCttt                                                                                          <  2:3384322/100‑3 (MQ=255)
            gAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttg                                                                                         <  2:3570544/100‑4 (MQ=255)
            gAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCttt                                                                                          <  2:874153/99‑3 (MQ=255)
             aGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgt                                                                                        <  2:289239/100‑5 (MQ=255)
             aGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGCAAACGTCCTCGCATTTGTTATGCAAAATGCAACAACGCCAGTCATTTCACTGGCttt                                                                                          <  1:3287929/98‑3 (MQ=255)
              ggCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCGAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGt                                                                                       >  1:2141669/1‑100 (MQ=255)
              ggCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgtt                                                                                       <  2:3015067/100‑6 (MQ=255)
              ggCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgtt                                                                                       <  1:129924/100‑6 (MQ=255)
              ggCGAGGAAATGTTCATAACGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgtt                                                                                       <  1:3095900/100‑6 (MQ=255)
               gCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgttg                                                                                      <  2:2119867/100‑7 (MQ=39)
               gCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgttg                                                                                      <  2:3713289/100‑7 (MQ=39)
               gCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTc                                                                                      <  2:2170242/100‑1 (MQ=255)
               gCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTc                                                                                      >  2:1844973/1‑100 (MQ=255)
               gCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGAATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTc                                                                                      <  1:980416/100‑1 (MQ=255)
                cGAGGAAATGTTCATACCGTATGGATTGTGGTATGTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgttgc                                                                                     <  1:2046240/100‑8 (MQ=38)
                cGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgttgc                                                                                     <  1:2863734/100‑8 (MQ=38)
                  aGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgttgcat                                                                                   <  2:3687859/100‑10 (MQ=38)
                   ggAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTcc                                                                                  <  1:3495352/100‑1 (MQ=255)
                   ggAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTcc                                                                                  <  2:113749/100‑1 (MQ=255)
                      aaTGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGGGAAATCACTGGCTCGCGTCTTCCGaa                                                                               <  2:2358908/100‑1 (MQ=255)
                          ttCATACCGTATGGACTACGGTACCTGGAAACGTCCTCGCATTTGTTATTCATCATGCAACAAAGCCAGTGAAATCACTCGCTCACCACTTACGAAGa                                                                             >  1:2528529/1‑98 (MQ=255)
                           tCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGt                                                                          <  2:1841636/100‑1 (MQ=255)
                           tCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGt                                                                          >  1:1967724/1‑100 (MQ=255)
                            cATACCGTATGGATTGTGGTATCTTGAAACGTCCTCGCATTTGTTATGCAAAATGCAACGAATCCAGTGAAATCACTGGCTCTCGTCTTCCGAAGATGtt                                                                         <  1:366937/100‑1 (MQ=255)
                              tACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTc                                                                       <  1:2842025/100‑1 (MQ=255)
                                 cGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCaaa                                                                    >  2:464134/1‑100 (MQ=255)
                                    aTGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCg                                                                 >  2:716954/1‑100 (MQ=255)
                                     tGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGc                                                                >  1:1016596/1‑100 (MQ=255)
                                       gATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAGATCggaa                                                              >  1:334146/1‑97 (MQ=255)
                                          tgtgGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAg                                                           >  1:2799526/1‑100 (MQ=255)
                                           gtgGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAgg                                                          >  1:2865075/1‑100 (MQ=255)
                                             ggTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGcc                                                        >  1:1319020/1‑100 (MQ=255)
                                             ggTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGcc                                                        >  2:2444002/1‑100 (MQ=255)
                                              gTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCa                                                       >  1:3060641/1‑100 (MQ=255)
                                              gTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCa                                                       >  1:3017916/1‑100 (MQ=255)
                                              gTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCa                                                       <  2:2938207/100‑1 (MQ=255)
                                                    ggAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAggcggc                                                 >  1:3257637/1‑100 (MQ=255)
                                                     gAAACGTCCTCGCATTTGTTATGCAAAATGCAAAAAAACCAGTGAAATCACTGGCTAGCGTATTCCGAAGATGTTTCAAATCGCACAAAGCCAgg                                                     >  1:171160/1‑95 (MQ=255)
                                                       aaCGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAAt                                              <  2:2886374/100‑1 (MQ=255)
                                                       aaCGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAAt                                              >  1:2295615/1‑100 (MQ=255)
                                                         cGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGc                                            >  1:2989988/1‑100 (MQ=255)
                                                         cGTCCTCGAATTTGTTATGCAAAATGAAAAAAAGCCAGTGAAATCACTGGCTCGCGTATTCCGAAGATGTTTCAAATCGCACAAGGCAAgacggtga                                               >  1:734855/1‑94 (MQ=255)
                                                          gTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGcc                                           >  1:1988902/1‑100 (MQ=255)
                                                          gTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGcc                                           <  2:158649/100‑1 (MQ=255)
                                                           tCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCa                                          <  2:361007/100‑1 (MQ=255)
                                                            ccTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCAc                                         <  2:3414141/100‑1 (MQ=255)
                                                                gCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCtt                                     >  2:1268632/1‑100 (MQ=255)
                                                                 cATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTa                                    >  2:2791379/1‑100 (MQ=255)
                                                                 cATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTACGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTa                                    >  1:1483408/1‑100 (MQ=255)
                                                                   tttGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTag                                   <  2:2111382/99‑2 (MQ=255)
                                                                   tttGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAg                                  <  1:3543318/100‑1 (MQ=255)
                                                                   tttGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAg                                  >  2:556866/1‑100 (MQ=255)
                                                                    ttGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAga                                 <  1:1428057/100‑1 (MQ=255)
                                                                     tGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgat                                <  2:2087771/100‑1 (MQ=255)
                                                                     tGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgat                                >  1:771997/1‑100 (MQ=255)
                                                                     tGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgat                                >  2:1042745/1‑100 (MQ=255)
                                                                     tGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgat                                >  2:2801024/1‑100 (MQ=255)
                                                                     tGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgat                                >  1:3573981/1‑100 (MQ=255)
                                                                       ttATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgatga                              >  2:2336648/1‑100 (MQ=255)
                                                                           gCAAAATGCAACAAAGCCAGTGAAATCACTGGCGCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCCGGCGGCAATGCCACGCTTAAGATGATCgg                          >  1:2337519/1‑100 (MQ=255)
                                                                            cAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGt                         >  2:331832/1‑100 (MQ=255)
                                                                               aaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCTACGATTAAGATGATCGGtt                        >  1:907473/1‑98 (MQ=255)
                                                                               aaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGc                      <  2:304589/100‑1 (MQ=255)
                                                                                aTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGcc                     <  2:3139663/100‑1 (MQ=255)
                                                                                 tGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCa                    <  2:3796846/100‑1 (MQ=255)
                                                                                   caacaaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCAtt                  >  1:2786084/1‑100 (MQ=255)
                                                                                   caacaaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCAtt                  <  2:1711974/100‑1 (MQ=255)
                                                                                    aacaaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTc                 <  2:39559/100‑1 (MQ=255)
                                                                                    aacaaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTc                 >  2:1665371/1‑100 (MQ=255)
                                                                                     cgcgaGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAAGCCAGGCGGCAATGCAAAGCTTAAGACGATCGGTTGCCAtt                  >  1:2682589/5‑98 (MQ=255)
                                                                                      caaAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTg               <  2:255327/100‑1 (MQ=255)
                                                                                       cgaGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGc              >  2:1583350/3‑100 (MQ=255)
                                                                                       cgaGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGc              >  2:2235863/3‑100 (MQ=255)
                                                                                       cgaGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGc              >  2:393054/3‑100 (MQ=255)
                                                                                       aaaGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGc              >  2:1490102/1‑100 (MQ=255)
                                                                                       aaaGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGc              >  2:3362275/1‑100 (MQ=255)
                                                                                         aGCCAGTGATTTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGccc            >  2:3148421/1‑100 (MQ=255)
                                                                                           ccAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGGTGATCGGTTGCCATTCTGCCCgg          <  1:526845/100‑1 (MQ=255)
                                                                                            cAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCggg         <  1:2443985/100‑1 (MQ=255)
                                                                                            cAGTGAAATCACTGGCTCGCGTCGTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCgg          <  1:625339/99‑1 (MQ=255)
                                                                                             aGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgc        <  2:1186589/100‑1 (MQ=255)
                                                                                             aGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgc        <  1:1091371/100‑1 (MQ=255)
                                                                                               tGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgcag      <  2:1283389/100‑1 (MQ=255)
                                                                                                gAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCCCAAGGCCAGGCGGCAATGCCACGCTTAAGATGATGGGTTGCCAGTCTGCCCGGgcagc     >  1:852391/1‑100 (MQ=255)
                                                                                                gAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgcagc     <  2:1646688/100‑1 (MQ=255)
                                                                                                 aaaTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGgcagc     >  2:771807/1‑99 (MQ=255)
                                                                                                   aTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCaaa  >  1:1196457/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAA  >  NC_000913/2707222‑2707420

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: