Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 524,368 C→T 6.6% intergenic (+232/‑584) tyrS → / → PP_16SC tyrosine‑‑tRNA ligase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947524,3680CT 6.6% 168.2 / 2.8 61intergenic (+232/‑584)tyrS/PP_16SCtyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (31/25);  new base T (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

GTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAA  >  NC_002947/524219‑524514
                                                                                                                                                      |                                                                                                                                                  
gtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtc                                                                                                                                                    <  1:1682328/150‑1 (MQ=255)
      ggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAg                                                                                                                                              <  1:249046/150‑1 (MQ=25)
       gagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAgg                                                                                                                                             <  1:1452477/150‑1 (MQ=255)
          ggggCTGCTTTGAGAGGGTGTTGGTTGTGAGATTTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATACGCGGCCTATAAAACGGTCCGCCTCGCGGGCGGCGCTCGGTCTCATAGgcgc                                                                                                                                          >  2:1729043/1‑150 (MQ=21)
           gggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCt                                                                                                                                         <  1:1729043/150‑1 (MQ=255)
            gtcTGTTTTGAGGGGGTGTTGGCGGTCAGATGCGTCCGCGACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCCCCGCGCGGGCGGCCCTCGGTCTCACAGgcgcgg                                                                                                                                        >  2:41171/3‑148 (MQ=25)
            ggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGgcg                                                                                                                                           >  1:667835/1‑147 (MQ=255)
            ggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGgcg                                                                                                                                           <  2:667835/147‑1 (MQ=255)
               tGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGaaa                                                                                                                                     >  2:1585390/1‑150 (MQ=255)
                          ggTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCg                                                                                                                          >  2:2265061/1‑150 (MQ=255)
                          ggTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGAGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCg                                                                                                                          >  2:198987/1‑150 (MQ=255)
                           gtgtTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGa                                                                                                                         >  1:1973670/1‑150 (MQ=255)
                              ttGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAaca                                                                                                                      >  1:345127/1‑150 (MQ=255)
                                 gCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACtt                                                                                                                   <  1:1724781/150‑1 (MQ=255)
                                   ggTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                                                 <  2:67586/150‑1 (MQ=255)
                                      gagaTGTGTCCGCCACAACCCTTGCAGCGCCTGTGATATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAg                                                                                                              >  2:477077/1‑150 (MQ=255)
                                              tCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATa                                                                                                      >  2:1360336/1‑150 (MQ=255)
                                                 gCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCACATCTGAAGCGCTGCCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGACACTGTCGAGGCGAACACTTGGCAGCCTCGATTCgg                                                                                                   >  2:1791813/1‑150 (MQ=21)
                                                         ccTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTc                                                                                                                                <  2:352987/113‑1 (MQ=255)
                                                         ccTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTc                                                                                                                                >  1:352987/1‑113 (MQ=255)
                                                            tGCAGCGCCTGTGAGATCGCGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCt                                                                                        <  2:1163888/150‑1 (MQ=255)
                                                            tGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCt                                                                                        <  1:3866/150‑1 (MQ=255)
                                                               aGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTcccc                                                                                         >  2:406443/1‑146 (MQ=255)
                                                               aGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTcccc                                                                                         <  1:406443/146‑1 (MQ=255)
                                                               aGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTtat                                                                                     >  2:815894/1‑150 (MQ=255)
                                                                          gaTCGAGCGCCGCCCGCGTCCGTCTTGAGACATTCACCGCCTATGAAGCCGAGCGCCGCGCGCGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCCCTTTc                                                                          <  1:1479151/150‑1 (MQ=255)
                                                                          gaTCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCCCCGAGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTc                                                                          >  2:934964/1‑150 (MQ=255)
                                                                           aTCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCt                                                                         >  1:2219591/1‑150 (MQ=255)
                                                                                      cccgcGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAg                                                              <  1:2265061/150‑1 (MQ=255)
                                                                                         gcgggcGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAACCTGTACTGCCGAACACTTGGCAGCCTCGATACGGTGTTTTCCCATTCTATGTCACTTGCTGGAAATTAAAGttt                                                           >  1:418977/6‑150 (MQ=18)
                                                                                         gcgTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTAGAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGtt                                                           <  1:1496814/150‑1 (MQ=255)
                                                                                            tCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGAc                                                        <  1:584649/150‑1 (MQ=255)
                                                                                            ggcGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGAc                                                        <  1:1940705/148‑1 (MQ=255)
                                                                                                      agaCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCtt                                                                        <  1:607702/125‑1 (MQ=255)
                                                                                                      agaCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCtt                                                                        >  2:607702/1‑125 (MQ=255)
                                                                                                        aCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGaa                                             <  1:41171/150‑1 (MQ=255)
                                                                                                          aTTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATc                                           >  2:19456/1‑150 (MQ=255)
                                                                                                           ttCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATcc                                          >  1:1332296/1‑150 (MQ=255)
                                                                                                            tCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                         >  1:144072/1‑150 (MQ=255)
                                                                                                            tCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                         <  1:1585390/150‑1 (MQ=255)
                                                                                                              aCCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTtat                                                                                     <  1:168513/104‑1 (MQ=255)
                                                                                                              aCCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTtat                                                                                     >  2:168513/1‑104 (MQ=255)
                                                                                                              aCCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt                                       >  2:1662156/1‑150 (MQ=255)
                                                                                                               gcGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGGTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATa                                                                                                      <  1:2177372/85‑1 (MQ=11)
                                                                                                               gcGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGGTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATa                                                                                                      >  2:2177372/2‑86 (MQ=11)
                                                                                                                   cTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATaa                                  <  1:1360336/150‑1 (MQ=255)
                                                                                                                     aTGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATg                                <  1:815894/150‑1 (MQ=255)
                                                                                                                     aTGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATg                                >  2:359341/1‑150 (MQ=255)
                                                                                                                          aCCGAGCGCCGCCCGCGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt                                       <  2:330392/138‑1 (MQ=255)
                                                                                                                          aCCGAGCGCCGCCCGCGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt                                       >  1:330392/1‑138 (MQ=255)
                                                                                                                           ccGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGcccc                          >  1:772657/1‑150 (MQ=255)
                                                                                                                                cgccgcCCGCGCGGCGCTCGGTTTCATAGGCGGTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCACCCTATAATTCCCCCCACTTc                     >  2:2242950/1‑150 (MQ=11)
                                                                                                                                 gccgcGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTcc                    >  1:1019339/1‑150 (MQ=255)
                                                                                                                                  ccgcGCGGGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCt                                                                         <  1:883119/96‑1 (MQ=0)
                                                                                                                                  ccgcGCGGGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCt                                                                         >  2:883119/1‑96 (MQ=0)
                                                                                                                                    gcgcgGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                                                 >  2:371934/1‑54 (MQ=255)
                                                                                                                                    gcgcgGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                                                 <  1:371934/54‑1 (MQ=255)
                                                                                                                                    gcgcgGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc                 <  2:1973670/150‑1 (MQ=255)
                                                                                                                                    gcgcgGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc                 <  1:359341/150‑1 (MQ=255)
                                                                                                                                          gcggcgCTCGGTCTCATAGGCGCTGAAACTGTCCTGGCGAACACTTGGCAGCCTCGATACGGCCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCTCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATc           >  1:718798/1‑150 (MQ=255)
                                                                                                                                              cgcTCGGTATCATAAGCGCTGAAACTGTCGTGGCTAACACTTGGCAGGCTCCATACGCTCTTTGACCCTTATATTTCACTTTATTGAAATTAAAGGTTGACGCCCTCTCGATTCCACTTACAATGCGCCCCACCTCCAGCGACATCgata       <  2:718798/150‑4 (MQ=14)
                                                                                                                                               gcTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAAc      >  1:1762815/1‑150 (MQ=255)
                                                                                                                                                   gggctcATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATATGGCTATAATGCGCCCCACTTCCAGCGACATCGGAACGaga  >  1:2308695/4‑148 (MQ=14)
                                                                                                                                                     tctcATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGaaa  >  2:1358311/1‑148 (MQ=255)
                                                                                                                                                     tctcATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGaaa  <  1:1358311/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAA  >  NC_002947/524219‑524514

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