Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,020,958 A→T 17.7% intergenic (‑39/‑180) sulA ← / → sxy SOS cell division inhibitor/CRP‑S‑dependent promoter expression factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,020,9580AT17.7% 312.5 / 56.2 129intergenic (‑39/‑180)sulA/sxySOS cell division inhibitor/CRP‑S‑dependent promoter expression factor
Reads supporting (aligned to +/- strand):  ref base A (61/45);  new base T (12/11);  total (73/56)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 154 total aligned reads displayed.

CAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAA  >  NC_000913/1020859‑1021057
                                                                                                   |                                                                                                   
cAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtaca                                                                                                     <  1:1173722/100‑1 (MQ=255)
  aTTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagt                                                                                                   <  1:448422/100‑1 (MQ=255)
  aTTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagt                                                                                                   <  1:569101/100‑1 (MQ=255)
   ttttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagta                                                                                                  <  2:1867607/100‑1 (MQ=255)
    tttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagtac                                                                                                 >  1:332498/1‑100 (MQ=255)
    tttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagtac                                                                                                 <  1:2275186/100‑1 (MQ=255)
     ttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagtaca                                                                                                <  2:2093484/100‑1 (MQ=255)
        cTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATcc                                                                                             >  2:461721/1‑100 (MQ=255)
          tgctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCat                                                                                           >  2:1296975/1‑99 (MQ=255)
          tgctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCat                                                                                           <  1:1476075/100‑2 (MQ=255)
           gctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGt                                                                                          >  2:2011910/1‑100 (MQ=255)
            ctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatac                                                                                         <  2:1893875/100‑4 (MQ=255)
            ctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatac                                                                                         <  1:1296970/100‑4 (MQ=255)
            ctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGGTACTGTATGGATGTACTGTACATCCatac                                                                                         <  1:1936459/100‑4 (MQ=255)
             tgcGGATGAGCGCGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGGTACTGTATGGATGTACTGTACATCCataca                                                                                        <  1:367124/100‑5 (MQ=255)
             tgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTg                                                                                         <  1:3548338/99‑1 (MQ=255)
              gcGGATGAGCACGACGAAGAACGATTTTTTTTGCCTGTAGTGTACATAATCACTCCAGCCCCTGTGACTTACTTTATGGATGTACAGTACCTCCATTGac                                                                                       <  1:2221041/100‑1 (MQ=255)
              gcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGGAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGac                                                                                       <  2:3506472/100‑1 (MQ=255)
              gcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatacag                                                                                       <  2:4353227/100‑6 (MQ=255)
              gcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatacag                                                                                       <  2:378757/100‑6 (MQ=255)
               cGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGaca                                                                                      >  1:112642/1‑100 (MQ=255)
                 gATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatacagtaa                                                                                    <  1:739098/100‑9 (MQ=39)
                    gagaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGATACTGTATAGATATACAGTACATCCAGTGACAACaaa                                                                                 >  1:4211392/1‑100 (MQ=255)
                     agaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAg                                                                                >  2:1021154/1‑100 (MQ=255)
                       aaCCACGAATAACGATGTGCATAGCCTGAAGTGTACATAATCAACCCAGCCCCTGTGATTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAg                                                                                <  1:2036041/98‑1 (MQ=255)
                        acgacgAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATc                                                                             >  1:3635875/1‑100 (MQ=255)
                        acgacgAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATc                                                                             <  2:4168431/100‑1 (MQ=255)
                           acgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAAc                                                                          <  1:2912875/100‑1 (MQ=255)
                           acgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAAc                                                                          >  2:537126/1‑100 (MQ=255)
                            cgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAAcc                                                                         >  2:3136165/1‑100 (MQ=255)
                                 aaCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTAtt                                                                    >  1:3155887/1‑100 (MQ=255)
                                     aTGTGCATATCCTGAAGTGTACATGATCCATACCGCCCATGTGAGTTACTGTATGTATGTACAGTACATCCATTTACAACAAAGATCAACCCTATTTTcn                                                                <  1:2275994/100‑2 (MQ=255)
                                     aTGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCg                                                                <  1:1404820/100‑1 (MQ=255)
                                      tgtgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCgg                                                               <  2:1785694/100‑1 (MQ=255)
                                          cATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAag                                                           <  1:350038/100‑1 (MQ=255)
                                          cATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAag                                                           >  2:1959171/1‑100 (MQ=255)
                                                 tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCg                                                    <  1:3270896/100‑1 (MQ=255)
                                                 tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCg                                                    <  1:461720/100‑1 (MQ=255)
                                                  gAAGTGTACATAATAAATCCAGCCCCTGTGAGTTACTGTATGGCTGTACAGTACATCCAGTGACAACAAAGATCAACCCTTTTTTCGGAAATAGCTTCg                                                    <  1:1131777/99‑1 (MQ=255)
                                                    aGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCaa                                                 >  1:3416773/1‑100 (MQ=255)
                                                    aGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCaa                                                 >  1:1193682/1‑100 (MQ=255)
                                                     gtgtACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCaaa                                                >  1:4099450/1‑100 (MQ=255)
                                                      tgtACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAt                                               >  2:3214031/1‑100 (MQ=255)
                                                        tACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCaaact                                              >  2:1126311/1‑97 (MQ=255)
                                                          cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                           <  1:4336950/100‑1 (MQ=255)
                                                          cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                           >  2:2334378/1‑100 (MQ=255)
                                                          cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                           <  2:657918/100‑1 (MQ=255)
                                                          cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                           <  1:94731/100‑1 (MQ=255)
                                                           aTAATCAAGCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGTTCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                           <  1:2802350/99‑1 (MQ=255)
                                                            tAATCAATCCAGCCCCTGTGAGTTACTTTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGt                                          <  1:1131327/99‑1 (MQ=255)
                                                            tAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTc                                         <  1:694264/100‑1 (MQ=255)
                                                            tAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTc                                         >  2:453478/1‑100 (MQ=255)
                                                            tAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTc                                         <  2:3059773/100‑1 (MQ=255)
                                                              atcaatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGt                                       >  2:2716217/1‑100 (MQ=255)
                                                              atcaatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGt                                       >  1:4187688/1‑100 (MQ=255)
                                                                caatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGtt                                      >  2:3851279/1‑99 (MQ=255)
                                                                caatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTg                                     >  1:2408288/1‑100 (MQ=255)
                                                                caatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTg                                     >  2:3323729/1‑100 (MQ=255)
                                                                 aatcCAACCCGTGTTGGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTgg                                    <  2:812458/100‑1 (MQ=255)
                                                                  atcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGt                                   >  1:1772935/1‑100 (MQ=255)
                                                                   tcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACCTCCAGTGACAACAACGATCAACCCTCTTTTCGTAAAGAGCCTCGCAAATTTTGTCGTTGGTg                                  >  2:1314557/1‑100 (MQ=255)
                                                                     cAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGAc                                >  2:3041054/1‑100 (MQ=255)
                                                                      aGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACg                               <  2:142557/100‑1 (MQ=255)
                                                                       gCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACgg                              >  2:4244935/1‑100 (MQ=255)
                                                                         cccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGa                            <  2:336708/100‑1 (MQ=255)
                                                                         cccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGa                            <  1:592550/100‑1 (MQ=255)
                                                                           cTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaaa                          >  2:2104961/1‑100 (MQ=255)
                                                                           cTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaaa                          >  1:1789457/1‑100 (MQ=255)
                                                                                aGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATa                     >  2:1826807/1‑100 (MQ=255)
                                                                                 gTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATaa                    >  2:3605932/1‑100 (MQ=255)
                                                                                  gttgtcactGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATaaa                   >  2:2169755/9‑100 (MQ=39)
                                                                                   tACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAt                  <  2:4288375/100‑1 (MQ=255)
                                                                                    aCTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAtt                 <  1:661416/100‑1 (MQ=255)
                                                                                    aCTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAtt                 <  2:3819826/100‑1 (MQ=255)
                                                                                    aCTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAt                  <  2:1772942/99‑1 (MQ=255)
                                                                                     gtcactGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTa                >  2:615592/6‑100 (MQ=255)
                                                                                     cTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTa                >  1:3808303/1‑100 (MQ=255)
                                                                                     cTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTa                >  1:4242567/1‑100 (MQ=255)
                                                                                      tcactGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa               >  2:2302145/5‑100 (MQ=255)
                                                                                      tcactGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa               >  1:3932332/5‑100 (MQ=255)
                                                                                      tcactGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa               >  2:4339975/5‑100 (MQ=255)
                                                                                      tGTATGGATGTACAGTACATCCAGTGACAACAAAGATGAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa               >  1:2357348/1‑100 (MQ=255)
                                                                                      tGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa               >  2:1931763/1‑100 (MQ=255)
                                                                                      tGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTCGGTGACGGGAAAACATAAATTaa               >  2:1863273/1‑100 (MQ=255)
                                                                                       cactGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt              >  1:4155333/4‑100 (MQ=255)
                                                                                        tATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc             <  1:1561515/100‑1 (MQ=255)
                                                                                        actGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc             >  2:3265702/3‑100 (MQ=255)
                                                                                        actGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc             >  1:1537766/3‑100 (MQ=255)
                                                                                           ggATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTg          >  1:270477/1‑100 (MQ=255)
                                                                                           ggATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTg          >  2:3976696/1‑100 (MQ=255)
                                                                                            gATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc         >  2:54769/1‑100 (MQ=255)
                                                                                            gATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc         >  1:4330483/1‑100 (MQ=255)
                                                                                             aTGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGcc        <  2:1516982/100‑1 (MQ=255)
                                                                                                tacagtagaTCCAGTGATAACAAGTATCAATCCTACTTTCGGAAAGAGCCTCGCAAATTGTGGCGTTGGTGACGGGAAAACATAAATTAATCTTGcc        <  1:4034343/97‑1 (MQ=255)
                                                                                                tacagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCt     >  1:3855109/1‑100 (MQ=255)
                                                                                                 acagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCtt    >  1:706943/1‑100 (MQ=255)
                                                                                                  cagtacaTCCAGTGACAACAAAGATCAACCCTGTTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTa   >  2:4180768/1‑100 (MQ=255)
                                                                                                   agtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTaa  >  2:4296646/1‑100 (MQ=255)
                                                                                                   agtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTaa  >  2:1745349/1‑100 (MQ=255)
                                                                                                   agtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTaa  >  2:2575987/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAA  >  NC_000913/1020859‑1021057

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: