Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,307,453 C→T 7.8% intergenic (‑154/+32) 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,4530CT 7.8% 219.1 / 6.6 77intergenic (‑154/+32)recC/PP_5SGmultifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (39/32);  new base T (3/3);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

AGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGAT  >  NC_002947/5307304‑5307571
                                                                                                                                                     |                                                                                                                      
aGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTAc                                                                                                                        <  1:2399373/150‑1 (MQ=255)
aGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTAc                                                                                                                        <  1:182434/150‑1 (MQ=255)
  gTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCt                                                                                                                      >  1:1962913/1‑150 (MQ=255)
  gTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCACGGCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCt                                                                                                                      >  1:1935193/1‑150 (MQ=255)
    ccttGGAGGCGCTGATCGCTGAGCGGGGCACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGCAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCATGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCtgc                                                                                                                    <  2:378187/147‑1 (MQ=39)
     tctGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCtgcg                                                                                                                   >  1:1009909/1‑150 (MQ=255)
       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                                                                                                       <  2:1056722/147‑1 (MQ=255)
                    ggCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAgg                                                                                                       >  1:1056722/1‑147 (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                                                                                           <  1:2313158/150‑1 (MQ=255)
                             gggAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaaa                                                                                           <  2:1009909/150‑1 (MQ=255)
                             gagAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaaa                                                                                           <  2:1363247/148‑1 (MQ=255)
                              ggAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGTGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAAt                                                                                          <  2:1374377/150‑1 (MQ=255)
                               gAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaa                                                                                            <  1:2132559/147‑1 (MQ=255)
                               gAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaa                                                                                            >  2:2132559/1‑147 (MQ=255)
                                      atCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAggg                                                                                                      >  2:1286984/1‑130 (MQ=255)
                                      atCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAggg                                                                                                      <  1:1286984/130‑1 (MQ=255)
                                      atCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTg                                                                                  <  1:1908285/150‑1 (MQ=255)
                                       ccACGCGTACATGAAACGGCCGAATGCCACAAGGCTGCGGAAGCGGCCTTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGa                                                                                 <  1:1096097/149‑1 (MQ=255)
                                        cACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAggg                                                                                                      <  1:445927/128‑1 (MQ=255)
                                        cACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAggg                                                                                                      >  2:445927/1‑128 (MQ=255)
                                         aCGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACg                                                                               >  2:1809178/1‑150 (MQ=255)
                                         aCGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACg                                                                               >  1:1726395/1‑150 (MQ=255)
                                                 aTGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCtgc                                                                                                                    >  1:2164716/1‑105 (MQ=255)
                                                 aTGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCtgc                                                                                                                    <  2:2164716/105‑1 (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                                                                          >  1:1503543/1‑124 (MQ=255)
                                                                        gCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGAc                                                                          <  2:1503543/124‑1 (MQ=255)
                                                                         cTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCacac                                               <  2:1962913/150‑1 (MQ=255)
                                                                           gCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTa                                             <  2:2179790/150‑1 (MQ=255)
                                                                                   ggCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc                                     <  1:305962/150‑1 (MQ=255)
                                                                                   ggCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc                                     >  1:1426063/1‑150 (MQ=255)
                                                                                    gCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAAt                                                                                          >  1:700894/1‑96 (MQ=255)
                                                                                    gCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAAt                                                                                          <  2:700894/96‑1 (MQ=255)
                                                                                     ccGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGa                                   >  1:441463/1‑150 (MQ=255)
                                                                                     ccGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGa                                   <  1:2158693/150‑1 (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                                                           <  1:1145711/110‑1 (MQ=255)
                                                                                                     aGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATggg                                                           >  2:1145711/1‑110 (MQ=255)
                                                                                                      gCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactact                  >  2:1664020/1‑150 (MQ=255)
                                                                                                      gCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactact                  <  1:1664020/150‑1 (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)
                                                                                                                                                     |                                                                                                                      
AGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGAT  >  NC_002947/5307304‑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: