Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,020,958 A→T 15.3% 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,9580AT15.3% 317.1 / 43.4 123intergenic (‑39/‑180)sulA/sxySOS cell division inhibitor/CRP‑S‑dependent promoter expression factor
Reads supporting (aligned to +/- strand):  ref base A (58/46);  new base T (10/9);  total (68/55)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

Only 100 of 148 total aligned reads displayed.

CAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAA  >  NC_000913/1020859‑1021057
                                                                                                   |                                                                                                   
cAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtaca                                                                                                     >  2:1811398/1‑100 (MQ=255)
  aTTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagt                                                                                                   <  1:1176056/100‑1 (MQ=255)
  aTTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagt                                                                                                   <  2:1975474/100‑1 (MQ=255)
  aTTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagt                                                                                                   >  1:2047589/1‑100 (MQ=255)
   ttttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtactgta                                                                                                  <  2:1123011/100‑1 (MQ=255)
   ttttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagta                                                                                                  <  1:467633/100‑1 (MQ=255)
   ttttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagta                                                                                                  <  2:535757/100‑1 (MQ=255)
   ttttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagta                                                                                                  <  1:2494705/100‑1 (MQ=255)
   ttttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagta                                                                                                  <  2:3597916/100‑1 (MQ=255)
        cTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATcc                                                                                             <  1:950441/100‑1 (MQ=255)
          taccgcGAATGAGAACGACGAAGAACCATGTGCGTGGCCCCAAGTCTACTTATTCCATCCAGGCCCTGTGACTTACTAGATGGATGTACTGTACATcctt                                                                                           <  1:541318/96‑3 (MQ=255)
           gctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGt                                                                                          >  2:892861/1‑100 (MQ=255)
           gctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGt                                                                                          >  1:335431/1‑100 (MQ=255)
           gctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGt                                                                                          >  2:830062/1‑100 (MQ=255)
            ctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatac                                                                                         <  2:1398188/100‑4 (MQ=255)
            ctgcGGATGAGAACGACGAAAAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTg                                                                                         <  1:1578584/100‑1 (MQ=255)
             tgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCataca                                                                                        <  1:2303159/100‑5 (MQ=255)
             tgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCataca                                                                                        <  2:2127798/100‑5 (MQ=255)
               cGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatacagt                                                                                      <  1:2478521/100‑7 (MQ=255)
                ggATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGacaa                                                                                     <  1:1530532/100‑1 (MQ=255)
                  aTGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatacagtaac                                                                                   <  1:3392098/100‑10 (MQ=38)
                  aTGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGacaaca                                                                                   <  1:1344545/100‑1 (MQ=255)
                  aTGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGacaac                                                                                    <  2:1443095/99‑1 (MQ=255)
                   tGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGAcaacaa                                                                                  >  1:3543876/1‑100 (MQ=255)
                    gaggACGACGAAGAACGATGTGCATAGGCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACATTACATCCAGTGaca                                                                                      <  1:1262596/95‑1 (MQ=255)
                     agaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAg                                                                                <  1:345785/100‑1 (MQ=255)
                     agaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAg                                                                                <  2:3118581/100‑1 (MQ=255)
                     agaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAg                                                                                >  2:1952707/1‑100 (MQ=255)
                        acgacgAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATc                                                                             <  2:457118/100‑1 (MQ=255)
                         cgacgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCACTGTGAGTTACTTTATGGATGTACAGTACATCCAGTGACAACAAAGATCa                                                                            <  1:2182299/100‑1 (MQ=255)
                           acgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAAc                                                                          <  1:2862104/100‑1 (MQ=255)
                            cgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAAcc                                                                         <  1:2973339/100‑1 (MQ=255)
                             gaagaaCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAAccc                                                                        <  1:3365253/100‑1 (MQ=255)
                                  aCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTAttt                                                                   >  1:2629125/1‑100 (MQ=255)
                                  aCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTAttt                                                                   >  1:3642836/1‑100 (MQ=255)
                                  aCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTAttt                                                                   >  2:524044/1‑100 (MQ=255)
                                    gATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTc                                                                 >  2:166960/1‑100 (MQ=255)
                                     aTGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCg                                                                <  2:799533/100‑1 (MQ=255)
                                       gtgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGa                                                              <  2:2781320/100‑1 (MQ=255)
                                       gtgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGa                                                              <  1:524038/100‑1 (MQ=255)
                                         gCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaaa                                                            >  2:3500621/1‑100 (MQ=255)
                                          cATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAag                                                           >  2:331136/1‑100 (MQ=255)
                                            tAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAagag                                                         <  1:1424999/100‑1 (MQ=255)
                                            tAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACAGCCAGTGACAACAAAGATCAACCCTATTTTCGGAAagag                                                         >  2:1174988/1‑100 (MQ=255)
                                              gcttGAAGTGTACATAATCAATCTAGCCGCTGTGAGTTACTGTATGGATGTAGAGTACATCAAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGc                                                        >  1:2453669/4‑99 (MQ=255)
                                              gCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGcc                                                       >  1:1169493/1‑100 (MQ=255)
                                                cTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTc                                                     >  2:2859451/1‑100 (MQ=255)
                                                cTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTc                                                     <  2:2476361/100‑1 (MQ=255)
                                                cTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTc                                                     >  1:1361743/1‑100 (MQ=255)
                                                     gtgtACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCaaa                                                <  1:1781293/100‑1 (MQ=255)
                                                       gtACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAtt                                              >  1:1030037/1‑100 (MQ=255)
                                                       gtACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAtt                                              >  1:1027607/1‑100 (MQ=255)
                                                          cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                           >  2:3770723/1‑100 (MQ=255)
                                                          cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                           <  2:2277889/100‑1 (MQ=255)
                                                          cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                           >  1:3659791/1‑100 (MQ=255)
                                                          cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                           >  1:3547645/1‑100 (MQ=255)
                                                          cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                           <  2:941189/100‑1 (MQ=255)
                                                           aTAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGt                                          >  2:437227/1‑100 (MQ=255)
                                                              atcaatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGt                                       >  1:2154711/1‑100 (MQ=255)
                                                                 aatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTgg                                    >  2:3035068/1‑100 (MQ=255)
                                                                 aatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTgg                                    >  2:1928981/1‑100 (MQ=255)
                                                                  atcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGt                                   >  2:3390410/1‑100 (MQ=255)
                                                                       gCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTTGCAAATTTTGTCGTTGGTGACgg                              >  1:895881/1‑100 (MQ=255)
                                                                       gCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACgg                              <  1:2816341/100‑1 (MQ=255)
                                                                       gCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACgg                              <  1:196027/100‑1 (MQ=255)
                                                                        ccccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGACCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACggg                             >  1:2978501/1‑100 (MQ=255)
                                                                         cccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGa                            <  1:2570474/100‑1 (MQ=255)
                                                                           cTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaaa                          >  2:1156430/1‑100 (MQ=255)
                                                                            tgtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaaaa                         <  2:2176368/100‑1 (MQ=255)
                                                                               gAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACAt                      >  1:3806224/1‑100 (MQ=255)
                                                                               gAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACAt                      >  2:774860/1‑100 (MQ=255)
                                                                                aGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATa                     <  2:998615/100‑1 (MQ=255)
                                                                                 gTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATaa                    >  1:272742/1‑100 (MQ=255)
                                                                                  ttACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATaaa                   <  1:1464371/100‑1 (MQ=255)
                                                                                      tGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa               >  2:678211/1‑100 (MQ=255)
                                                                                       gTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt              <  1:3166845/100‑1 (MQ=255)
                                                                                       gTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTCGTGACGGAAAAACATaaataaat              >  1:2695549/1‑100 (MQ=255)
                                                                                       ctctGGATGTACTGTACATCCAGTGACAACAGAGATCAACCCAATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt              >  2:2579872/4‑100 (MQ=255)
                                                                                       cactGGATGTACTGTACATCCAGTGAAAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt              >  2:2318685/4‑100 (MQ=255)
                                                                                        actGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc             >  1:1931773/3‑100 (MQ=255)
                                                                                        actGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc             >  1:1396280/3‑100 (MQ=255)
                                                                                         ctGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATaaaagatcgg            >  1:783791/2‑93 (MQ=255)
                                                                                         ctGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCt            >  1:2635457/2‑100 (MQ=255)
                                                                                         aTGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCt            <  2:2929882/100‑1 (MQ=255)
                                                                                          tGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCtt           >  1:3151332/1‑100 (MQ=255)
                                                                                          tGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCtt           >  1:1768278/1‑100 (MQ=255)
                                                                                           ggATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTg          >  1:248366/1‑100 (MQ=255)
                                                                                            gATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCACATTTTCTCGTTGGTGACGGGAACACATAAATTAATCt            >  1:658657/1‑97 (MQ=255)
                                                                                            gATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc         >  2:3736037/1‑100 (MQ=255)
                                                                                             aTGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGcc        <  2:904915/100‑1 (MQ=255)
                                                                                              tGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGTCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGccc       <  1:3750878/100‑1 (MQ=255)
                                                                                               gtacagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGcccc      >  1:1132464/1‑100 (MQ=255)
                                                                                                tactgtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCt     >  1:1295917/1‑100 (MQ=255)
                                                                                                tacagtacaTCCAGTGACAAGAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCt     >  1:2325032/1‑100 (MQ=255)
                                                                                                tacagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCt     >  2:2724317/1‑100 (MQ=255)
                                                                                                tacagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCt     >  2:1688423/1‑100 (MQ=255)
                                                                                                tacagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCt     >  1:890603/1‑100 (MQ=255)
                                                                                                 acagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCtt    >  1:2260985/1‑100 (MQ=255)
                                                                                                   agtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTaa  >  1:1685605/1‑100 (MQ=255)
                                                                                                   agtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTaa  >  2:2715077/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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: