Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,541,456 A→G 60.4% 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,4560AG60.4% 101.0 / 135.0 131I75T (ATT→ACT) cysWsulfate/thiosulfate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (31/21);  new base G (45/34);  total (76/55)
Fisher's exact test for biased strand distribution p-value = 8.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 131 total aligned reads displayed.

GTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCAC  >  NC_000913/2541359‑2541553
                                                                                                 |                                                                                                 
gTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATc                                                                                                 >  1:200872/1‑100 (MQ=255)
 tCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCa                                                                                                >  2:1885329/1‑100 (MQ=255)
         cGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatc                                                                                        <  2:1649864/100‑1 (MQ=255)
          gTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatca                                                                                       <  2:2770642/100‑1 (MQ=255)
            cagcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatcatc                                                                                     >  1:3139086/1‑100 (MQ=255)
             agcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGAtcatca                                                                                    >  1:1272686/1‑100 (MQ=255)
              gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAc                                                                                   >  1:984196/1‑100 (MQ=255)
              gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAc                                                                                   >  2:2853076/1‑100 (MQ=255)
              gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAc                                                                                   >  2:1237936/1‑100 (MQ=255)
              gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAc                                                                                   <  1:2556353/100‑1 (MQ=255)
              gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAc                                                                                   >  2:1934057/1‑100 (MQ=255)
              gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAc                                                                                   >  2:1882624/1‑100 (MQ=255)
              gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCAc                                                                                   >  2:2142567/1‑100 (MQ=255)
              gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCAc                                                                                   >  2:188288/1‑100 (MQ=255)
               cagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAcc                                                                                  >  2:2193737/1‑100 (MQ=255)
               cagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAcc                                                                                  >  2:2243694/1‑100 (MQ=255)
                agTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCg                                                                                 >  1:795583/1‑100 (MQ=255)
                agTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCg                                                                                 >  2:2694441/1‑100 (MQ=255)
                 gTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATTCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGt                                                                                >  1:2875537/1‑100 (MQ=255)
                    tcTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAg                                                                             <  1:2112013/99‑1 (MQ=255)
                    aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAg                                                                             <  2:3102309/100‑1 (MQ=255)
                     cTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGc                                                                            >  1:516814/1‑100 (MQ=255)
                     cTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGc                                                                            >  1:601562/1‑100 (MQ=255)
                     cTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGc                                                                            >  2:772624/1‑100 (MQ=255)
                         cGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGa                                                                        <  2:86891/100‑1 (MQ=255)
                          gTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGa                                                                        <  2:462193/99‑1 (MQ=255)
                            ccAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATgg                                                                     >  1:1390501/1‑100 (MQ=255)
                            ccAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATgg                                                                     <  2:3139183/100‑1 (MQ=255)
                              aaggAAGTTAAAGCGCGTCACCAGCAAGGCCAGCAGAATGCCCAACAACAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCg                                                                   <  1:152161/98‑1 (MQ=255)
                              aGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCg                                                                   >  2:1559032/1‑100 (MQ=255)
                               gggAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGt                                                                  >  1:2588420/1‑100 (MQ=255)
                                 gAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgc                                                                <  2:74111/100‑1 (MQ=255)
                                  aaGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgca                                                               >  2:2680108/1‑100 (MQ=255)
                                   aGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcag                                                              >  1:1375233/1‑100 (MQ=255)
                                      tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagca                                                            <  1:2723521/99‑1 (MQ=255)
                                      tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                           <  2:2827647/100‑1 (MQ=255)
                                      tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                           >  2:1388757/1‑100 (MQ=255)
                                       aaaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg                                                          >  1:1599694/1‑100 (MQ=255)
                                        aaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGt                                                         >  1:25094/1‑100 (MQ=255)
                                           cgcgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCg                                                      <  1:1885317/100‑1 (MQ=255)
                                           cgcgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCg                                                      <  2:473160/100‑1 (MQ=255)
                                           cgcgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCg                                                      <  2:2153851/100‑1 (MQ=255)
                                             cgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGa                                                    >  2:350265/1‑100 (MQ=255)
                                               tCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCTGGATc                                                  >  1:2391694/1‑100 (MQ=255)
                                               tCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATc                                                  <  1:1660515/100‑1 (MQ=255)
                                                 aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGTATGTCCGGa                                                    >  1:1069500/1‑96 (MQ=255)
                                                    agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                                             >  1:725421/1‑100 (MQ=255)
                                                    agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                                             >  2:1896478/1‑100 (MQ=255)
                                                     gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAg                                            >  2:2876305/1‑100 (MQ=255)
                                                     gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAg                                            >  1:1929792/1‑100 (MQ=255)
                                                     gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAg                                            >  2:934593/1‑100 (MQ=255)
                                                         ggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTc                                        <  2:1225180/100‑1 (MQ=255)
                                                         ggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTc                                        >  2:2847174/1‑100 (MQ=255)
                                                         ggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTc                                        <  1:1543513/100‑1 (MQ=255)
                                                           ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTg                                      >  1:945981/1‑100 (MQ=255)
                                                           ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTg                                      <  2:1339030/100‑1 (MQ=255)
                                                           ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTg                                      <  1:2236887/100‑1 (MQ=255)
                                                           ccagCAGAAGGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCACCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGCTCGGCCAGATTCTg                                      <  1:671765/100‑1 (MQ=255)
                                                                 gAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAAc                                >  2:1217738/1‑100 (MQ=255)
                                                                 gAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAAc                                >  2:1454102/1‑100 (MQ=255)
                                                                   aTGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCg                              >  2:649815/1‑100 (MQ=255)
                                                                     gCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGc                            >  1:853249/1‑100 (MQ=255)
                                                                     gCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGc                            >  2:183519/1‑100 (MQ=255)
                                                                     gCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCAGCACCGTCAGCCAGATGTCGTTCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGc                            <  1:1679819/100‑1 (MQ=255)
                                                                      ccGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGAGTCTGTAAAACCGGCa                           >  2:3334017/1‑100 (MQ=255)
                                                                      ccGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                           <  2:2144599/100‑1 (MQ=255)
                                                                      ccGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCAACACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCACAGTCTGTAAAACCGGCa                           <  1:3171964/100‑1 (MQ=255)
                                                                      ccGAACACCAGGTCTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                           <  1:2410041/100‑1 (MQ=255)
                                                                       cGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAt                          <  1:1029062/100‑1 (MQ=255)
                                                                        gAACGCCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                           >  1:289197/1‑98 (MQ=255)
                                                                         aaCACCAGGTTTACCGGTACGGCAATCAGCGCGATCACCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATc                         <  1:932685/99‑1 (MQ=255)
                                                                          acacCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAg                       >  1:3142642/1‑100 (MQ=255)
                                                                          acacCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAg                       <  2:627842/100‑1 (MQ=255)
                                                                           cacCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGc                      <  2:929139/100‑1 (MQ=255)
                                                                            acCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcc                     >  2:627621/1‑100 (MQ=255)
                                                                               aGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCt                  <  1:680304/100‑1 (MQ=255)
                                                                                ggTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcccc                   >  1:1503510/1‑98 (MQ=255)
                                                                                   ttACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCt              <  1:2606164/100‑1 (MQ=255)
                                                                                    tACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg             <  2:2412300/100‑1 (MQ=255)
                                                                                    tACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg             >  2:2357402/1‑100 (MQ=255)
                                                                                     aCCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGa            <  1:384801/100‑1 (MQ=255)
                                                                                      ccGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGaa           <  2:2595012/100‑1 (MQ=255)
                                                                                      ccGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGaa           <  2:3021847/100‑1 (MQ=255)
                                                                                      ccGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGaa           >  2:333845/1‑100 (MQ=255)
                                                                                      ccGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGTCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGaa           >  2:1133773/1‑100 (MQ=255)
                                                                                      ccGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGaa           >  1:1358973/1‑100 (MQ=255)
                                                                                        ggTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg         >  1:100014/1‑100 (MQ=255)
                                                                                          tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc       >  2:813771/1‑100 (MQ=255)
                                                                                          tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc       <  2:15854/100‑1 (MQ=255)
                                                                                           aCGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCt      >  1:2316554/1‑100 (MQ=255)
                                                                                            cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc       <  1:2594822/98‑1 (MQ=255)
                                                                                            cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg     >  2:1013518/1‑100 (MQ=255)
                                                                                            cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg     >  2:1197129/1‑100 (MQ=255)
                                                                                            cGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg     >  1:1889069/1‑100 (MQ=255)
                                                                                             ggCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGc    >  2:1046284/1‑100 (MQ=255)
                                                                                             ggCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGc    >  1:2465827/1‑100 (MQ=255)
                                                                                             ggCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGc    <  1:758029/100‑1 (MQ=255)
                                                                                             ggCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGc    <  1:817840/100‑1 (MQ=255)
                                                                                               cAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCAc  <  1:26028/100‑1 (MQ=255)
                                                                                               cAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATCTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCAc  <  1:2172796/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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: