Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,457,903 G→A 6.6% L545L (CTG→CTA leuS → leucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,457,9030GA 6.6% 177.9 / 2.6 60L545L (CTG→CTAleuSleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base A (2/2);  total (28/32)
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.25e-01

TGGACACCTTCGTCGAGTCGTCCTGGTACTTCGCCCGCTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAACGGCGGCAAGGACTGGTTCAACC  >  NC_002947/5457758‑5458051
                                                                                                                                                   |                                                                                                                                                     
tGGACACCTTCGTCGAGTCGTCCTGGTACTTCGCCCGCTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACctgctg                                                                                                                                                    >  3:218958/1‑149 (MQ=255)
       cTTCGTCGAGTCGTCCTGGTACTTCGCCCGCTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTAcgcgc                                                                                                                                             >  1:297906/1‑149 (MQ=255)
       cTTCGTCGAGTCGTCCTGGTACTTCGCCCGCTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTACTGTAcgcgc                                                                                                                                             >  6:29286/1‑149 (MQ=255)
              gAGTCGTCCTGGTACTTCGCCCGCTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGcttctt                                                                                                                                      >  7:96627/1‑149 (MQ=255)
               aGTCGTCCTGGTACTTCGCCCGCTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCttcttc                                                                                                                                     <  8:64587/149‑1 (MQ=255)
                         gTACTTCGCCCGCTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGa                                                                                                                           <  5:248820/149‑1 (MQ=255)
                         gTACTTCGCCCGCTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGa                                                                                                                           >  2:24464/1‑149 (MQ=255)
                               cGCCCGCTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGt                                                                                                                      <  2:297906/148‑1 (MQ=255)
                                  ccGCTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTACTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACg                                                                                                                  <  5:29286/149‑1 (MQ=255)
                                     cTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAg                                                                                                               <  4:218958/149‑1 (MQ=255)
                                        cGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGc                                                                                                             <  1:24464/148‑1 (MQ=255)
                                         gCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCt                                                                                                           >  2:209134/1‑149 (MQ=255)
                                          ccTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTg                                                                                                          >  3:2245/1‑149 (MQ=255)
                                              gacAAACTACGACAAAGGCCTTGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCTGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTc                                                                                                       <  5:310551/146‑1 (MQ=255)
                                              gCCAAACTACGACAAGGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTg                                                                                                          >  1:39535/1‑145 (MQ=255)
                                              gCCAAACTACGACAAGGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTg                                                                                                          <  2:39535/145‑1 (MQ=255)
                                                cAAACTACGACAAAGGCCTGGTCGACCCCAAATCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCATCCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGc                                                                                                                        <  1:207456/130‑1 (MQ=255)
                                                cAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGc                                                                                                                        >  2:207456/1‑129 (MQ=255)
                                                cAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCAcc                                                                                                    >  2:91585/1‑149 (MQ=255)
                                                  aaCTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTc                                                                                                  >  4:135620/1‑149 (MQ=255)
                                                     tACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAGACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCg                                                                                                 >  1:202953/1‑148 (MQ=255)
                                                     tACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGaa                                                                                               >  4:248000/1‑149 (MQ=255)
                                                        gACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAAc                                                                                              >  5:243407/1‑147 (MQ=255)
                                                          cAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGc                                                                                          >  3:109169/1‑149 (MQ=255)
                                                               gCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTc                                                                                    <  4:2245/149‑1 (MQ=255)
                                                                 cTGGTCGACACCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTc                                                                                    >  5:95494/1‑147 (MQ=255)
                                                                        aCCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACctg                                                                           <  4:266285/149‑1 (MQ=255)
                                                                           ccAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTc                                                                                    <  4:135245/137‑1 (MQ=255)
                                                                           ccAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTc                                                                                    >  3:135245/1‑137 (MQ=255)
                                                                            cAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCt                                                                                                   >  3:286253/1‑122 (MQ=255)
                                                                            cAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCt                                                                                                   <  4:286253/122‑1 (MQ=255)
                                                                            cAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGa                                                                       <  3:135620/149‑1 (MQ=255)
                                                                                  cgccAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAgg                                                                 >  5:309319/1‑149 (MQ=255)
                                                                                   gccAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGAGCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAggg                                                                <  2:228381/149‑1 (MQ=255)
                                                                                     cAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGAccc                                                                    <  6:126875/143‑1 (MQ=255)
                                                                                     cAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGAccc                                                                    >  5:126875/1‑143 (MQ=255)
                                                                                         ccACTGGCTGCCGGTAGATCATTACATCGGCGGTATCGAACACGCGATACTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATTCGTGACGAAGTCCTGGTCACCTCGAACGAGCAGTTCAAGAACCTGCTGACCCAGGGCAttg                                                           <  6:158810/149‑3 (MQ=255)
                                                                                          cACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTc                                                                                                       <  6:105769/105‑1 (MQ=255)
                                                                                          cACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTc                                                                                                       >  5:105769/1‑105 (MQ=255)
                                                                                           aCTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCAtggtg                                                         <  7:266799/149‑1 (MQ=255)
                                                                                               gCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTc                                                      <  8:96627/148‑1 (MQ=255)
                                                                                                     ggTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGaaa                                               <  6:309319/149‑1 (MQ=255)
                                                                                                      gTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAgg                                                                 >  5:21597/1‑130 (MQ=255)
                                                                                                      gTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAgg                                                                 <  6:21597/130‑1 (MQ=255)
                                                                                                      gTACATCACTACCTCGGCGGTAGCGAACACGCGATCCAGCACCTGCTGGACGCGCCCTTCTTCCACAAGCTGCTGCGGGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAAc                                              <  4:230531/149‑1 (MQ=255)
                                                                                                         gATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACcta                                           <  1:209134/149‑1 (MQ=255)
                                                                                                                 cATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACcgcg                                   <  6:95494/149‑1 (MQ=255)
                                                                                                                 cATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACcgc                                    <  4:109169/148‑1 (MQ=255)
                                                                                                                  aTCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCAcc                                                                                                    <  8:49237/84‑1 (MQ=255)
                                                                                                                  aTCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCAcc                                                                                                    >  7:49237/1‑84 (MQ=255)
                                                                                                                     ggcggTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGc                               >  3:226491/1‑149 (MQ=255)
                                                                                                                      gcggTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGc                               >  1:56125/1‑148 (MQ=255)
                                                                                                                              gAACACGCGATCCTGCACCTACTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAAcgg                      >  3:116291/1‑149 (MQ=255)
                                                                                                                               aaCACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAAcggc                     >  2:81366/1‑149 (MQ=255)
                                                                                                                                acacGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAAcggcg                    <  2:202953/149‑1 (MQ=255)
                                                                                                                                 cacGCGATCCTGCACCTACTGTACGCGCGCTTCTTCCACAAGCTGATGCGTTACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAAcggcgg                   <  4:116291/149‑1 (MQ=255)
                                                                                                                                   cgcgATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACctac                                          >  7:230644/1‑124 (MQ=255)
                                                                                                                                   cgcgATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACctac                                          <  8:230644/124‑1 (MQ=255)
                                                                                                                                     cgATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAACGGCGGCAAg               <  3:248000/149‑1 (MQ=255)
                                                                                                                                           tGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAACGGCGGCAAGGACTgg         >  3:196543/1‑149 (MQ=255)
                                                                                                                                             cACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAACGGCGGCAAGGACTGGtt       <  4:226491/149‑1 (MQ=255)
                                                                                                                                                 ctgctgTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAACGGCGGCAAGGACTGGTTCaa    >  7:253844/1‑149 (MQ=255)
                                                                                                                                                  tgctgTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCAAGTTCAAGAACATGATGACCCAGGGCTGGGTGTTCGCCGAACCCTACTACCGCGTTGCCAGCAACGCCGGCAAGGACTGGTTCaa    >  1:252617/1‑149 (MQ=255)
                                                                                                                                                   gctgTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAACGGCGGCAAGGACTGGTTCAAcc  >  5:58958/1‑149 (MQ=255)
                                                                                                                                                   gctgTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAACGGCGGCAAGGACTGGTTCAAcc  >  1:69530/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGGACACCTTCGTCGAGTCGTCCTGGTACTTCGCCCGCTACGCCTCGCCAAACTACGACAAAGGCCTGGTCGACCCCAAAGCCGCCAACCACTGGCTGCCGGTAGATCAGTACATCGGCGGTATCGAACACGCGATCCTGCACCTGCTGTACGCGCGCTTCTTCCACAAGCTGATGCGTGACGAAGGCCTGGTCACCTCGAACGAGCCGTTCAAGAACCTGCTGACCCAGGGCATGGTGGTCGCCGAAACCTACTACCGCGTTGCCAGCAACGGCGGCAAGGACTGGTTCAACC  >  NC_002947/5457758‑5458051

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