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,9120TA100.0% 383.5 / ‑6.7 115I117I (ATT→ATAtehBtellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (52/62);  total (52/63)

Only 100 of 115 total aligned reads displayed.

TGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGT  >  NC_000913/1501816‑1502006
                                                                                                |                                                                                              
tGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGcc                                                                                             >  1:314561/1‑100 (MQ=255)
 gatattaaTAACCTCACATTTGATAGGCAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCa                                                                                            <  2:865629/94‑1 (MQ=255)
 gATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCa                                                                                            >  1:509140/1‑100 (MQ=255)
 gATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCa                                                                                            >  2:29727/1‑100 (MQ=255)
 gATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCa                                                                                            >  1:1776660/1‑100 (MQ=255)
   tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAat                                                                                          >  2:679891/1‑100 (MQ=255)
    cTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAata                                                                                         >  1:4087191/1‑100 (MQ=255)
     tGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAatat                                                                                        >  2:2283489/1‑100 (MQ=255)
     tGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAata                                                                                         <  2:1039369/99‑1 (MQ=255)
     tGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTCTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAACCCATCCCCGGGTTGATAGCCActat                                                                                        >  2:2308968/1‑100 (MQ=255)
      gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGGCTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                       <  2:2639336/100‑1 (MQ=255)
      gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTACAACCATCCCCGGGTTGATAGCCAATATg                                                                                       <  1:1959370/100‑1 (MQ=255)
      gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                       >  1:1971810/1‑100 (MQ=255)
      gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                       <  1:2067726/100‑1 (MQ=255)
       aaTAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGc                                                                                      <  1:3869259/100‑1 (MQ=255)
       aaTAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGc                                                                                      <  1:1956758/100‑1 (MQ=255)
           aCCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACg                                                                                  >  2:3462204/1‑100 (MQ=255)
            ccTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGt                                                                                 >  2:2921491/1‑100 (MQ=255)
             cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGtt                                                                                >  1:1368484/1‑100 (MQ=255)
                 caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                            <  1:3227897/100‑1 (MQ=255)
                 caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                            <  1:2647135/100‑1 (MQ=255)
                 caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                            <  1:2268903/100‑1 (MQ=255)
                  aTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACt                                                                           >  1:2594457/1‑100 (MQ=255)
                    ttGATAGACCGTACGATTTTATTCTTTCGACTGTGGGGCTGATGTTCTTTGAGGCTAAAACAGTCCCCGGGTTGATAGCTAATATGCAACGTTACACt                                                                           >  2:3498080/1‑98 (MQ=255)
                    ttGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTaa                                                                         <  2:1758664/100‑1 (MQ=255)
                     tGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTaaa                                                                        >  1:2290408/1‑100 (MQ=255)
                      gATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAc                                                                       <  2:1675864/100‑1 (MQ=255)
                        tAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAcc                                                                      >  1:1906187/1‑99 (MQ=255)
                          gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                   <  1:1043887/100‑1 (MQ=255)
                          gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                   <  2:4087164/100‑1 (MQ=255)
                          gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                   <  1:2226487/100‑1 (MQ=255)
                          gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCAGGTTGTTAGCCAATATGCAACGTTGCACTAAACCtgg                                                                   >  2:2784192/1‑100 (MQ=255)
                           aCAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggt                                                                  <  1:2593964/100‑1 (MQ=255)
                            cAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtg                                                                 <  2:1099095/100‑1 (MQ=255)
                             aGTACGACTTTATGCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCGAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtgg                                                                <  1:1013476/100‑1 (MQ=255)
                              gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGCGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTggtggt                                                               >  1:1356778/1‑100 (MQ=255)
                              gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggt                                                                  >  2:2820195/1‑97 (MQ=255)
                              gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTggtggt                                                               >  2:3701777/1‑100 (MQ=255)
                               tACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGtt                                                              <  2:1718216/100‑1 (MQ=255)
                               tACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGtt                                                              >  1:4042980/1‑100 (MQ=255)
                                aCGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTa                                                             <  2:2352948/100‑1 (MQ=255)
                                   aTTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACaa                                                          <  2:2996726/100‑1 (MQ=255)
                                    ttttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAc                                                         <  1:3984234/100‑1 (MQ=255)
                                    ttttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAc                                                         <  1:2596766/100‑1 (MQ=255)
                                    ttttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAc                                                         <  1:3053574/100‑1 (MQ=255)
                                     tttATTCTTTCGACTGTGGTGCTGATGTTCCTTTAGTCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTATAAcc                                                        >  2:1304840/1‑100 (MQ=255)
                                     tttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCGATATGCAACGTTGCACTAAACCTGGTGGTTACAAcc                                                        <  1:2921507/100‑1 (MQ=255)
                                      ttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCt                                                       <  2:2586480/100‑1 (MQ=255)
                                      ttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCt                                                       >  2:4014399/1‑100 (MQ=255)
                                       tATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTg                                                      >  2:3391274/1‑100 (MQ=255)
                                         ttCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGAt                                                    <  1:3990739/100‑1 (MQ=255)
                                         ttCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCCAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGAt                                                    <  2:2174554/100‑1 (MQ=255)
                                           cTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtg                                                  <  2:1132568/100‑1 (MQ=255)
                                           cTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtg                                                  >  2:453564/1‑100 (MQ=255)
                                           cTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtg                                                  <  1:3707891/100‑1 (MQ=255)
                                            tttCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACATGGTGGTTACAACCTGATtgt                                                 >  2:806203/1‑100 (MQ=255)
                                              tCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTgg                                               <  1:1991994/100‑1 (MQ=255)
                                              tCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTCGCACTAAACCTGGTGGTTACAACCTGATTGTgg                                               <  2:2694125/100‑1 (MQ=255)
                                               cGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggc                                              <  2:1524456/100‑1 (MQ=255)
                                                gACTTTGGTGCTGATGTTCCTTGAGGCTAAAACCATCTCCGGGTTGATAGCAAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcg                                             <  2:1109436/100‑1 (MQ=255)
                                                gACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcg                                             >  1:3596859/1‑100 (MQ=255)
                                                 aCTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcgg                                            <  1:1205497/100‑1 (MQ=255)
                                                  cTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcggc                                           >  2:3966546/1‑100 (MQ=255)
                                                    gttgtgCTGATGTTCCTGGAGGCTAAAACCATCCCCTGGTTGATTGCCAATATGCAACGTGGCACTAAACCTGGTGGTTACAACCTGATTGTGgcggcg                                          <  2:1799156/96‑1 (MQ=255)
                                                    gtggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGa                                         >  2:3871959/1‑100 (MQ=255)
                                                    gtggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGa                                         >  1:2957482/1‑100 (MQ=255)
                                                      ggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGAt                                        <  2:2900132/99‑1 (MQ=255)
                                                      ggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATg                                       <  1:2283502/100‑1 (MQ=255)
                                                      ggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATg                                       >  1:2311113/1‑100 (MQ=255)
                                                      ggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATg                                       <  2:3135002/100‑1 (MQ=255)
                                                      ggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATg                                       <  1:1107205/100‑1 (MQ=255)
                                                         gCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGAt                                    <  1:2576064/100‑1 (MQ=255)
                                                           tGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAc                                  <  1:1514535/100‑1 (MQ=255)
                                                           tGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCAAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAc                                  <  2:3221981/100‑1 (MQ=255)
                                                            gATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAACTGGTGGTTACAAACTGATTGTGGCGGCGATGGATa                                   >  2:3339461/1‑98 (MQ=255)
                                                            cttGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAcc                                 <  2:2483498/98‑1 (MQ=255)
                                                             aTGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCg                                >  1:599969/1‑100 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAccag                               >  1:2483512/1‑98 (MQ=255)
                                                              tGTACCTTGAGTCTAAAACCATCCCCGGGTTGATAGACAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                               <  2:1838788/100‑1 (MQ=255)
                                                               gTTCCTTGAGGCTACAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTTGTGGTTACAACCTGATTATGGCGGCGATGGATAccaat                              >  2:1290950/1‑97 (MQ=255)
                                                               gTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCt                              >  2:167338/1‑100 (MQ=255)
                                                                ttCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTg                             <  1:3706294/100‑1 (MQ=255)
                                                                  ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGTCGGCGATGGATACCGCTGAt                           <  2:2519837/100‑1 (MQ=255)
                                                                  ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGCCGGCGATGGATACCGCTGAt                           <  2:2515736/100‑1 (MQ=255)
                                                                      gAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATc                       <  1:3670536/100‑1 (MQ=255)
                                                                      gAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATc                       <  1:2741224/100‑1 (MQ=255)
                                                                        ggCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCa                     >  2:831888/1‑100 (MQ=255)
                                                                        ggCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCa                     >  1:2930147/1‑100 (MQ=255)
                                                                         gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                    <  1:167339/100‑1 (MQ=255)
                                                                            aaaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTa                 <  2:2290395/100‑1 (MQ=255)
                                                                            aaaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTa                 >  1:78767/1‑100 (MQ=255)
                                                                             aagccATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAGCCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAc                <  2:3223703/97‑1 (MQ=255)
                                                                              agccgtCCCCGGGTTGATAGCCTACATGCAACGTGGCATGAAGCCTGGTGGTTACAACCTGATGGTGGCGGCGAAGGATACCGCCGATTATCCATgg                  <  2:2265818/92‑2 (MQ=255)
                                                                              aaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCGCTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAcc               >  2:3515140/1‑100 (MQ=255)
                                                                               aCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCg              >  1:3239477/1‑100 (MQ=255)
                                                                                 cATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTc            >  2:3202885/1‑100 (MQ=255)
                                                                                      ccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCtt       >  1:3794027/1‑100 (MQ=255)
                                                                                        gggTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTcc     >  2:3254101/1‑100 (MQ=255)
                                                                                         ggTTGATAGCCAATATGCAACGTTGCAGTAAACCTGGTGGGTAGAGGCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTccc    <  2:2176889/100‑1 (MQ=255)
                                                                                         ggTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATTTACCGTCGGCTTccc    >  2:1996518/1‑100 (MQ=255)
                                                                                                |                                                                                              
TGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGT  >  NC_000913/1501816‑1502006

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