Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,541,456 A→G I75T (ATT→ACT)  cysW ← sulfate/thiosulfate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,541,4560AG99.3% 475.1 / ‑4.9 136I75T (ATT→ACT) cysWsulfate/thiosulfate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (65/70);  total (66/70)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 136 total aligned reads displayed.

TCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGC  >  NC_000913/2541360‑2541551
                                                                                                |                                                                                               
tCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCa                                                                                              <  2:2671289/100‑1 (MQ=255)
tCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGAAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTc                                                                                               >  2:3299725/1‑99 (MQ=255)
 ccAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAg                                                                                             >  2:2437190/1‑100 (MQ=255)
 ccAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAg                                                                                             >  2:2413699/1‑100 (MQ=255)
 ccAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAg                                                                                             >  2:1105235/1‑100 (MQ=255)
  cAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAgc                                                                                            >  1:279931/1‑100 (MQ=255)
  cAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAgc                                                                                            >  1:457072/1‑100 (MQ=255)
   aGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAgcg                                                                                           >  2:843151/1‑100 (MQ=255)
   aGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAgcg                                                                                           >  1:448576/1‑100 (MQ=255)
      aGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGa                                                                                        >  2:2081775/1‑100 (MQ=255)
       gCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGat                                                                                       <  2:1880351/100‑1 (MQ=255)
       gCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGat                                                                                       <  1:2466165/100‑1 (MQ=255)
       gCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGat                                                                                       <  2:2672321/100‑1 (MQ=255)
         gTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatca                                                                                     <  2:1063899/100‑1 (MQ=255)
         gTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatca                                                                                     >  1:1221868/1‑100 (MQ=255)
         gTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTGTACCGGTACGGCAGTCAGCGCGatca                                                                                     <  1:2066445/100‑1 (MQ=255)
            agcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGAtcatca                                                                                  >  2:1013233/1‑100 (MQ=255)
             gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAc                                                                                 <  1:763611/100‑1 (MQ=255)
             gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAc                                                                                 >  1:1778729/1‑100 (MQ=255)
             gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAc                                                                                 >  2:3029254/1‑100 (MQ=255)
               agTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCg                                                                               <  2:884199/100‑1 (MQ=255)
                   aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAg                                                                           <  2:1612727/100‑1 (MQ=255)
                       gCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAg                                                                       <  2:2265497/100‑1 (MQ=255)
                           ccAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATgg                                                                   <  1:2081744/100‑1 (MQ=255)
                           ccAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATgg                                                                   >  2:1318615/1‑100 (MQ=255)
                            cAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGc                                                                  >  2:1657893/1‑100 (MQ=255)
                             aGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCg                                                                 <  1:1120968/100‑1 (MQ=255)
                              gggAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGt                                                                >  2:1727158/1‑100 (MQ=255)
                                 aaGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgca                                                             <  2:2459475/100‑1 (MQ=255)
                                  aGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcag                                                            <  2:2450325/100‑1 (MQ=255)
                                   gTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACTGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagc                                                           <  2:1947832/100‑1 (MQ=255)
                                    ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagca                                                          <  2:88416/100‑1 (MQ=255)
                                     tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                         >  2:2889692/1‑100 (MQ=255)
                                     tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTCACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                         >  1:1933395/1‑100 (MQ=255)
                                       aaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagca                                                          <  1:257650/97‑1 (MQ=255)
                                              tCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGCTc                                                >  1:1686223/1‑100 (MQ=255)
                                              tCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATc                                                >  2:1032166/1‑100 (MQ=255)
                                               cACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCg                                               >  1:1956101/1‑100 (MQ=255)
                                               cACCAGCCAGGCCAGCAGAATGCCGAACACCACGTTTACGGGCACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCg                                               <  1:2437157/100‑1 (MQ=255)
                                                aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                              >  2:1022544/1‑100 (MQ=255)
                                                aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                              <  1:499919/100‑1 (MQ=255)
                                                aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                              <  1:1213233/100‑1 (MQ=255)
                                                aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                              >  1:43668/1‑100 (MQ=255)
                                                aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGGTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                              >  2:297737/1‑100 (MQ=255)
                                                 ccagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGc                                             >  1:1571110/1‑100 (MQ=255)
                                                  cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGCTCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                            <  1:2946257/100‑1 (MQ=255)
                                                  cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                            <  2:2399683/100‑1 (MQ=255)
                                                  cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                            <  1:1164310/100‑1 (MQ=255)
                                                  cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                            >  1:3430127/1‑100 (MQ=255)
                                                  cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                            <  1:628931/100‑1 (MQ=255)
                                                      caggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGAt                                        <  1:2726759/100‑1 (MQ=255)
                                                          ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTg                                    <  1:1520685/100‑1 (MQ=255)
                                                          ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTg                                    >  2:3475476/1‑100 (MQ=255)
                                                           ccgcagAATGCCGGACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGt                                   >  2:2337605/3‑100 (MQ=255)
                                                           cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGt                                   <  1:610580/100‑1 (MQ=255)
                                                           cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGt                                   >  2:1396260/1‑100 (MQ=255)
                                                            gacgctcttccGATCTCCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTa                                  <  1:2805120/91‑1 (MQ=255)
                                                               agAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaaaa                               >  1:8546/1‑100 (MQ=255)
                                                                 aaTGCCGAACACCGGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAAcc                             <  2:3408939/100‑1 (MQ=255)
                                                                    gCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGc                          >  1:2394813/1‑100 (MQ=255)
                                                                    gCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGc                          >  1:1725032/1‑100 (MQ=255)
                                                                    gCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGc                          <  1:764237/100‑1 (MQ=255)
                                                                     ccGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                         >  1:241281/1‑100 (MQ=255)
                                                                       gAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATc                       <  1:135511/100‑1 (MQ=255)
                                                                       gAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATc                       >  2:608941/1‑100 (MQ=255)
                                                                        aaCACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCa                      <  1:1009241/100‑1 (MQ=255)
                                                                          cacCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGc                    >  2:1858348/1‑100 (MQ=255)
                                                                           acCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcc                   <  1:3202890/100‑1 (MQ=255)
                                                                            ccAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGccc                  >  2:2178558/1‑100 (MQ=255)
                                                                            ccAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGccc                  >  2:1208466/1‑100 (MQ=255)
                                                                            ccAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGccc                  >  1:906719/1‑100 (MQ=255)
                                                                             cAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcccc                 >  2:217460/1‑100 (MQ=255)
                                                                              aGGTTTAGCGGTACGGCAGTCAGCGCGATCATCACCGTCAGGCAGATGGCGTGCAGTATGTCAGGCTCGCCCAGCTTCTGTAAAACCGGCATCAGCCCCt                >  1:1710917/1‑100 (MQ=255)
                                                                               ggTTTATCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCtt               >  2:648753/1‑100 (MQ=255)
                                                                               ggTTTAGCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCACCATGGCCGGATCGGCCAGCTCCTGTAGAACCGGCATCAGcccc                 >  1:1025719/1‑98 (MQ=255)
                                                                               ggTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCtt               >  2:1962781/1‑100 (MQ=255)
                                                                                 tttACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGc             <  2:241284/100‑1 (MQ=255)
                                                                                 tttACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGc             >  1:159118/1‑100 (MQ=255)
                                                                                   tACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg           <  2:1549088/100‑1 (MQ=255)
                                                                                   tACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg           <  1:1148556/100‑1 (MQ=255)
                                                                                   tACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg           <  2:1496178/100‑1 (MQ=255)
                                                                                     ccGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGaa         <  1:3003688/100‑1 (MQ=255)
                                                                                      cGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAAt        <  1:517034/100‑1 (MQ=255)
                                                                                      cGGTACGCCAGTCAGCGCGATCATCACCGTCAGCCAGACGGCGTGCAGCATGTCCGGATCGGCCAGATCCTGTACAACCGGCATCAGCCCCTTGCTGAAt        >  1:3274438/1‑100 (MQ=255)
                                                                                       ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg       <  1:609624/100‑1 (MQ=255)
                                                                                       ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg       <  2:537214/100‑1 (MQ=255)
                                                                                       ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg       >  2:709635/1‑100 (MQ=255)
                                                                                       ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg       <  2:1139430/100‑1 (MQ=255)
                                                                                         tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc     <  2:720638/100‑1 (MQ=255)
                                                                                         tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc     <  1:1008724/100‑1 (MQ=255)
                                                                                         tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc     <  2:2949037/100‑1 (MQ=255)
                                                                                         tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc     <  1:1755388/100‑1 (MQ=255)
                                                                                          aCGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCt    >  2:1690260/1‑100 (MQ=255)
                                                                                          aCGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGCCCGGATCGGCCAGATGCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCt    <  1:820890/100‑1 (MQ=255)
                                                                                           cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg   <  2:3341924/100‑1 (MQ=255)
                                                                                           cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg   >  1:1835798/1‑100 (MQ=255)
                                                                                           cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg   <  1:211641/100‑1 (MQ=255)
                                                                                           cGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg   >  2:2201809/1‑100 (MQ=255)
                                                                                            ggCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGc  <  2:43668/100‑1 (MQ=255)
                                                                                            ggCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGc  <  1:2985692/100‑1 (MQ=255)
                                                                                                |                                                                                               
TCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGC  >  NC_000913/2541360‑2541551

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