Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,013 G→A 6.5% intergenic (+29/‑31) 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,0130GA 6.5% 189.7 / 5.4 78intergenic (+29/‑31)secA/mutTprotein translocation ATPase/8‑oxo‑dGTP diphosphatase
Reads supporting (aligned to +/- strand):  ref base G (35/37);  new base A (2/3);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-02

CGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCG  >  NC_000913/110876‑111155
                                                                                                                                         |                                                                                                                                              
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)
                                                                                                                                         |                                                                                                                                              
CGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCG  >  NC_000913/110876‑111155

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: