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,0110CT99.1% 803.3 / ‑5.3 226intergenic (‑154/‑6)yceF/yceDhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (133/91);  total (134/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

Only 100 of 231 total aligned reads displayed.

AGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCT  >  NC_000913/1145921‑1146100
                                                                                          |                                                                                         
cgcaaCCCGCGATGTGTGCGAATTATTGGCAACCAGCACCCACAGGCTGCCTATTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgtgc                                                                                         <  2:7039952/90‑1 (MQ=255)
agcaaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgtgc                                                                                         >  2:6070958/4‑93 (MQ=255)
 gcaaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGcc                                                                                        <  2:5090781/91‑1 (MQ=255)
 gcaaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGcc                                                                                        >  2:3352466/3‑93 (MQ=255)
 gcaaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGcc                                                                                        >  1:7515756/3‑93 (MQ=255)
   aaCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCt                                                                                      <  1:3392802/93‑1 (MQ=255)
    aCCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTa                                                                                     >  2:991474/1‑93 (MQ=255)
     cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTAt                                                                                    >  1:5690874/1‑93 (MQ=255)
     cccGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTAt                                                                                    >  1:4228249/1‑93 (MQ=255)
       cGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGc                                                                                  >  1:7410989/1‑93 (MQ=255)
        gCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCa                                                                                 >  2:2850793/1‑93 (MQ=255)
        gCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCa                                                                                 >  1:3085462/1‑93 (MQ=255)
         cAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCaa                                                                                >  1:7802118/1‑93 (MQ=255)
          aaTGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCaaa                                                                               >  1:3397952/1‑93 (MQ=255)
           aTGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCaaaa                                                                              <  1:1177888/93‑1 (MQ=255)
           aTGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCaaaa                                                                              <  1:7916242/93‑1 (MQ=255)
           aTGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCaaaa                                                                              >  2:6358765/1‑93 (MQ=255)
            tgtgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAg                                                                             <  1:4623069/93‑1 (MQ=255)
            tgtgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAg                                                                             <  2:5235388/93‑1 (MQ=255)
              tgtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGt                                                                           <  1:416138/93‑1 (MQ=255)
               gtgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTa                                                                          >  1:2550268/1‑93 (MQ=255)
                tgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTaa                                                                         <  2:7122066/93‑1 (MQ=255)
                tgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTaa                                                                         <  2:1404007/93‑1 (MQ=255)
                tgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTaa                                                                         >  2:2456183/1‑93 (MQ=255)
                tgCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTaa                                                                         >  2:4650938/1‑93 (MQ=255)
                 gCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTaaa                                                                        <  2:4480779/93‑1 (MQ=255)
                  cGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTaaaa                                                                       <  1:2135358/93‑1 (MQ=255)
                   gAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAAt                                                                      >  1:4260364/1‑93 (MQ=255)
                    aaTTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAAtt                                                                     <  1:967772/93‑1 (MQ=255)
                       tattGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTAcc                                                                  >  1:6831891/1‑93 (MQ=255)
                       tattGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTAcc                                                                  >  1:5329927/1‑93 (MQ=255)
                         ttGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTAcccc                                                                >  2:3849590/1‑93 (MQ=255)
                            gCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGa                                                             <  1:5551767/93‑1 (MQ=255)
                              aaaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGAct                                                           >  1:6497664/1‑93 (MQ=255)
                              aaaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGAct                                                           >  2:3688369/1‑93 (MQ=255)
                               aaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGActc                                                          >  2:7104935/1‑93 (MQ=255)
                               aaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGActc                                                          >  1:2742612/1‑93 (MQ=255)
                               aaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGActc                                                          >  2:2717140/1‑93 (MQ=255)
                               aaaGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGActc                                                          <  2:4972952/93‑1 (MQ=255)
                                  ggCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCg                                                       >  1:2293819/1‑93 (MQ=255)
                                  ggCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCg                                                       >  2:7806482/1‑93 (MQ=255)
                                   gCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGa                                                      <  1:4443873/93‑1 (MQ=255)
                                    cAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGAt                                                     <  1:3328419/93‑1 (MQ=255)
                                     aaCCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATc                                                    <  2:5567060/93‑1 (MQ=255)
                                     aaCCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATc                                                    >  1:2667795/1‑93 (MQ=255)
                                      aCCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATcc                                                   <  1:1459422/93‑1 (MQ=255)
                                      aCCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCCGACTCTCGATcc                                                   <  1:7656913/93‑1 (MQ=255)
                                       ccACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCg                                                  >  2:169200/1‑93 (MQ=255)
                                        cacaGCCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCgg                                                 <  2:4331449/93‑1 (MQ=255)
                                         acaGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGt                                                <  2:5393869/93‑1 (MQ=255)
                                         acaGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGt                                                >  1:579117/1‑93 (MQ=255)
                                         acaGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGt                                                >  2:7074744/1‑93 (MQ=255)
                                           aGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTc                                              <  2:5028083/93‑1 (MQ=255)
                                           aGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTc                                              <  2:491141/93‑1 (MQ=255)
                                           aGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTc                                              <  1:3202530/93‑1 (MQ=255)
                                            ggatGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCg                                             >  2:6853226/4‑93 (MQ=255)
                                            ggCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCg                                             >  2:5505154/1‑93 (MQ=255)
                                             gCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGt                                            >  1:7102231/1‑93 (MQ=255)
                                              cTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTa                                           >  1:7071153/1‑93 (MQ=255)
                                              cTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTa                                           >  2:2684593/1‑93 (MQ=255)
                                              cTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTa                                           <  2:3162306/93‑1 (MQ=255)
                                              cTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTa                                           <  2:7532828/93‑1 (MQ=255)
                                               tGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTAc                                          >  2:5827011/1‑93 (MQ=255)
                                               tGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTAc                                          >  2:5267255/1‑93 (MQ=255)
                                               tGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTAc                                          <  1:6049660/93‑1 (MQ=255)
                                               tGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTAc                                          >  2:4866360/1‑93 (MQ=255)
                                                gCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACg                                         >  1:745456/1‑93 (MQ=255)
                                                 ccTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACgg                                        >  1:2067594/1‑93 (MQ=255)
                                                 ccTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACgg                                        >  1:4013194/1‑93 (MQ=255)
                                                  cTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGc                                       <  1:1377931/93‑1 (MQ=255)
                                                  cTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGc                                       <  1:1300965/93‑1 (MQ=255)
                                                  cTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGc                                       >  1:3491148/1‑93 (MQ=255)
                                                  cTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGc                                       >  1:1098974/1‑93 (MQ=255)
                                                   tttttCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACggat                                      >  2:1375282/1‑91 (MQ=255)
                                                   tttttCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCt                                      >  2:1500791/1‑93 (MQ=255)
                                                   tttttCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCt                                      >  1:3254750/1‑93 (MQ=255)
                                                    ttttCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTc                                     <  1:5153617/93‑1 (MQ=255)
                                                     tttCTTTGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCa                                    <  2:3000061/93‑1 (MQ=255)
                                                          ttGACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAAc                               >  1:4518628/1‑93 (MQ=255)
                                                            gACTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGc                             <  2:4487528/93‑1 (MQ=255)
                                                             aCTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGcc                            >  1:1849147/1‑93 (MQ=255)
                                                             aCTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGcc                            <  1:5406366/93‑1 (MQ=255)
                                                             aCTCTATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGcc                            <  2:4466021/93‑1 (MQ=255)
                                                              ctctATGACGTTACAAAGTTAATGTGCGTGCCCTATGCAAAAGGTAAAATTCCCCCTGCCTCTCGATCCGGTTCGTACGGCTCAAAAACGCCt                           <  1:1489180/93‑1 (MQ=255)
                                                               tctATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCtt                          >  2:4770570/1‑93 (MQ=255)
                                                               tctATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCtt                          >  2:3750913/1‑93 (MQ=255)
                                                                 tATGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGa                        <  1:4366753/93‑1 (MQ=255)
                                                                  aTGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTCCCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGAt                       <  2:5812336/93‑1 (MQ=255)
                                                                  aTGACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGAt                       >  2:3132443/1‑93 (MQ=255)
                                                                    gACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTa                     <  1:5859651/93‑1 (MQ=255)
                                                                    gACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTa                     <  1:3894338/93‑1 (MQ=255)
                                                                    gACGTTACAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTa                     >  2:459759/1‑93 (MQ=255)
                                                                           cAAAGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGt              >  2:5905274/1‑93 (MQ=255)
                                                                            aaaGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTa             <  1:3936077/93‑1 (MQ=255)
                                                                            aaaGTTAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTa             >  2:1999319/1‑93 (MQ=255)
                                                                                 tAATATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCtat        >  1:6974089/1‑93 (MQ=255)
                                                                                  aaTATGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCtata       >  2:859603/1‑93 (MQ=255)
                                                                                   atatGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAc      >  1:3656012/1‑93 (MQ=255)
                                                                                   atatGCGCGCCCTATGCAAAAGGTAAAATTACCACTGACTCTCGATCCGGTTCGTACGGTTCAAAAACGACTTGATTACCAGGGTATCTATAc      >  1:6711820/1‑93 (MQ=255)
                                                                                    tatGCGTGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATAcc     >  1:7109761/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCGCCCTATGCAAAAGGTAAAATTACCCCTGACTCTCGATCCGGTTCGTACGGCTCAAAAACGCCTTGATTACCAGGGTATCTATACCCCT  >  NC_000913/1145921‑1146100

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