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,4560AG97.9% 315.6 / ‑2.7 92I75T (ATT→ACT) cysWsulfate/thiosulfate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (41/49);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCA  >  NC_000913/2541361‑2541552
                                                                                               |                                                                                                
ccAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAg                                                                                              >  2:1324416/1‑100 (MQ=255)
 cAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAgc                                                                                             >  1:1831481/1‑100 (MQ=255)
   gTAGCGTCGGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAgcgc                                                                                           <  1:107446/100‑1 (MQ=255)
   gTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAgcgc                                                                                           <  1:1890218/100‑1 (MQ=255)
    tAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGcgcg                                                                                          <  2:866942/100‑1 (MQ=255)
    tAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGcgcg                                                                                          <  2:1320721/100‑1 (MQ=255)
    tAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGcgcg                                                                                          <  1:2371864/100‑1 (MQ=255)
      gCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGat                                                                                        >  2:772936/1‑100 (MQ=255)
          cagcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCCAACACCGGGTTTACCGGTACGGCAGTCAGCGCGatcatc                                                                                    >  2:1795962/1‑100 (MQ=255)
            gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAc                                                                                  >  2:1547373/1‑100 (MQ=255)
            gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAc                                                                                  <  2:1871352/100‑1 (MQ=255)
             cagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAcc                                                                                 >  1:484986/1‑100 (MQ=255)
             cagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCAcc                                                                                 >  1:1837693/1‑100 (MQ=255)
                  aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAg                                                                            <  1:441957/100‑1 (MQ=255)
                   cTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGc                                                                           >  1:1647890/1‑100 (MQ=255)
                    tGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGcc                                                                          >  1:69685/1‑100 (MQ=255)
                     ggCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCa                                                                         <  1:533881/100‑1 (MQ=255)
                     ggCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCa                                                                         >  2:725507/1‑100 (MQ=255)
                       cGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGa                                                                       <  2:1622386/100‑1 (MQ=255)
                        gTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGAt                                                                      <  2:2501000/100‑1 (MQ=255)
                           cAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCGGTCAGCGCGATCATCACCGTCAGCCAGATGGc                                                                   <  1:2428144/100‑1 (MQ=255)
                           cAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGc                                                                   >  2:2012809/1‑100 (MQ=255)
                           cAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGc                                                                   >  2:1613224/1‑100 (MQ=255)
                            aGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCg                                                                  <  1:992777/100‑1 (MQ=255)
                             gggAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGt                                                                 >  1:931652/1‑100 (MQ=255)
                             gggAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGt                                                                 >  2:2060609/1‑100 (MQ=255)
                               gAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgc                                                               <  1:743626/100‑1 (MQ=255)
                               gAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgc                                                               <  2:1955415/100‑1 (MQ=255)
                                aaGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgca                                                              <  2:2209540/100‑1 (MQ=255)
                                   ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagca                                                           >  2:995176/1‑100 (MQ=255)
                                   ttAAAGCACGTCACCAGCCAGGCCAGCAGAAGGCGGAACGCCAGGTGTACCTGTACGGCAGTCAGCGCGATCATCACCGGCAGCTAGAGGGCGTgcagca                                                           <  2:1510817/100‑1 (MQ=255)
                                    tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                          <  2:1839205/100‑1 (MQ=255)
                                    tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                          >  2:1975260/1‑100 (MQ=255)
                                    tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                          >  1:1410806/1‑100 (MQ=255)
                                     aaaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg                                                         >  1:280284/1‑100 (MQ=255)
                                     aaaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGACAGATGGCGTGCAGCATg                                                         >  1:2158036/1‑100 (MQ=255)
                                     aaaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg                                                         >  1:2441314/1‑100 (MQ=255)
                                        gcgcGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTc                                                       <  2:2527014/99‑1 (MQ=255)
                                          gcgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCgg                                                    <  2:1133195/100‑1 (MQ=255)
                                            gTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGAt                                                  >  1:1506004/1‑100 (MQ=255)
                                            gTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGAt                                                  >  1:1851166/1‑100 (MQ=255)
                                              cACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCg                                                <  1:427790/100‑1 (MQ=255)
                                               aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                               <  1:1044922/100‑1 (MQ=255)
                                               aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                               >  1:1115932/1‑100 (MQ=255)
                                               aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                               <  2:2518450/100‑1 (MQ=255)
                                               aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                               <  2:2496690/100‑1 (MQ=255)
                                               aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                               >  2:1742170/1‑100 (MQ=255)
                                               aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                               >  2:1102016/1‑100 (MQ=255)
                                                ccagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGc                                              >  1:1418154/1‑100 (MQ=255)
                                                 cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                             <  2:1287444/100‑1 (MQ=255)
                                                 cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                             >  2:1036987/1‑100 (MQ=255)
                                                  agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                                            <  1:1613235/100‑1 (MQ=255)
                                                  agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                                            <  1:2296875/100‑1 (MQ=255)
                                                  agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                                            <  1:2254401/100‑1 (MQ=255)
                                                  agccagGCCAGCAGAATGCCGAACACCAGATTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                                            <  1:2209379/100‑1 (MQ=255)
                                                   gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAg                                           <  1:1714401/100‑1 (MQ=255)
                                                   gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAg                                           <  2:969381/100‑1 (MQ=255)
                                                     caggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGAt                                         <  2:764102/100‑1 (MQ=255)
                                                     caggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGAt                                         >  2:264360/1‑100 (MQ=255)
                                                      aggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGAtt                                        <  1:1036995/100‑1 (MQ=255)
                                                       ggccagCAGAATGCCGAACACCAGGTGTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTc                                       >  2:1460254/1‑100 (MQ=255)
                                                       gatctgcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTc                                       <  2:1119433/95‑1 (MQ=255)
                                                         ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTg                                     <  2:2386590/100‑1 (MQ=255)
                                                         ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTg                                     <  2:1213825/100‑1 (MQ=255)
                                                          cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGt                                    <  2:1831467/100‑1 (MQ=255)
                                                           agcagAATGCCGAACCCCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTa                                   <  2:2233558/100‑1 (MQ=255)
                                                            gcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTcagatc                                  >  1:1119441/1‑95 (MQ=255)
                                                            gcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaa                                  >  1:2366079/1‑100 (MQ=255)
                                                             cagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCACCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCCGCCAGATTCTGTaaa                                 <  2:2469251/100‑1 (MQ=255)
                                                              agAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaaaa                                >  1:1196253/1‑100 (MQ=255)
                                                                aaTGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAAcc                              <  2:1411113/100‑1 (MQ=255)
                                                                 aTGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCg                             >  2:1073135/1‑100 (MQ=255)
                                                                   gCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGc                           >  2:829513/1‑100 (MQ=255)
                                                                    ccGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                          >  1:1659128/1‑100 (MQ=255)
                                                                     cGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAt                         <  1:877157/100‑1 (MQ=255)
                                                                     cGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAt                         >  1:1837418/1‑100 (MQ=255)
                                                                     cGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAt                         >  2:1175794/1‑100 (MQ=255)
                                                                     cGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAt                         >  1:2226644/1‑100 (MQ=255)
                                                                      gAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATc                        <  1:2311707/100‑1 (MQ=255)
                                                                      gAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATc                        >  1:184294/1‑100 (MQ=255)
                                                                       aaCACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCa                       <  2:1785888/100‑1 (MQ=255)
                                                                           ccAGTTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGAGCGGCCAGATTCTGTAAAACCGGCATCAGccc                   <  2:488311/100‑1 (MQ=255)
                                                                            cAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcccc                  <  2:1584280/100‑1 (MQ=255)
                                                                              ggTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCtt                >  2:334806/1‑100 (MQ=255)
                                                                               gTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTg               >  1:1254703/1‑100 (MQ=255)
                                                                                tttACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGTCGTGCCGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGc              <  2:3171/100‑1 (MQ=255)
                                                                                  tACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg            <  2:2262045/100‑1 (MQ=255)
                                                                                  tACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg            <  1:1102024/100‑1 (MQ=255)
                                                                                      ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg        >  2:1765576/1‑100 (MQ=255)
                                                                                      ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg        >  2:2116730/1‑100 (MQ=255)
                                                                                          cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg    <  1:1653092/100‑1 (MQ=255)
                                                                                            gCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCa  <  2:581749/100‑1 (MQ=255)
                                                                                               |                                                                                                
CCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCA  >  NC_000913/2541361‑2541552

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: