Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,501,912 T→A 100% I117I (ATT→ATA)  tehB → tellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,501,9120TA99.3% 458.1 / ‑5.2 133I117I (ATT→ATA) tehBtellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (70/62);  total (71/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

TGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTT  >  NC_000913/1501816‑1502008
                                                                                                |                                                                                                
tGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGcc                                                                                               >  2:561397/1‑100 (MQ=255)
 gATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCa                                                                                              >  2:951962/1‑100 (MQ=255)
 gATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCa                                                                                              >  1:3167020/1‑100 (MQ=255)
   tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAat                                                                                            <  2:724658/100‑1 (MQ=255)
   tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAat                                                                                            >  1:2522202/1‑100 (MQ=255)
   tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAat                                                                                            <  2:3426909/100‑1 (MQ=255)
    cTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAata                                                                                           >  1:1551447/1‑100 (MQ=255)
    cTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAata                                                                                           >  2:500410/1‑100 (MQ=255)
    cTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAata                                                                                           >  2:1630646/1‑100 (MQ=255)
     tGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAatat                                                                                          <  1:150561/100‑1 (MQ=255)
      gAATAACCTCACATTTGATAGACAGTACGATTTTGTTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAat                                                                                            <  1:3463968/97‑1 (MQ=255)
      gAATAACATCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGAGCTGATGTTCCTTGAGGCTAAAACCATCTCCGGGTTGATAGCCAATATg                                                                                         >  2:3699455/1‑100 (MQ=255)
        aTAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCa                                                                                       <  2:240675/100‑1 (MQ=255)
          aaCCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAAc                                                                                     <  2:2971788/100‑1 (MQ=255)
          aaCCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAAc                                                                                     <  1:1325672/100‑1 (MQ=255)
            ccTCACATTTGATAGACAGTACGATTTTATTCTTTCTACTGTGGTGATGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCaa                                                                                      >  1:2791229/1‑97 (MQ=255)
            ccTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGt                                                                                   <  2:2755199/100‑1 (MQ=255)
            ccTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGt                                                                                   >  2:3347508/1‑100 (MQ=255)
             cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGtt                                                                                  >  2:2584499/1‑100 (MQ=255)
             cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGCTTGATAGCCAATATGCAAc                                                                                     >  1:3152525/1‑97 (MQ=255)
              tCACATTTGATAGAGAGTACGATTTTATTCGTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTg                                                                                 >  1:1422462/1‑100 (MQ=255)
              tCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTg                                                                                 <  2:1465512/100‑1 (MQ=255)
               cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                                <  1:2529745/100‑1 (MQ=255)
                 taTTTGATAGACCGAACGCTTTTATTCTTTCGACTGTGGTGCTGGTGTTGCTTGGCGCTAAAACCCTGCCCGGGTCGATAGCCAATATGGAACGtgg                                                                                 >  1:3005738/2‑95 (MQ=255)
                 caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                                <  1:2781585/98‑1 (MQ=255)
                 caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                              <  2:727647/100‑1 (MQ=255)
                 caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                              >  2:68117/1‑100 (MQ=255)
                  aTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACt                                                                             >  2:3449420/1‑100 (MQ=255)
                     tGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTaaa                                                                          <  2:2790210/100‑1 (MQ=255)
                      gATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAc                                                                         <  2:1343233/100‑1 (MQ=255)
                          gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                     >  1:3832806/1‑100 (MQ=255)
                            cAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtg                                                                   >  2:1144522/1‑100 (MQ=255)
                             aGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtgg                                                                  <  2:4125411/100‑1 (MQ=255)
                             aGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtgg                                                                  <  1:1790623/100‑1 (MQ=255)
                              gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTggtggt                                                                 >  1:2655485/1‑100 (MQ=255)
                              gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTggtggt                                                                 <  1:1144519/100‑1 (MQ=255)
                              gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTggtggt                                                                 >  2:1643438/1‑100 (MQ=255)
                                aCGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTa                                                               >  1:3220042/1‑100 (MQ=255)
                                 cGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTAc                                                              >  1:221518/1‑100 (MQ=255)
                                   aTTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACaa                                                            >  1:905052/1‑100 (MQ=255)
                                     tttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAcc                                                          <  1:3598294/100‑1 (MQ=255)
                                     tttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAcc                                                          >  2:395514/1‑100 (MQ=255)
                                     tttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAcc                                                          >  2:396014/1‑100 (MQ=255)
                                     tttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAcc                                                          >  2:2749689/1‑100 (MQ=255)
                                     tttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAcc                                                          <  1:981520/100‑1 (MQ=255)
                                     tttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAc                                                           <  1:3667177/99‑1 (MQ=255)
                                        aTTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGa                                                       <  1:68117/100‑1 (MQ=255)
                                        aTTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGa                                                       <  1:1488267/100‑1 (MQ=255)
                                          tCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGAtt                                                     >  2:2749903/1‑100 (MQ=255)
                                           cTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtg                                                    >  2:1792717/1‑100 (MQ=255)
                                           cTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtg                                                    >  1:2726889/1‑100 (MQ=255)
                                            tttCGACTGTGGTGCTGATGTTCCTTGAGGCTAAACCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTACACCTGGTGGTTACAACCTGATtgt                                                   >  1:1284566/1‑100 (MQ=255)
                                             ttCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtgtg                                                  >  1:1079852/1‑100 (MQ=255)
                                               cGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggc                                                <  2:2010761/100‑1 (MQ=255)
                                                  cTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcggc                                             >  2:968635/1‑100 (MQ=255)
                                                   tgtgGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGgcggcg                                            <  2:366839/100‑1 (MQ=255)
                                                     tggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGAt                                          <  2:1431993/100‑1 (MQ=255)
                                                      ggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGa                                           >  1:711719/1‑98 (MQ=255)
                                                       gtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATgg                                        <  1:1792713/100‑1 (MQ=255)
                                                         gCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGAt                                      <  2:2083407/100‑1 (MQ=255)
                                                          cTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATa                                     >  1:3606145/1‑100 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTACACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATCCCGc                                 >  2:1350040/1‑100 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 <  2:2520450/100‑1 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 <  1:817649/100‑1 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 <  2:3139613/100‑1 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 >  1:2230590/1‑100 (MQ=255)
                                                                ttCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTg                               <  1:500408/100‑1 (MQ=255)
                                                                ttCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTg                               <  1:1630642/100‑1 (MQ=255)
                                                                 tCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTg                               >  1:3905419/1‑99 (MQ=255)
                                                                   cTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAtt                            >  2:1808059/1‑100 (MQ=255)
                                                                    ttGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAttn                           <  1:1807412/100‑2 (MQ=255)
                                                                    ttGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTa                           <  1:4188845/100‑1 (MQ=255)
                                                                     tGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTAt                          >  2:371968/1‑100 (MQ=255)
                                                                       acgCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATcc                        <  2:2205813/98‑1 (MQ=255)
                                                                       aGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATcc                        <  2:3723549/100‑1 (MQ=255)
                                                                        ggCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATAGATACCGCTGATTATCCa                       >  1:3348257/1‑100 (MQ=255)
                                                                        ggCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAACTGATTGTGGCGGCGCTGGATACCGCTGATTATCCa                       >  1:1059809/1‑100 (MQ=255)
                                                                        ggCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAATAAACCTGGTGGTTACAACCTGATTGTGGCGGCGCTGGATACCGCTGATTATc                         >  1:780245/1‑98 (MQ=255)
                                                                         gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                      <  1:3408742/100‑1 (MQ=255)
                                                                           tAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGt                    >  1:867995/1‑100 (MQ=255)
                                                                           tAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGt                    >  1:3814495/1‑100 (MQ=255)
                                                                             aaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAc                  >  2:3287037/1‑100 (MQ=255)
                                                                             aaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAc                  >  1:3953441/1‑100 (MQ=255)
                                                                             aaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAc                  <  2:4252501/100‑1 (MQ=255)
                                                                             aaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACATGGTGGTTACAACCTGATTGTGGCGGCGCTGGATACCGCTGATTATCCATGTAc                  >  2:560928/1‑100 (MQ=255)
                                                                              aaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAcc                 >  1:3970406/1‑100 (MQ=255)
                                                                                ccATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGt               >  2:2561850/1‑100 (MQ=255)
                                                                                 cATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTc              <  2:92296/100‑1 (MQ=255)
                                                                                 cATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTc              >  1:4039226/1‑100 (MQ=255)
                                                                                    ccccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGc           >  2:3740681/1‑100 (MQ=255)
                                                                                    ccccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGc           <  1:3130405/100‑1 (MQ=255)
                                                                                    ccccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGc           <  1:429640/100‑1 (MQ=255)
                                                                                       cGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTc        <  1:3230950/100‑1 (MQ=255)
                                                                                        gggTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTcc       <  2:1696574/100‑1 (MQ=255)
                                                                                        gggTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTcc       <  1:2386194/100‑1 (MQ=255)
                                                                                         ggTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTccc      <  2:3465885/100‑1 (MQ=255)
                                                                                           ttGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGt    <  2:3362800/100‑1 (MQ=255)
                                                                                             gATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGttt  >  2:2626580/1‑100 (MQ=255)
                                                                                             gATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGttt  >  2:3780416/1‑100 (MQ=255)
                                                                                             gATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACAGTCGGCTTCCCGttt  >  2:4208045/1‑100 (MQ=255)
                                                                                                |                                                                                                
TGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTT  >  NC_000913/1501816‑1502008

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: ‑