Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,146,011 C→T intergenic (‑154/‑6) yceF ← / → yceD hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,146,0110CT100.0% 546.0 / NA 150intergenic (‑154/‑6)yceF/yceDhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (76/74);  total (76/74)

Only 100 of 155 total aligned reads displayed.

AGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCT  >  NC_000913/1145921‑1146100
                                                                                          |                                                                                         
agcaaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgtgc                                                                                         >  1:100719/4‑93 (MQ=255)
 gcaaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGcc                                                                                        <  1:3016495/91‑1 (MQ=255)
 gcaaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGcc                                                                                        >  2:4639156/3‑93 (MQ=255)
  caaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGccc                                                                                       <  2:3956596/92‑1 (MQ=255)
   aaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCt                                                                                      <  2:743562/93‑1 (MQ=255)
   aaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCt                                                                                      >  2:319656/1‑93 (MQ=255)
    aCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACTAAGTTAATATGCGTGCCCTa                                                                                     <  2:2666496/93‑1 (MQ=255)
    aCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTa                                                                                     <  2:1421379/93‑1 (MQ=255)
     cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTAt                                                                                    <  1:2602409/93‑1 (MQ=255)
     cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTAt                                                                                    >  1:249647/1‑93 (MQ=255)
     cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTAt                                                                                    <  1:3039809/93‑1 (MQ=255)
     cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTAt                                                                                    <  2:292561/93‑1 (MQ=255)
     cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTAt                                                                                    >  2:1527861/1‑93 (MQ=255)
      ccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATg                                                                                   >  2:2880671/1‑93 (MQ=255)
       cGCAATGTATGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGc                                                                                  >  1:3185851/1‑93 (MQ=255)
        gCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCa                                                                                 >  2:1526573/1‑93 (MQ=255)
         cAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCaa                                                                                >  1:4455541/1‑93 (MQ=255)
           aTGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCaaaa                                                                              <  2:3872047/93‑1 (MQ=255)
            tgtgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAg                                                                             <  2:3246827/93‑1 (MQ=255)
            tgtgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAg                                                                             <  2:1160874/93‑1 (MQ=255)
            tgtgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAg                                                                             >  2:1650137/1‑93 (MQ=255)
            tgtgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAg                                                                             <  1:1574466/93‑1 (MQ=255)
             gtgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAgg                                                                            <  1:4556539/93‑1 (MQ=255)
             gtgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAgg                                                                            >  2:2514834/1‑93 (MQ=255)
              tgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGt                                                                           <  1:3209787/93‑1 (MQ=255)
               gtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTa                                                                          <  1:4263018/93‑1 (MQ=255)
               gtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTa                                                                          <  1:1559176/93‑1 (MQ=255)
               gtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTa                                                                          <  1:596153/93‑1 (MQ=255)
                  cGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTaaaa                                                                       >  2:2121043/1‑93 (MQ=255)
                  cGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTaaaa                                                                       <  2:2539185/93‑1 (MQ=255)
                  cGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTaaaa                                                                       <  2:2301761/93‑1 (MQ=255)
                  cGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTaaaa                                                                       <  2:653435/93‑1 (MQ=255)
                    aaTTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAAtt                                                                     >  1:883232/1‑93 (MQ=255)
                      ttattGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTAc                                                                   >  1:486794/1‑93 (MQ=255)
                      ttattGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTAc                                                                   >  1:1823938/1‑93 (MQ=255)
                      ttattGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTAc                                                                   <  1:4009831/93‑1 (MQ=255)
                      ttattGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTAc                                                                   >  1:2205975/1‑93 (MQ=255)
                         ttGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTAcccc                                                                <  1:1474458/93‑1 (MQ=255)
                            gCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGa                                                             <  2:590185/93‑1 (MQ=255)
                             cAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGAc                                                            >  1:1291813/1‑93 (MQ=255)
                             cAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGAc                                                            <  1:982202/93‑1 (MQ=255)
                              aaaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGAct                                                           >  2:3486693/1‑93 (MQ=255)
                               aaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGActc                                                          <  2:1111902/93‑1 (MQ=255)
                                 aGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACtctc                                                        >  1:2212006/1‑93 (MQ=255)
                                    cAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGAt                                                     <  1:833119/93‑1 (MQ=255)
                                    cAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGAt                                                     <  1:999170/93‑1 (MQ=255)
                                    cAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGAt                                                     <  2:4232835/93‑1 (MQ=255)
                                     aaCCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATc                                                    <  1:4606/93‑1 (MQ=255)
                                       ccACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTCAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCg                                                  >  2:513982/1‑93 (MQ=255)
                                        cacaGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCgg                                                 >  1:1840317/1‑93 (MQ=255)
                                        cacaGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCgg                                                 <  2:2296015/93‑1 (MQ=255)
                                           aGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAACGGTAAAATTAACTCTGACTCTCGTTCCGGTTc                                              >  2:2315857/1‑93 (MQ=255)
                                           aGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTc                                              >  2:473518/1‑93 (MQ=255)
                                           aGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGCAAAATTACCCCTGACTCTCGATCCGGTTc                                              >  1:3636481/1‑93 (MQ=255)
                                            ggCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCg                                             >  2:3077229/1‑93 (MQ=255)
                                             gCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGt                                            <  2:4301246/93‑1 (MQ=255)
                                                gCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACg                                         >  2:141861/1‑93 (MQ=255)
                                                gCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACg                                         >  1:182438/1‑93 (MQ=255)
                                                 ccTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACgg                                        >  2:3618933/1‑93 (MQ=255)
                                                 ccTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACgg                                        <  1:1030637/93‑1 (MQ=255)
                                                 ccTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACgg                                        >  1:2883241/1‑93 (MQ=255)
                                                   tttttCTTTGACTCTATGAGGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTCCCCCTGACTCACGATCCGGTTCGTACGGCt                                      <  2:4550369/93‑1 (MQ=255)
                                                   tttttCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGcc                                      >  1:37146/1‑92 (MQ=255)
                                                   tttttCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGcc                                      >  1:1127243/1‑92 (MQ=255)
                                                   tttttCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCt                                      >  2:1009026/1‑93 (MQ=255)
                                                      ttCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCaa                                   >  1:1435103/1‑93 (MQ=255)
                                                       tCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCaaa                                  <  1:486883/93‑1 (MQ=255)
                                                        cTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCaaaa                                 >  1:278806/1‑93 (MQ=255)
                                                         tttGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCaaaaa                                >  1:1052162/1‑93 (MQ=255)
                                                         tttGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCaaaaa                                <  2:3837631/93‑1 (MQ=255)
                                                         tttGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCaaaaa                                <  1:1768061/93‑1 (MQ=255)
                                                         tttGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCaaaaa                                <  1:4032884/93‑1 (MQ=255)
                                                          ttGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAAc                               >  2:946119/1‑93 (MQ=255)
                                                              ctctATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCt                           <  2:3246552/93‑1 (MQ=255)
                                                               tctATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCtt                          >  1:564064/1‑93 (MQ=255)
                                                               tctATGACGTTACAAAGTTAATATGAGTGACCTATGCAAAAGGTAAAATTAACACTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCtt                          >  1:3885865/1‑93 (MQ=255)
                                                                ctATGACGTTACAAAGTTAATTTGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCCATCCGGTTCGTACGGCTCCAAAACGCCTTg                         >  2:415483/1‑93 (MQ=255)
                                                                 tATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGa                        <  2:1817696/93‑1 (MQ=255)
                                                                 tATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACCGCTCAAAAACGCCTTGa                        >  1:714048/1‑93 (MQ=255)
                                                                  aTGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGAt                       <  1:3596418/93‑1 (MQ=255)
                                                                  aTGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGAt                       <  1:2325177/93‑1 (MQ=255)
                                                                     aCGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTAc                    <  1:3722509/93‑1 (MQ=255)
                                                                     aCGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTAc                    <  1:2111939/93‑1 (MQ=255)
                                                                        ttACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAg                 <  2:821510/93‑1 (MQ=255)
                                                                        ttACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTACGTACGGCTCAAAAACGCCTTGATTACCAg                 >  1:3914306/1‑93 (MQ=255)
                                                                         tACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAgg                >  2:4161790/1‑93 (MQ=255)
                                                                          aCAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAggg               <  2:3988945/93‑1 (MQ=255)
                                                                          aCAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAggg               >  2:2336725/1‑93 (MQ=255)
                                                                           cAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGt              >  1:3533985/1‑93 (MQ=255)
                                                                           cAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGt              >  2:2168880/1‑93 (MQ=255)
                                                                            aaaGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGGTCAAAAACGCCTTGATTACCAGGGTa             >  1:3135082/1‑93 (MQ=255)
                                                                            aaaGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTa             <  2:1728042/93‑1 (MQ=255)
                                                                               gTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCt          <  2:2927217/93‑1 (MQ=255)
                                                                                ttAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCta         >  1:4156486/1‑93 (MQ=255)
                                                                                ttAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCta         <  1:2538289/93‑1 (MQ=255)
                                                                                  aaTATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCtata       <  2:793119/93‑1 (MQ=255)
                                                                                    tatGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAcc     >  1:755748/1‑93 (MQ=255)
                                                                                      tGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAcccc   <  1:1448620/93‑1 (MQ=255)
                                                                                      tGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAcccc   >  2:3558857/1‑93 (MQ=255)
                                                                                       gcgtgcCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCt  >  2:2804015/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCT  >  NC_000913/1145921‑1146100

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