Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,853 T→C 17.5% intergenic (+64/‑280) bolA → / → tig stationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913454,8530TC17.5% 321.3 / 60.7 150intergenic (+64/‑280)bolA/tigstationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)
Reads supporting (aligned to +/- strand):  ref base T (60/63);  new base C (15/11);  total (75/75)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

Only 100 of 151 total aligned reads displayed.

CTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATA  >  NC_000913/454755‑454951
                                                                                                  |                                                                                                  
ctcctcCCTGTCGTGGGGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTg                                                                                                   >  2:3781593/1‑100 (MQ=255)
 tcctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTTCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGc                                                                                                  <  2:2117616/100‑1 (MQ=255)
 tcctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGc                                                                                                  <  2:2118530/100‑1 (MQ=255)
  cctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCg                                                                                                 >  1:706517/1‑100 (MQ=255)
    tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                               >  1:1126617/1‑100 (MQ=255)
     cccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGc                                                                                              >  2:90198/1‑100 (MQ=255)
       cTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCg                                                                                            >  1:461909/1‑100 (MQ=255)
            gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttttt                                                                                       <  2:3679014/100‑2 (MQ=255)
              ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTTTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                     <  2:929444/100‑1 (MQ=255)
              ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                     <  1:3566894/100‑4 (MQ=255)
              ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                     <  2:2010960/100‑4 (MQ=255)
              ggAGCACGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                     <  2:3904880/100‑4 (MQ=255)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAACAACGGCCTGCGGGCCgttttgttt                                                                                    >  2:1441117/1‑100 (MQ=255)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                    <  1:3729456/100‑5 (MQ=39)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                    <  1:605262/100‑5 (MQ=255)
                aGCGGGAAGCATGGCGTCACAACGCATTTGCACCTGTCGGCGCTTTTCCAGTATGTTGTTAAAGATTTTATGAAAAACAGCCAGCGGGCCgttgtgtttt                                                                                   <  2:1166892/100‑1 (MQ=255)
                aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                   <  1:3192164/100‑1 (MQ=255)
                aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                   >  2:2576964/1‑100 (MQ=255)
                 gCAGGAAGCATCTCGTAAAAACGCATTTGAAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcata                                                                                   <  2:496051/99‑6 (MQ=38)
                 gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAACAACGGCCTGCGGGCCGttttgttttg                                                                                  >  2:656483/1‑100 (MQ=255)
                 gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                  >  1:996136/1‑100 (MQ=255)
                 gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                  <  1:3771528/100‑1 (MQ=255)
                 gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                  <  2:640682/100‑1 (MQ=255)
                    ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                               >  2:2948870/1‑100 (MQ=255)
                    ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGACGTTTTGTTTTGTc                                                                                >  1:2700869/1‑99 (MQ=255)
                    ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGACCGTTTTGTTTTGTCt                                                                               >  2:2303289/1‑100 (MQ=255)
                      aaGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                             <  1:2168075/100‑1 (MQ=255)
                       aGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGa                                                                            <  2:1264273/100‑1 (MQ=255)
                        gCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAt                                                                           <  1:1131447/100‑1 (MQ=255)
                              cgTAAAACCGCATTTGCAACTGTCGGCGCTTTTCCAGTCTGTTGCTAAAGATTTTATGAAAAACGACCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                      >  2:2073796/1‑99 (MQ=255)
                              cgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                      >  1:3344618/1‑99 (MQ=255)
                               gTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgcgc                                                                    <  2:125026/100‑1 (MQ=255)
                                tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCt                                                                   <  1:365651/100‑1 (MQ=255)
                                 aaaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCtt                                                                  <  1:48508/100‑1 (MQ=255)
                                      cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                             <  2:3419210/100‑1 (MQ=255)
                                      cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                             <  1:352225/100‑1 (MQ=255)
                                      cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                             >  2:1086182/1‑100 (MQ=255)
                                        cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGCTTTTGCGCTTTTTGccc                                                           >  1:2705265/1‑100 (MQ=255)
                                        cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                           <  1:1854469/100‑1 (MQ=255)
                                        cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                           <  1:137776/100‑1 (MQ=255)
                                         aTTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCa                                                          <  1:2417155/100‑1 (MQ=255)
                                          tttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAg                                                         <  1:2933994/100‑1 (MQ=255)
                                             gCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAt                                                      >  1:2857717/1‑100 (MQ=255)
                                               aaCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTc                                                    <  1:3377178/100‑1 (MQ=255)
                                               aaCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTc                                                    <  2:1905588/100‑1 (MQ=255)
                                                  tGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGa                                                 <  1:847851/100‑1 (MQ=255)
                                                     cGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGa                                              >  2:1103860/1‑100 (MQ=255)
                                                     cGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGa                                              >  1:1961899/1‑100 (MQ=255)
                                                         gcTTTTCCAGTATGTTGCTAAAGATTTTATGAGAATCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAAt                                          <  2:2391162/100‑1 (MQ=255)
                                                         gcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAAt                                          <  2:3615701/100‑1 (MQ=255)
                                                           ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                        <  2:2969437/100‑1 (MQ=255)
                                                           ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                        >  2:340601/1‑100 (MQ=255)
                                                            tttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGc                                       >  2:1675989/1‑100 (MQ=255)
                                                             ttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGcc                                      >  1:675924/1‑100 (MQ=255)
                                                             ttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGcc                                      <  1:256960/100‑1 (MQ=255)
                                                              tCCAGTAAGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGccc                                     >  1:2076826/1‑100 (MQ=255)
                                                               ccAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCTGCATTCAGACGAAAATTGCCCg                                    >  2:3285477/1‑100 (MQ=255)
                                                                 aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCggg                                  <  2:1088125/100‑1 (MQ=255)
                                                                 aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCggg                                  <  1:1763247/100‑1 (MQ=255)
                                                                     tGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                              >  1:3864682/1‑100 (MQ=255)
                                                                     ngTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCAAGCATTCAGACGAAAATTGCCCGGGAAtt                              >  2:2816480/2‑100 (MQ=255)
                                                                       ttGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgt                            <  1:432606/100‑1 (MQ=255)
                                                                        tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                           <  1:2902219/100‑1 (MQ=255)
                                                                         gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGGTTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                          >  1:1267438/1‑100 (MQ=255)
                                                                         gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                          >  2:2409788/1‑100 (MQ=255)
                                                                          cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaa                         <  1:3296330/100‑1 (MQ=255)
                                                                          cTAAAGATTTTATGAAAAACGGCATGCGGGCCGTTTTGTTTTGTCTGGATTTTCCGCTTTTTGACAAGCATTCAGACGAAAATTGCTCGGGAATTGTGaa                         >  2:1918330/1‑100 (MQ=255)
                                                                            aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                       <  1:2267313/100‑1 (MQ=255)
                                                                            aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTCCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                       >  2:2282534/1‑100 (MQ=255)
                                                                             aaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaa                      <  1:3115300/100‑1 (MQ=255)
                                                                              aGATTTTATGAAAATCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaaa                     <  1:2349663/100‑1 (MQ=255)
                                                                               gATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAAt                    <  1:8556/100‑1 (MQ=255)
                                                                               gATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAAt                    <  2:773504/100‑1 (MQ=255)
                                                                                  tttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGTCGAAAATTGCCCGGGAATTGTGAAAAAATAcg                 <  1:3162594/100‑1 (MQ=255)
                                                                                  tttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGTCTGGGTATTGTGAAAATATAcg                 <  2:3344596/100‑1 (MQ=255)
                                                                                  tttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcg                 >  1:766611/1‑100 (MQ=255)
                                                                                  tttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcg                 <  2:921216/100‑1 (MQ=255)
                                                                                   caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCCGCATTAAGACGCAAATTGCCCGGGAATTGTGAAAAAATa                   >  2:1448988/7‑97 (MQ=37)
                                                                                   caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgc                >  2:1487212/7‑100 (MQ=38)
                                                                                     aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTTTTTTTTATTTAGACGAAAATTGCCCGGGAATTGTGAAAAAATAtgcg               >  2:2795083/5‑99 (MQ=19)
                                                                                     aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa              >  1:1211274/5‑100 (MQ=39)
                                                                                     aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa              >  1:2565954/5‑100 (MQ=39)
                                                                                      tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGCCGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa              >  1:1337853/1‑99 (MQ=255)
                                                                                      tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc             >  1:3066393/1‑100 (MQ=255)
                                                                                      aacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc             >  1:3578256/4‑100 (MQ=255)
                                                                                      aacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGAGTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc             >  2:2106383/4‑100 (MQ=255)
                                                                                       gAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa            >  1:3100669/1‑100 (MQ=255)
                                                                                       gAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa            >  2:1006368/1‑100 (MQ=255)
                                                                                       acaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa            >  1:1415609/3‑100 (MQ=255)
                                                                                        caaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCTGACGAAAATTGCACGGGAATTGTGAAAAAATACGCGACAg           >  1:622/2‑100 (MQ=255)
                                                                                        caaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAg           >  1:2126619/2‑100 (MQ=255)
                                                                                        aaaaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAg           >  1:1183048/1‑100 (MQ=255)
                                                                                         aaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgc          >  2:2048667/1‑100 (MQ=255)
                                                                                          aaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcg         >  2:948965/1‑100 (MQ=255)
                                                                                           naCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgc        >  1:2283504/2‑100 (MQ=255)
                                                                                           aaCGGCCTGCGGGCCGTGTTGTTTTGTCTGGATTTTGCGTTTTTTGCCCAGCATTCAGCCGAAAATTGCCCGGGAATTTTGAAAAAATACGCGACAgcg         <  2:757586/99‑1 (MQ=255)
                                                                                              ggCCTGGGGGCCGTGTTGTTTTGTCTGGATTTTGCGCTTTTTGCCGAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgcgc      <  2:2740899/99‑1 (MQ=255)
                                                                                              ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCGCAGCATTAAGGCGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgcg       <  1:1770746/98‑1 (MQ=255)
                                                                                              ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa     <  1:3023609/100‑1 (MQ=255)
                                                                                                 cTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATa  >  1:1874792/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATA  >  NC_000913/454755‑454951

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: