Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,541,456 A→G 100% 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.1% 365.0 / ‑4.8 108I75T (ATT→ACT) cysWsulfate/thiosulfate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (51/56);  total (51/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

Only 100 of 109 total aligned reads displayed.

GTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCAC  >  NC_000913/2541359‑2541553
                                                                                                 |                                                                                                 
gTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTc                                                                                                 <  2:242022/100‑1 (MQ=255)
 tCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTc                                                                                                 <  1:2400838/99‑1 (MQ=255)
 tCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCa                                                                                                >  2:2403929/1‑100 (MQ=255)
 tCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCa                                                                                                >  1:1646993/1‑100 (MQ=255)
  ccAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCAGAACACCAGGTTTACCGGTACGGCAGTCa                                                                                                >  1:1983758/1‑99 (MQ=255)
    aGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAgcg                                                                                             <  1:2527857/100‑1 (MQ=255)
       aGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGa                                                                                          >  2:1288486/1‑100 (MQ=255)
        gCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGat                                                                                         <  1:1199467/100‑1 (MQ=255)
         cGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGTCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatc                                                                                        <  1:2310900/100‑1 (MQ=255)
         cGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatc                                                                                        >  2:2814785/1‑100 (MQ=255)
          gTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACAAGATTTACCAATACGGCAGTCAGCGCGataa                                                                                       >  1:1893547/1‑98 (MQ=255)
           tCAGCAGTAACTGGCGTCCAGGGAAGTTAAATCGCTTCACCAGCTAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatca                                                                                       <  1:2335366/99‑1 (MQ=255)
           ctcttccgatcTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatcat                                                                                      <  1:2203903/90‑1 (MQ=255)
             agcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGAtcatca                                                                                    >  1:2626283/1‑100 (MQ=255)
               cagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCGGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAcc                                                                                  >  2:1610040/1‑100 (MQ=255)
                 gTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGt                                                                                <  1:855522/100‑1 (MQ=255)
                 gTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGt                                                                                >  2:1783240/1‑100 (MQ=255)
                 gTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGt                                                                                >  2:1891187/1‑100 (MQ=255)
                     cTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTAAGCGCGATCATCACCGTCAGc                                                                            >  2:1041867/1‑100 (MQ=255)
                      tGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGcc                                                                           <  1:1876673/100‑1 (MQ=255)
                        gCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCACCAGAATGCCGAACACCAGGATTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAg                                                                         <  2:1646974/100‑1 (MQ=255)
                         cGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGa                                                                        >  2:1911698/1‑100 (MQ=255)
                            ccAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATgg                                                                     <  1:1891209/100‑1 (MQ=255)
                             cAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGc                                                                    >  1:3045293/1‑100 (MQ=255)
                             cAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGc                                                                    >  2:2996249/1‑100 (MQ=255)
                             cAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGc                                                                    >  1:1217061/1‑100 (MQ=255)
                              aGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATgg                                                                     <  2:137342/98‑1 (MQ=255)
                              aGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCg                                                                   <  2:254244/100‑1 (MQ=255)
                               gggAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGt                                                                  >  1:190108/1‑100 (MQ=255)
                                  aaGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgca                                                               <  1:1901344/100‑1 (MQ=255)
                                   aGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcag                                                              >  1:2474928/1‑100 (MQ=255)
                                    gTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagc                                                             >  2:278740/1‑100 (MQ=255)
                                    gTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagc                                                             >  2:88651/1‑100 (MQ=255)
                                     ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagca                                                            <  2:647609/100‑1 (MQ=255)
                                     ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagca                                                            >  1:1770612/1‑100 (MQ=255)
                                     ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagca                                                            >  1:2904083/1‑100 (MQ=255)
                                      tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGATACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                           >  1:947474/1‑100 (MQ=255)
                                        aaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGt                                                         <  1:2570509/100‑1 (MQ=255)
                                        aaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGATACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGt                                                         >  1:761292/1‑100 (MQ=255)
                                         aGCGCGTCACTATCTAGTCTAGCAGAATGTTGATCTCCAGGTTTACCGGTACGGCAGTCCGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTc                                                        <  1:1646485/100‑1 (MQ=255)
                                           cgcgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCg                                                      <  2:311330/100‑1 (MQ=255)
                                            gcgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCAACACCATCAGCCAGATGGCGTGCAGCATGTcc                                                       >  1:1657167/1‑98 (MQ=255)
                                               tCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATc                                                  >  2:1246823/1‑100 (MQ=255)
                                                  ccagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGc                                               <  2:2696280/100‑1 (MQ=255)
                                                  ccagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTACGGATCGGc                                               >  2:731442/1‑100 (MQ=255)
                                                   cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGACAGCCAGATGGCGTGCAGAATGTCCGGATCGGcc                                              >  1:248572/1‑100 (MQ=255)
                                                    agccagGCCAGCAGAATTCCGAACACCAGTTTTACCGGTACGTCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                                <  2:1311256/97‑1 (MQ=255)
                                                    agccagGCCAGCAGAATGCCGAACACCAGGTTTATCGGTACGGTAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                                             <  1:2331435/100‑1 (MQ=255)
                                                     gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAg                                            <  1:530929/100‑1 (MQ=255)
                                                         ggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTc                                        <  2:2296776/100‑1 (MQ=255)
                                                         ggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCAACACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTc                                        >  1:900290/1‑100 (MQ=255)
                                                          gccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCt                                       >  1:525005/1‑100 (MQ=255)
                                                            cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCTATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGt                                     <  1:38794/100‑1 (MQ=255)
                                                             agcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGTGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTa                                    <  1:653930/100‑1 (MQ=255)
                                                             agcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTa                                    <  1:3013319/100‑1 (MQ=255)
                                                             agcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTa                                    <  2:2763643/100‑1 (MQ=255)
                                                              gcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCATCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTa                                    <  2:1657147/99‑1 (MQ=255)
                                                              gcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaa                                   >  2:3010872/1‑100 (MQ=255)
                                                                   aTGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCG                              >  1:357611‑M1/1‑100 (MQ=255)
                                                                     gCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGC                            <  1:1911721‑M1/100‑1 (MQ=255)
                                                                      ccGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCA                           <  1:1099292‑M1/100‑1 (MQ=255)
                                                                      ccGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCA                           <  1:2723218‑M1/100‑1 (MQ=255)
                                                                      ccGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCA                           <  1:1246835‑M1/100‑1 (MQ=255)
                                                                          acacCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAg                       <  2:718528/100‑1 (MQ=255)
                                                                           cacCAGGCTCACCCGTACGGCAGACAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATATCCGGATCGGCCAGAGTCTGTAAAACCGGCATCAGc                      >  2:114245/1‑100 (MQ=255)
                                                                            acCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCTTCATCCAGATGGCATACAGCATATCATAAACGGCCAGATTCTGTAAAACCAGCa                           >  1:9342/1‑94 (MQ=255)
                                                                            acCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcc                     >  1:1070866/1‑100 (MQ=255)
                                                                             ccAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGccc                    >  1:617912/1‑100 (MQ=255)
                                                                             ccAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTACGGATCGGCCAGATTCTGTAAAAACGGCATCAGcc                     >  1:2837128/1‑99 (MQ=255)
                                                                              cAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCTAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcccc                   >  2:1018228/1‑100 (MQ=255)
                                                                               aGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGccccn                  <  2:1737403/100‑2 (MQ=255)
                                                                               aGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCt                  <  1:2245443/100‑1 (MQ=255)
                                                                                  tttACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAGCAGCCCCTTGc               <  2:949080/100‑1 (MQ=255)
                                                                                   ttACCGGTACGGCAGTCAGCACGATCATCACCGTCAGCCAGATGGCGTGCAGCATCTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCt              >  2:591553/1‑100 (MQ=255)
                                                                                    tACCGGTACGGCAGTCAGCGCGATCATCATCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCAAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg             <  2:1826593/100‑1 (MQ=255)
                                                                                    tACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg             <  2:1562341/100‑1 (MQ=255)
                                                                                     aCCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGa            >  1:2029019/1‑100 (MQ=255)
                                                                                     aCCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCAGATCAGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGa            >  1:882744/1‑100 (MQ=255)
                                                                                      ccGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAa           <  1:2562573‑M2/100‑2 (MQ=255)
                                                                                        ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAatg         >  1:1041874‑M2/1‑97 (MQ=255)
                                                                                          tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATtcagcaaggggctgatgcc       <  1:2314459‑M1/100‑20 (MQ=255)
                                                                                          tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATtcagcaaggggctgatgcc       <  2:2631022‑M1/100‑20 (MQ=255)
                                                                                          tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAatgcc       >  2:1343087‑M2/1‑95 (MQ=255)
                                                                                          tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAatgcc       >  1:721638‑M2/1‑95 (MQ=255)
                                                                                          tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAatgc        >  1:1973717‑M2/1‑95 (MQ=255)
                                                                                          tACGGCAGTCAGAGCGATCATTACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCAAGATTCTGTACAACCGGTATtcagcaaggggctgatgcc       <  2:133822‑M1/100‑20 (MQ=255)
                                                                                          tACGGCAGGCAGCGCGATCACCACCGACAGCCAGATGGCGTGCAGCATGTCCGTATCGGCCAGATTCTGTAAAACCGGCATtcagcaaggggctgatgcc       <  2:1159528‑M1/100‑20 (MQ=255)
                                                                                           aCGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATtcagcaaggggctgatgccg      <  1:1606168‑M1/100‑21 (MQ=255)
                                                                                           aCGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCt      <  2:1770592/100‑1 (MQ=255)
                                                                                           aCGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCt      <  2:2604117/100‑1 (MQ=255)
                                                                                            cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGTGTTCAGCATTTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg     <  1:783504/100‑1 (MQ=255)
                                                                                            cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAatagatc     <  2:658372‑M2/100‑8 (MQ=255)
                                                                                             ggCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGc    <  2:112751/100‑1 (MQ=255)
                                                                                             ggCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGc    <  2:1204983/100‑1 (MQ=255)
                                                                                             ggCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGc    <  2:183796/100‑1 (MQ=255)
                                                                                             ggCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGc    <  2:901967/100‑1 (MQ=255)
                                                                                              gCAGTCAGCGCGATCTTCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCa   <  1:1047039/100‑1 (MQ=255)
                                                                                              gCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCa   >  2:1257052/1‑100 (MQ=255)
                                                                                              gCAGTCAGCGCGATCATCACCGTCAGCAAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATtcagcaaggggctgatgccggt    <  1:723930‑M1/99‑23 (MQ=255)
                                                                                               cAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCAc  <  2:2541509/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
GTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCAC  >  NC_000913/2541359‑2541553

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