Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 970,240 T→G 6.2% F29C (TTC→TGC)  tgt → queuine tRNA‑ribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947970,2400TG 6.2% 169.7 / 4.0 65F29C (TTC→TGC) tgtqueuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (26/35);  new base G (2/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

TTCATCACCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGCAC  >  NC_002947/970093‑970385
                                                                                                                                                   |                                                                                                                                                 
ttCATCACCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCa                                                                                                                                                 >  2:1033734/1‑150 (MQ=255)
     caCCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCg                                                                                                                                            <  2:602273/150‑1 (MQ=255)
      aCCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCgg                                                                                                                                           <  1:1109044/150‑1 (MQ=255)
       cccGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATGAATGCCggg                                                                                                                                          <  2:658776/150‑2 (MQ=255)
        ccGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTg                                                                                                                                         >  2:217724/1‑150 (MQ=255)
        ccGCAATCCGGCGCCACGCGGGCCCGAGGAACCACTTTGGGGCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTg                                                                                                                                         <  2:1066488/150‑1 (MQ=255)
         cGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTgg                                                                                                                                        <  1:1091945/150‑1 (MQ=255)
         cGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTgg                                                                                                                                        >  2:18858/1‑150 (MQ=255)
         cGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTACTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTgg                                                                                                                                        >  1:443589/1‑150 (MQ=255)
         cGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATCAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCGTGCCGGTgg                                                                                                                                        >  1:197535/1‑150 (MQ=255)
          gCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTggg                                                                                                                                       >  1:638217/1‑150 (MQ=255)
           cAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTgggg                                                                                                                                      <  1:1033734/150‑1 (MQ=255)
             aTCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGAc                                                                                                                                    <  1:133877/150‑1 (MQ=255)
                      cACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCAc                                                                                                                           <  1:41341/150‑1 (MQ=255)
                            ggCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCCCTGTcca                                                                                                                     <  2:638217/150‑3 (MQ=255)
                             gCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAg                                                                                                                    <  1:901386/150‑1 (MQ=255)
                              cccGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAgg                                                                                                                   <  1:217724/150‑1 (MQ=255)
                                cGAGTATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGc                                                                                                                 <  1:449267/150‑1 (MQ=255)
                                     aTCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCAtgct                                                                                                            >  1:784673/1‑150 (MQ=255)
                                          gCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCAc                                                                                                       <  1:1260941/150‑1 (MQ=255)
                                              gAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGa                                                                                                   >  2:1177478/1‑150 (MQ=255)
                                                   gCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCg                                                                                              >  2:189815/1‑150 (MQ=255)
                                                     aCCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAg                                                                                            <  1:18858/150‑1 (MQ=255)
                                                       cTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGc                                                                                          >  1:344039/1‑150 (MQ=255)
                                                          tCGAATGTCCTTCGAACTGCTGGCCGCCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCAt                                                                                       >  1:29140/1‑150 (MQ=255)
                                                                 tCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCAcgcg                                                                                                    <  2:936779/130‑1 (MQ=255)
                                                                 tCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCAcgcg                                                                                                    >  1:936779/1‑130 (MQ=255)
                                                                  ccTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAgg                                                                                           >  2:1245895/1‑138 (MQ=255)
                                                                  ccTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAgg                                                                                           <  1:1245895/138‑1 (MQ=255)
                                                                   cTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCga                                                                              >  1:734497/1‑150 (MQ=255)
                                                                    ttCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCAc                                                                                                                           <  2:894553/104‑1 (MQ=255)
                                                                    ttCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCAc                                                                                                                           >  1:894553/1‑104 (MQ=255)
                                                                    ttCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCgag                                                                             >  2:592649/1‑150 (MQ=255)
                                                                       gAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAgatg                                                                          >  2:1071008/1‑150 (MQ=255)
                                                                         aCTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAgatgat                                                                        >  1:929499/1‑150 (MQ=255)
                                                                              tgGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCg                                                                   <  2:344039/150‑1 (MQ=255)
                                                                              tgGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCg                                                                   <  1:1071008/150‑1 (MQ=255)
                                                                                  cACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCaa                                                               >  2:801049/1‑150 (MQ=255)
                                                                                     cGACGGCAAGGCCCGTCGGGGCCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAacac                                                            <  1:801049/150‑1 (MQ=255)
                                                                                          gCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAaca                                                             <  1:661329/144‑1 (MQ=255)
                                                                                          gCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAaca                                                             >  2:661329/1‑144 (MQ=255)
                                                                                            aaGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCAc                                                     >  2:1116427/1‑150 (MQ=255)
                                                                                                  cGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCga                                                                              <  1:530360/119‑1 (MQ=255)
                                                                                                  cGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCga                                                                              >  2:530360/1‑119 (MQ=255)
                                                                                                   gTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCgg                                                                  <  2:1181008/130‑1 (MQ=255)
                                                                                                   gTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCgg                                                                  >  1:1181008/1‑130 (MQ=255)
                                                                                                         gTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCTATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTc                                        <  2:784673/150‑1 (MQ=255)
                                                                                                           cGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCg                                      <  1:426864/150‑1 (MQ=255)
                                                                                                           cGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCg                                      <  1:1116427/150‑1 (MQ=255)
                                                                                                            gCATCACCTTCCCCCGTGGCGCCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCgg                                     <  2:197535/150‑1 (MQ=255)
                                                                                                            gCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCGGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCgg                                     <  2:929499/150‑1 (MQ=255)
                                                                                                            gCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCgg                                     >  2:752265/1‑150 (MQ=255)
                                                                                                            gCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCgg                                     >  2:1027618/1‑150 (MQ=255)
                                                                                                                    ttCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAg                             <  2:989867/150‑1 (MQ=255)
                                                                                                                     tcccgggTGGCACCGGGGAACCACCGGTATTCATGCAGTTGTGGACCGATGGCACGGTCATGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAgg                            <  2:209884/144‑1 (MQ=255)
                                                                                                                     tCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACAt                                                                                                <  2:554867/82‑1 (MQ=255)
                                                                                                                     tCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACAt                                                                                                >  1:554867/1‑82 (MQ=255)
                                                                                                                     tCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAgg                            <  1:189815/150‑1 (MQ=255)
                                                                                                                         aCGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGAt                        <  2:29140/150‑1 (MQ=255)
                                                                                                                                     ggAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACaa            <  2:562554/150‑1 (MQ=255)
                                                                                                                                         ccaccGGCATGAATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCCCCCGCGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACGATCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGc        >  1:554272/4‑150 (MQ=255)
                                                                                                                                          caccGGCATGAATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTc                                                                    >  2:767301/3‑89 (MQ=255)
                                                                                                                                          caccGGCATGAATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTc                                                                    <  1:767301/87‑1 (MQ=255)
                                                                                                                                             cGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGc    <  1:1177478/150‑1 (MQ=255)
                                                                                                                                               gCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGCAc  <  1:752265/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTCATCACCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGCAC  >  NC_002947/970093‑970385

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