Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,307,460 G→A 8.2% intergenic (‑161/+25) recC ← / ← PP_5SG multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit/5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,307,4600GA 8.2% 198.6 / 5.9 70intergenic (‑161/+25)recC/PP_5SGmultifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (34/30);  new base A (3/3);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGAT  >  NC_002947/5307311‑5307571
                                                                                                                                                     |                                                                                                               
tGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCtgcgtg                                                                                                                 >  2:2313158/1‑150 (MQ=255)
    ggCGCGGATGGCTGGGTGGGGAACGATATCACGCGTACATGAATCGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCa                                                                                                             <  2:423427/150‑1 (MQ=255)
      cgcTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGa                                                                                                           <  1:802911/150‑1 (MQ=255)
         tGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAg                                                                                                        >  2:462896/1‑150 (MQ=255)
             ggCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAgg                                                                                                       >  1:1056722/1‑147 (MQ=255)
             ggCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAgg                                                                                                       <  2:1056722/147‑1 (MQ=255)
                tGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGGCCTCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAGGGGtttt                                                                                                 <  2:686339/150‑1 (MQ=38)
                  agagGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGc                                                                                               >  2:1386071/1‑150 (MQ=255)
                  agagGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGc                                                                                               >  1:2179790/1‑150 (MQ=255)
                   gagGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAGGGGTTTTgtc                                                                                              >  2:1103866/1‑148 (MQ=39)
                      gggAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaaa                                                                                           <  2:1009909/150‑1 (MQ=255)
                      gggAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaaa                                                                                           <  1:2313158/150‑1 (MQ=255)
                      gagAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaaa                                                                                           <  2:1363247/148‑1 (MQ=255)
                       ggAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGTGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAAt                                                                                          <  2:1374377/150‑1 (MQ=255)
                        gAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaa                                                                                            >  2:2132559/1‑147 (MQ=255)
                        gAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaa                                                                                            <  1:2132559/147‑1 (MQ=255)
                               atCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAggg                                                                                                      <  1:1286984/130‑1 (MQ=255)
                               atCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAggg                                                                                                      >  2:1286984/1‑130 (MQ=255)
                               atCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTg                                                                                  <  1:1908285/150‑1 (MQ=255)
                                ccACGCGTACATGAAACGGCCGAATGCCACAAGGCTGCGGAAGCGGCCTTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGa                                                                                 <  1:1096097/149‑1 (MQ=255)
                                 cACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAggg                                                                                                      >  2:445927/1‑128 (MQ=255)
                                 cACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAggg                                                                                                      <  1:445927/128‑1 (MQ=255)
                                  aCGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACg                                                                               >  1:1726395/1‑150 (MQ=255)
                                  aCGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACg                                                                               >  2:1809178/1‑150 (MQ=255)
                                          aTGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTa                                                                       >  2:1211447/1‑150 (MQ=255)
                                          aTGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTa                                                                       >  1:1524335/1‑150 (MQ=255)
                                            gAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTAct                                                                     <  2:1383026/150‑1 (MQ=255)
                                                     gAATGTCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgg                                                            <  2:711819/150‑1 (MQ=255)
                                                     gAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGACTCTTGACGATGACCTACTCTCACATgg                                                            <  1:2234779/150‑1 (MQ=255)
                                                     gAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgg                                                            >  1:136787/1‑150 (MQ=255)
                                                      aaTGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATggg                                                           <  2:1726395/150‑1 (MQ=255)
                                                           cATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAAc                                                      >  1:721331/1‑150 (MQ=255)
                                                            aTCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAAcc                                                     >  2:2158693/1‑150 (MQ=255)
                                                             tCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAAccc                                                    <  2:241539/150‑1 (MQ=255)
                                                                 gCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGAc                                                                          <  2:1503543/124‑1 (MQ=255)
                                                                 gCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGAc                                                                          >  1:1503543/1‑124 (MQ=255)
                                                                  cTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCacac                                               <  2:1962913/150‑1 (MQ=255)
                                                                    gCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTa                                             <  2:2179790/150‑1 (MQ=255)
                                                                            ggCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc                                     >  1:1426063/1‑150 (MQ=255)
                                                                            ggCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc                                     <  1:305962/150‑1 (MQ=255)
                                                                             gCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAAt                                                                                          <  2:700894/96‑1 (MQ=255)
                                                                             gCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAAt                                                                                          >  1:700894/1‑96 (MQ=255)
                                                                              ccGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGa                                   <  1:2158693/150‑1 (MQ=255)
                                                                              ccGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGa                                   >  1:441463/1‑150 (MQ=255)
                                                                                ctGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTAATGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATg                                 <  1:508379/149‑1 (MQ=255)
                                                                                 tGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGc                                <  1:462896/150‑1 (MQ=255)
                                                                                  gCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTGTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCa                               >  1:1697310/1‑150 (MQ=255)
                                                                                    ccGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATc                             >  1:249208/1‑150 (MQ=255)
                                                                                    ccGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATc                             >  2:893701/1‑150 (MQ=255)
                                                                                     cGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCg                            >  2:2100921/1‑150 (MQ=255)
                                                                                      gCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGt                           >  2:980675/1‑150 (MQ=255)
                                                                                      gCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGt                           >  2:1456331/1‑150 (MQ=255)
                                                                                          ccGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCa                       <  2:721331/150‑1 (MQ=255)
                                                                                            gCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATTAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCact                     >  1:2307153/1‑150 (MQ=255)
                                                                                            gCAGCGATTCACGCAGGTCGAGTTTTTTGCGTACAACGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCact                     >  2:758319/1‑150 (MQ=255)
                                                                                              aGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATggg                                                           >  2:1145711/1‑110 (MQ=255)
                                                                                              aGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATggg                                                           <  1:1145711/110‑1 (MQ=255)
                                                                                               gCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactact                  <  1:1664020/150‑1 (MQ=255)
                                                                                               gCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactact                  >  2:1664020/1‑150 (MQ=255)
                                                                                                   ttCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGt              >  1:767148/1‑150 (MQ=39)
                                                                                                      gcgcAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAggg                                                                                                      >  2:1262000/1‑59 (MQ=17)
                                                                                                      gcgcAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAggg                                                                                                      <  1:1262000/59‑1 (MQ=17)
                                                                                                        gcAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCggg         >  1:1419022/1‑150 (MQ=35)
                                                                                                            gTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGCCCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatgg     <  1:1053201/150‑1 (MQ=21)
                                                                                                             tCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactact                  <  1:1929105/136‑1 (MQ=33)
                                                                                                             tCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactact                  >  2:1929105/1‑136 (MQ=33)
                                                                                                             tCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatggg    >  1:168856/1‑150 (MQ=33)
                                                                                                             tCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatggg    >  1:1762069/1‑150 (MQ=33)
                                                                                                             tCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatggg    <  2:1524335/150‑1 (MQ=33)
                                                                                                               gAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatgggat  >  1:1681879/1‑150 (MQ=32)
                                                                                                               gAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCgggatgggat  >  1:1818325/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                               
TGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGAT  >  NC_002947/5307311‑5307571

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