Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 989,161 T→A 8.3% F338Y (TTT→TAT) ‡ hisS → histidine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947989,1610TA 8.3% 171.9 / 5.8 60F338Y (TTT→TAT) ‡hisShistidine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (26/29);  new base A (3/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

ACGACGGCCTGGTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAGCCGACAAGAGCGGTGCACTGTTCGC  >  NC_002947/989014‑989305
                                                                                                                                                    |                                                                                                                                                
gtaacgGCCTGGTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCttt                                                                                                                                                 <  5:272538/146‑1 (MQ=255)
 cgacgGCTTGGTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTg                                                                                                                                                <  2:279705/149‑1 (MQ=255)
 cgacgGCCTGGTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCttt                                                                                                                                                 <  5:122563/148‑1 (MQ=255)
   acgGCCTGGTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCttt                                                                                                                                                 >  3:60698/1‑146 (MQ=255)
   acgGCCTGGTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCttt                                                                                                                                                 <  4:60698/146‑1 (MQ=255)
        cTGGTCGAACAGATGGGTGGCAAGCCGACCAGCTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGaa                                                                                                                                           >  1:35876/1‑148 (MQ=255)
          ggTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAgg                                                                                                                                       >  4:279027/1‑149 (MQ=255)
           ggcGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTCGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGc                                                                                                                                      <  3:282659/147‑1 (MQ=255)
           gTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTTCGCCAAAGGCGCACAg                                                                                                                                        <  7:235404/147‑1 (MQ=255)
           gTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGc                                                                                                                                      >  7:30600/1‑149 (MQ=255)
            tCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCg                                                                                                                                     <  2:738/149‑1 (MQ=255)
                  aGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTg                                                                                                                               >  3:126040/1‑149 (MQ=255)
                        gTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTggccggc                                                                                                                         >  7:201104/1‑149 (MQ=255)
                         tGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTggccggcc                                                                                                                        >  3:192871/1‑149 (MQ=255)
                           gCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTg                                                                                                                      >  1:265955/1‑149 (MQ=255)
                             aaGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTgcg                                                                                                                    >  2:122050/1‑149 (MQ=255)
                                ccGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCt                                                                                                                 <  4:126040/149‑1 (MQ=255)
                                             gTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCtt                                                                                                     <  3:279027/149‑1 (MQ=255)
                                                ggTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTgc                                                                                                   <  2:59783/148‑1 (MQ=255)
                                                      gCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGAcc                                                                                             <  8:30600/148‑1 (MQ=255)
                                                        cATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCAGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCt                                                                                          >  2:201987/1‑149 (MQ=255)
                                                            ggCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGcc                                                                                      >  1:264638/1‑149 (MQ=255)
                                                              cATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCgg                                                                                    <  3:109396/149‑1 (MQ=255)
                                                                    gcgcCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTgc                                                                              <  1:132986/149‑1 (MQ=255)
                                                                        cTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCt                                                                          >  7:301970/1‑149 (MQ=255)
                                                                        cTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCt                                                                          >  4:271123/1‑149 (MQ=255)
                                                                        cTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCt                                                                          >  8:115822/1‑149 (MQ=255)
                                                                           aTCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTgg                                                                        <  2:265955/148‑1 (MQ=255)
                                                                            tCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCa                                                                      <  7:115822/149‑1 (MQ=255)
                                                                             ccTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGc                                                                                  <  7:267759/136‑1 (MQ=255)
                                                                             ccTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGc                                                                                  >  8:267759/1‑136 (MQ=255)
                                                                                 ctgTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCaa                                                                 <  4:192871/149‑1 (MQ=255)
                                                                                    ttGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGTTCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGc                                                              <  2:35876/149‑1 (MQ=255)
                                                                                    ttGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGc                                                              <  2:185933/149‑1 (MQ=255)
                                                                                        aaaCCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCggt                                                          >  6:115922/1‑149 (MQ=255)
                                                                                        aaaCCCTGGGCAAGGTGCCAGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGTCGGGCCTGCGCCTGGCAGTCAACGCCggt                                                          <  1:201987/149‑1 (MQ=255)
                                                                                           cccTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCggtggt                                                       <  8:154131/149‑1 (MQ=255)
                                                                                            ccTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCGTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGACGttgttg                                                      >  8:233706/1‑149 (MQ=255)
                                                                                                ggCAAGGAGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGc                                                    <  7:233706/147‑1 (MQ=255)
                                                                                                  cAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCt                                                >  8:197058/1‑149 (MQ=255)
                                                                                                  cAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCt                                                >  4:313695/1‑149 (MQ=255)
                                                                                                  cAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCt                                                <  1:122050/149‑1 (MQ=255)
                                                                                                    aGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTc                                              <  2:264638/149‑1 (MQ=255)
                                                                                                      gTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCGTc                                              >  3:28531/1‑147 (MQ=255)
                                                                                                       tGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCaa                                            >  8:314511/1‑148 (MQ=255)
                                                                                                        gCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCTAAGGCTTGCGCGACCGCTTGCCGTGCCTGCGCCTGGCAGTCAACGCCTGTGGTGGCAGCTTCAaga                                          <  8:201104/149‑1 (MQ=255)
                                                                                                         ccGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAagag                                         <  7:314511/149‑1 (MQ=255)
                                                                                                         ccGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAaga                                          <  3:271123/148‑1 (MQ=255)
                                                                                                             aGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTc                                                                   >  4:318028/1‑119 (MQ=255)
                                                                                                             aGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTc                                                                   <  3:318028/119‑1 (MQ=255)
                                                                                                               tCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGtt                                   >  1:3346/1‑149 (MQ=255)
                                                                                                                  aTCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCaa                                <  2:280526/149‑1 (MQ=255)
                                                                                                                   tCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTc                                  >  4:225673/1‑146 (MQ=255)
                                                                                                                      gccgccAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGaaa                            <  8:301970/149‑1 (MQ=255)
                                                                                                                       ccgccAGATCGCCGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAg                           >  8:291798/1‑149 (MQ=255)
                                                                                                                                 gACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAGCCGACAaga                  <  6:102250/148‑1 (MQ=255)
                                                                                                                                     tccgccTGTTCGCCAAAGGCGCACAGGTGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGtggtgg                                                     >  4:67738/5‑109 (MQ=38)
                                                                                                                                     tccgccTGTTCGCCAAAGGCGCACAGGTGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGtggtgg                                                     <  3:67738/105‑1 (MQ=38)
                                                                                                                                       cgccTGTTCGCCAAAGGCGCACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAGCCGACAAGAGCGGTg            >  4:126791/3‑148 (MQ=255)
                                                                                                                                        aCCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAGCCGACAAGAGCGGTGCa          >  1:283299/1‑149 (MQ=255)
                                                                                                                                        aCCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAGCCGACAAGAGCGGTGCa          >  6:299806/1‑149 (MQ=255)
                                                                                                                                         ccTGTTCGCCAAAGGCGCACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAGCCGACAAGAGCGGTGCAc         >  7:108611/1‑149 (MQ=255)
                                                                                                                                         ccTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAGCCGACAAGAGCGGTGCa          <  7:197058/148‑1 (MQ=255)
                                                                                                                                                gcCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAGCCGACAAGAGCGGTGCACTGTTCGc  <  2:283299/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ACGACGGCCTGGTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAGCCGACAAGAGCGGTGCACTGTTCGC  >  NC_002947/989014‑989305

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: ‑