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,9120TA98.5% 445.7 / ‑2.3 134I117I (ATT→ATAtehBtellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (75/57);  total (77/57)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

Only 100 of 134 total aligned reads displayed.

TGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTT  >  NC_000913/1501816‑1502008
                                                                                                |                                                                                                
tGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGcc                                                                                               >  1:4002666/1‑100 (MQ=255)
tGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGcc                                                                                               <  1:2244469/100‑1 (MQ=255)
    cTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAata                                                                                           <  1:685187/100‑1 (MQ=255)
    cTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAata                                                                                           >  2:4037457/1‑100 (MQ=255)
      gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                         >  1:625868/1‑100 (MQ=255)
          aaCCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAAc                                                                                     >  1:4009305/1‑100 (MQ=255)
           aCCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACg                                                                                    >  2:403444/1‑100 (MQ=255)
            ccTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGt                                                                                   >  1:2174293/1‑100 (MQ=255)
            ccTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGt                                                                                   >  1:2174308/1‑100 (MQ=255)
             cTCACATTTGATAGACAGTACGATTTTATTATTACGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAAGATGCAACGtt                                                                                  >  2:2738732/1‑100 (MQ=255)
               cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                                <  2:2199349/100‑1 (MQ=255)
               cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                                <  1:3052130/100‑1 (MQ=255)
               cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                                <  2:593874/100‑1 (MQ=255)
               cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                                <  2:3516000/100‑1 (MQ=255)
               cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                                >  2:3689427/1‑100 (MQ=255)
                 caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                              <  1:3129097/100‑1 (MQ=255)
                    ttGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTaa                                                                           <  1:1323522/100‑1 (MQ=255)
                     tGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTaaa                                                                          >  2:4011966/1‑100 (MQ=255)
                      gATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAc                                                                         >  2:2897130/1‑100 (MQ=255)
                      gATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAc                                                                         >  1:3029265/1‑100 (MQ=255)
                       aTAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCTCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAcc                                                                        >  2:1067570/1‑100 (MQ=255)
                        tAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCATTGAGGCTAAAACAATCCCCGGGTTGATAGCCAATATGCAACGTTGCCCTAAACCt                                                                       >  1:3671097/1‑100 (MQ=255)
                         aGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtg                                                                      <  2:4113948/100‑1 (MQ=255)
                          gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                     >  1:318332/1‑100 (MQ=255)
                           aCAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggt                                                                    <  2:3348844/100‑1 (MQ=255)
                            cAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtg                                                                   >  1:586521/1‑100 (MQ=255)
                            cAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtg                                                                   >  1:2179563/1‑100 (MQ=255)
                             aGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtgg                                                                  <  2:1653276/100‑1 (MQ=255)
                              gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTggtggt                                                                 <  1:517486/100‑1 (MQ=255)
                               tACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGtt                                                                <  1:2729970/100‑1 (MQ=255)
                               tACGATTTTACTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATACGCAACGTTGCACTAAACCTGGTGCtt                                                                >  2:3497740/1‑100 (MQ=255)
                                aCGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTa                                                               <  2:3122249/100‑1 (MQ=255)
                                 cGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTAc                                                              >  1:247826/1‑100 (MQ=255)
                                   aTTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTCCa                                                             >  2:2055783/1‑99 (MQ=255)
                                   aTTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACaa                                                            >  1:1952369/1‑100 (MQ=255)
                                   aTTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACaa                                                            <  2:625864/100‑1 (MQ=255)
                                   aTTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACaa                                                            >  2:3222644/1‑100 (MQ=255)
                                     tttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAcc                                                          <  1:2653058/100‑1 (MQ=255)
                                        aTTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGa                                                       >  2:3392028/1‑100 (MQ=255)
                                           cTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtg                                                    >  1:150808/1‑100 (MQ=255)
                                           cTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtg                                                    <  2:3550711/100‑1 (MQ=255)
                                             ttCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtgtg                                                  >  1:3521219/1‑100 (MQ=255)
                                               cGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggc                                                >  1:3756654/1‑100 (MQ=255)
                                                gACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcg                                               >  1:380946/1‑100 (MQ=255)
                                                gACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcg                                               <  2:1324276/100‑1 (MQ=255)
                                                  cTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcggc                                             >  1:435165/1‑100 (MQ=255)
                                                  cTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcggc                                             <  2:3972261/100‑1 (MQ=255)
                                                     tggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGAt                                          <  2:3600075/100‑1 (MQ=255)
                                                     tggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGAt                                          <  1:1673761/100‑1 (MQ=255)
                                                      ggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATg                                         <  1:4174016/100‑1 (MQ=255)
                                                      ggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATg                                         <  1:4174015/100‑1 (MQ=255)
                                                      ggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATg                                         <  1:3218985/100‑1 (MQ=255)
                                                      ggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATg                                         >  1:2837103/1‑100 (MQ=255)
                                                       gtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATgg                                        >  2:129474/1‑100 (MQ=255)
                                                       gtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATgg                                        <  2:1204448/100‑1 (MQ=255)
                                                       gtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATgg                                        >  2:2487044/1‑100 (MQ=255)
                                                       gtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATgg                                        <  2:913443/100‑1 (MQ=255)
                                                        tgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGa                                       <  1:1603125/100‑1 (MQ=255)
                                                           tGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAc                                    <  1:1270772/100‑1 (MQ=255)
                                                            gttGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGTCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAcc                                   <  2:1359026/98‑1 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 >  1:1065820/1‑100 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 >  2:3972052/1‑100 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 >  1:1006453/1‑100 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 >  1:35924/1‑100 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 <  2:1663489/100‑1 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACATGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 >  2:549205/1‑100 (MQ=255)
                                                                ttCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTg                               <  2:1679163/100‑1 (MQ=255)
                                                                  ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAt                             >  2:2694028/1‑100 (MQ=255)
                                                                   cTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAtt                            >  2:434527/1‑100 (MQ=255)
                                                                       aGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATcc                        <  2:3846078/100‑1 (MQ=255)
                                                                        ggCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCa                       <  1:2835477/100‑1 (MQ=255)
                                                                        ggCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCa                       >  2:371237/1‑100 (MQ=255)
                                                                        ggCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCa                       >  1:651343/1‑100 (MQ=255)
                                                                        ggCTAAAACCATCCCCGGGTTGATAGACAATATGCACCGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTCCTTCTCCa                       <  2:1649559/100‑1 (MQ=255)
                                                                         gCTAAAACCATCCCCGGGTTGCTAGCCAATATGCAACATTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGGTGATTAtgc                        >  2:3602846/1‑96 (MQ=255)
                                                                         gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                      >  2:1050563/1‑100 (MQ=255)
                                                                         gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                      >  1:1578496/1‑100 (MQ=255)
                                                                         gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                      <  1:3929567/100‑1 (MQ=255)
                                                                          cTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATg                     >  2:1690060/1‑100 (MQ=255)
                                                                          cTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATg                     >  1:738468/1‑100 (MQ=255)
                                                                           tAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGt                    <  2:2698202/100‑1 (MQ=255)
                                                                            aaaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTa                   >  1:361779/1‑100 (MQ=255)
                                                                             aaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAc                  <  1:549209/100‑1 (MQ=255)
                                                                               aCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCg                >  2:4103918/1‑100 (MQ=255)
                                                                               aCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCg                >  2:4007580/1‑100 (MQ=255)
                                                                                  aTCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCg             >  1:3789853/1‑100 (MQ=255)
                                                                                  aTCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCg             >  2:2973875/1‑100 (MQ=255)
                                                                                  aTCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCg             >  1:1096250/1‑100 (MQ=255)
                                                                                  aTCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCg             <  2:2556121/100‑1 (MQ=255)
                                                                                   tCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCgg            >  2:1816953/1‑100 (MQ=255)
                                                                                   tCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCgg            >  1:3227700/1‑100 (MQ=255)
                                                                                    ccccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGc           >  1:2906225/1‑100 (MQ=255)
                                                                                     cccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCt          >  1:3482972/1‑100 (MQ=255)
                                                                                     cccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCt          >  2:2282038/1‑100 (MQ=255)
                                                                                     cccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCt          >  2:774141/1‑100 (MQ=255)
                                                                                     cccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCt          >  1:929739/1‑100 (MQ=255)
                                                                                         ggTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTccc      <  2:246799/100‑1 (MQ=255)
                                                                                           ttGATAGACAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGt    >  1:2186112/1‑100 (MQ=255)
                                                                                             gATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGttt  <  2:2409408/100‑1 (MQ=255)
                                                                                             gATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGttt  >  1:3526235/1‑100 (MQ=255)
                                                                                                |                                                                                                
TGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTT  >  NC_000913/1501816‑1502008

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: