Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,016 T→C 6.7% intergenic (+32/‑28) secA → / → mutT protein translocation ATPase/8‑oxo‑dGTP diphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,0160TC 6.7% 205.9 / 5.0 80intergenic (+32/‑28)secA/mutTprotein translocation ATPase/8‑oxo‑dGTP diphosphatase
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base C (2/4);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

CGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTA  >  NC_000913/110876‑111158
                                                                                                                                            |                                                                                                                                              
cccAGCTGCACCGGCGGCGCAAACCGGCGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCt                                                                                                                                              <  2:645278/141‑1 (MQ=255)
 gCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCCTGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTg                                                                                                                                             >  2:534036/1‑143 (MQ=255)
      tgaaCTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGGTCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCt                                                                                                                                        <  1:772786/140‑1 (MQ=255)
      cgCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGAACATTCGCTTTACTGGTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCt                                                                                                                                        <  1:252462/142‑1 (MQ=255)
       gCACTGGAGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCtt                                                                                                                                       <  1:897918/143‑1 (MQ=255)
         aCTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCtttt                                                                                                                                     >  1:952593/1‑143 (MQ=255)
          cTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCtttta                                                                                                                                    <  2:872346/143‑2 (MQ=255)
               ggcgCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGCAGTAAAAGGCGCAGGATTCTGCGCCTTCTTTATAg                                                                                                                               >  2:83469/1‑143 (MQ=255)
                     aaCCGGAGAGCGCAAAGTAGGACGTAACGCTCCTTGCCCCTGCGGATCTGGTAAAAAATACCAGCAGTGCGCTGGCCGCATGCAATAAAAGCTCAGTTTTGAAGTAACAGGCGCAGGAGTCTGCGCCTTTTTTGTAGGTTCaa                                                                                                                         >  2:446383/1‑143 (MQ=255)
                         ggAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAACTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACAGTTGAAGTAAAAGGCGCCGGCTTCTGCGCCTTCTTTATAGGATTAAGACa                                                                                                                     >  2:18560/1‑143 (MQ=255)
                              gcgcAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaa                                                                                                                >  2:913964/1‑143 (MQ=255)
                               cgcAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCAGTGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaaa                                                                                                               <  2:663761/143‑1 (MQ=255)
                                gcAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaaaa                                                                                                              >  1:1050146/1‑143 (MQ=255)
                                 cAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaaaaa                                                                                                             >  2:684646/1‑143 (MQ=255)
                                       aGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTAATAAAACGTGCTGCATTACACGAAAACTGTGGAAGTAAAAGGCGCCGGATTCTGCGCCTTTTTTATAGGTTTACGACAATGAAAACGCTGCa                                                                                                       >  2:1077569/1‑143 (MQ=255)
                                        ggACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGTCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAACGAAAAAGCTGCaa                                                                                                      >  2:394982/1‑143 (MQ=255)
                                         gACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCaaa                                                                                                     <  2:462551/143‑1 (MQ=255)
                                            gTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAAAATGAAAAAGCTGCCAAtcg                                                                                                  >  2:1097694/1‑141 (MQ=255)
                                            gTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGACTTTTTTATAGGTTTAAGAAAATAAAAAAGCTGCAAATTg                                                                                                  >  2:937656/1‑143 (MQ=255)
                                             taacagtCCGTGGCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGc                                                                                                 <  1:822436/137‑1 (MQ=255)
                                              aaCGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                                                                <  1:611729/143‑1 (MQ=255)
                                              aaCGATCCCAGCCCGAGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                                                                <  1:796859/143‑1 (MQ=255)
                                               aCGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAAATGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAGAGCTGCAAATTGCgg                                                                                               >  1:25176/1‑143 (MQ=255)
                                                cGATCCTGGCCCGTGAGGTTCTGGTAAAAAATACTAGCGGGGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggt                                                                                              <  2:328557/143‑1 (MQ=255)
                                                 gATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATACAAGCTAACAGTTGAAGTAAAAGGCGCAGGACTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAACAGCTGCAACTTGCggta                                                                                             >  2:217712/1‑143 (MQ=255)
                                                 gATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggta                                                                                             >  1:993117/1‑143 (MQ=255)
                                                 gATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggta                                                                                             >  1:163982/1‑143 (MQ=255)
                                                 gATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTGTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggga                                                                                             >  2:45698/1‑141 (MQ=255)
                                                     cTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAACTTGCgggaggta                                                                                         >  2:922380/1‑143 (MQ=255)
                                                      ttccccGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGtat                                                                                        <  1:478187/140‑1 (MQ=255)
                                                       tGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCCGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAACGGCTGCCACTTGCGGTAGGtatt                                                                                       >  2:788796/1‑143 (MQ=255)
                                                          ccGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTattatt                                                                                    <  1:1088861/143‑1 (MQ=255)
                                                            gtagcgTTCTCTTAACAAATATAAGCAATGCCATGGCTGCCTGCAATAAAAGCTAGCTCTTGAAGTAAGAGACGCAGGATTCTGGGCTTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCATATTGCGGTAGGTATTATTCg                                                                                  <  2:977486/138‑1 (MQ=255)
                                                            ctGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCg                                                                                  <  1:775300/142‑1 (MQ=255)
                                                               cGGTTCTGGTAAAAAATACAAGCCGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGCCAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCaa                                                                               <  2:520414/143‑1 (MQ=255)
                                                               cGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCaa                                                                               <  1:346003/143‑1 (MQ=255)
                                                               cGGTGCTGGTAAAAAAGACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCaa                                                                               <  2:622116/143‑1 (MQ=255)
                                                                ggTTGTGGTAAAAAATAGAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGTGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAAc                                                                              <  1:986077/143‑1 (MQ=255)
                                                                  ttcctgTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACga                                                                            <  2:116979/138‑1 (MQ=255)
                                                                  ttCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACga                                                                            <  2:530900/143‑1 (MQ=255)
                                                                     tGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGaa                                                                         <  1:913964/143‑1 (MQ=255)
                                                                       gTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTCGGTATTATTCGCAACGAGGACa                                                                       >  2:1067313/1‑143 (MQ=255)
                                                                       gTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGACTCTGCTCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACa                                                                       >  1:695494/1‑143 (MQ=255)
                                                                       gTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGCTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGATAACa                                                                       >  2:259517/1‑143 (MQ=255)
                                                                        tAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACaa                                                                      >  1:780145/1‑143 (MQ=255)
                                                                          aaaaaTACAAGCAGTGCCATTGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTGTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATg                                                                    <  2:1019324/143‑1 (MQ=255)
                                                                                  aaGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                            >  1:661181/1‑143 (MQ=255)
                                                                                  aaGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                            <  1:684646/143‑1 (MQ=255)
                                                                                  aaGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                            <  2:163982/143‑1 (MQ=255)
                                                                                  aaGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                            <  2:1050146/143‑1 (MQ=255)
                                                                                     cggTGCCATGGCCGCCTGCAATAAAAGCTGACTGGTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTGATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtata                                                         <  2:106089/141‑1 (MQ=255)
                                                                                      aGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATaa                                                        >  1:593535/1‑143 (MQ=255)
                                                                                          ccATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATGTTTATGAcacg                                                    >  2:689962/1‑143 (MQ=255)
                                                                                              gtccgtcTGCAATAAAAGTTAACTGTGTAAGTAAAAGGCTCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTTCTATTCGCGTTTAGAACAATGAAATCTTTATAACGCGTcgc                                                <  2:341842/141‑1 (MQ=255)
                                                                                                ccgccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGCGAACAATGAAATCTTAATAACGCGTCgcgc                                              >  1:1090023/1‑143 (MQ=255)
                                                                                                  gcgtGCAATCAAAGCTAACTGATGAAGTAAGAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCgcag                                            <  2:176025/140‑1 (MQ=255)
                                                                                                      gCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTCTTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAACAGa                                        >  2:1072785/1‑143 (MQ=255)
                                                                                                       cAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAGGCTGCAAATTGCGGTAGGTAGAATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGAt                                       >  2:783645/1‑143 (MQ=255)
                                                                                                        aaTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATg                                      <  1:305722/143‑1 (MQ=255)
                                                                                                         aTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATgc                                     <  1:534036/143‑1 (MQ=255)
                                                                                                           aaaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATgcgc                                   <  1:216485/143‑1 (MQ=255)
                                                                                                                 tAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCCAGTTGCGGTGGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATgg                             <  1:755327/143‑1 (MQ=255)
                                                                                                                       ttGAACTAAAAGGCGCAGGATTCTGCGGCTTTAGTATAGGTTTAAGACAATGAAGAAGCTGCAAATTACAGTAGGTATTATTCGCCACGAGAACAATAACATCTTAACAACACGCCGCGCAGCAGATGCGAACAAGGCGaaca                       >  2:694457/1‑141 (MQ=255)
                                                                                                                        tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATaa                      >  1:360252/1‑143 (MQ=255)
                                                                                                                        ctataaaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATaa                      >  1:443343/7‑143 (MQ=255)
                                                                                                                         gAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAAGCTGTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATaaa                     >  1:300445/1‑143 (MQ=255)
                                                                                                                          aaGTAAACGCCGCAGTATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAAc                    <  2:378364/143‑1 (MQ=255)
                                                                                                                          aaGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACga                                                                            >  1:1058990/1‑87 (MQ=255)
                                                                                                                          aaGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAAc                    <  2:375865/143‑1 (MQ=255)
                                                                                                                          aaGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGCTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACga                                                                            <  2:1058990/87‑1 (MQ=255)
                                                                                                                              aaatGGCGAAGGCGTCTGCGACTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                <  2:589507/143‑1 (MQ=255)
                                                                                                                              aaaaGGCGGAGGATTCTGCGGCTTTTTTATAGGGTTAAGACCATGAAAGAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                <  2:36126/143‑1 (MQ=255)
                                                                                                                              aaaaGGCGCACGATTCTGCGCCTTTTTTACAGGTTTAAGACAAAGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                <  1:802920/143‑1 (MQ=255)
                                                                                                                               aaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTGGCAACGAGAACAATGAAATCTTGATAGCGCGTCGCGCAGCAGATGGGCACCGGGCGAATGACCTGGAg               >  2:58362/1‑143 (MQ=255)
                                                                                                                               aaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGGGCAGCAGATGCGCACATGGCGAATAAACTGGAg               >  1:627367/1‑143 (MQ=255)
                                                                                                                                   gcgcAGGATCCTGCGCCCTTTTCATAGGTCTAAGACAATGAAAAAGCCGCACATTGCGGTAGGTATTACTCGCAACGAGAACAATGATATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTc           <  1:446383/143‑1 (MQ=255)
                                                                                                                                     gcAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATCAAATCTTTATAACGCGTCGCGCAGCCGCCGCGCACAAGGCGAATACACTCGACTTcccc         >  2:116292/1‑143 (MQ=255)
                                                                                                                                      cAGGAGTCTGCGCCTTTTTGATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcg        <  1:992170/143‑1 (MQ=255)
                                                                                                                                         gATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcggcg     >  2:684778/1‑143 (MQ=255)
                                                                                                                                         gATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcggcg     >  1:1089773/1‑143 (MQ=255)
                                                                                                                                            tCTGCGCCTTTTTTATAGGTTTAACACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAGGCGTCGCGCAGCAGATGCGCACATGGCGAATCAAACGGCGTTTCCcggccgac  >  2:802323/1‑139 (MQ=255)
                                                                                                                                            |                                                                                                                                              
CGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTA  >  NC_000913/110876‑111158

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