Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 989,155 G→T 7.9% C336F (TGC→TTC)  hisS → histidine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947989,1550GT 7.9% 198.6 / 6.7 63C336F (TGC→TTC) hisShistidine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base T (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

GCCGTTACGACGGCCTGGTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAGCCGACAAGAGCGGTGCAC  >  NC_002947/989008‑989298
                                                                                                                                                    |                                                                                                                                               
gCCGTTACGACGGCCTGGTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTgc                                                                                                                                                <  5:228538/149‑1 (MQ=255)
      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                                                                   >  8:115822/1‑149 (MQ=255)
                                                                              cTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCt                                                                   >  4:271123/1‑149 (MQ=255)
                                                                                 aTCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTgg                                                                 <  2:265955/148‑1 (MQ=255)
                                                                                  tCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCa                                                               <  7:115822/149‑1 (MQ=255)
                                                                                   ccTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGc                                                                           >  8:267759/1‑136 (MQ=255)
                                                                                   ccTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGc                                                                           <  7:267759/136‑1 (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                                                            <  3:318028/119‑1 (MQ=255)
                                                                                                                   aGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTc                                                            >  4:318028/1‑119 (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)
                                                                                                                                                    |                                                                                                                                               
GCCGTTACGACGGCCTGGTCGAACAGATGGGTGGCAAGCCGACCA‑CTGGCGTTGGTTTCGCCATGGGCATCGAGCGCCTGATCCTGCTGTTGGAAACCCTGGGCAAGGTGCCGGAGTCCATCAGCCGCCAGATCGACGTTTACCTGTGCGCCTTTGGCGAACAGGCGGAGCTGGCCGGCCTGCGCCTGTCCGAAGGCTTGCGCGACCGCTTGCCGGGCCTGCGCCTGGCAGTCAACGCCGGTGGTGGCAGCTTCAAGAGCCAGTTCAAGAAAGCCGACAAGAGCGGTGCAC  >  NC_002947/989008‑989298

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