Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,896,809 T→C 6.3% L43P (CTG→CCG)  PP_2551 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,896,8090TC 6.3% 234.0 / 3.9 79L43P (CTG→CCG) PP_2551LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base C (2/3);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCGCCATCTTCGAC  >  NC_002947/2896665‑2896957
                                                                                                                                                |                                                                                                                                                    
cACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGc                                                                                                                                                   <  1:58451/148‑1 (MQ=255)
cACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCt                                                                                                                                                  >  3:246403/1‑149 (MQ=255)
 aCCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTc                                                                                                                                                 >  3:274311/1‑149 (MQ=255)
   ggCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCTCCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAAGCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCAc                                                                                                                                               <  7:243986/148‑1 (MQ=255)
    gCAGTGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCAGCTGCCGGCTCagg                                                                                                                                              <  8:290810/149‑3 (MQ=255)
    gCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCAcc                                                                                                                                              <  8:253865/149‑1 (MQ=255)
       gggAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGc                                                                                                                                           <  7:290831/149‑1 (MQ=255)
        ggAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCg                                                                                                                                          >  6:24428/1‑149 (MQ=255)
        ggAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCg                                                                                                                                          <  7:40680/149‑1 (MQ=255)
         gAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCTCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGt                                                                                                                                         <  5:43611/149‑1 (MQ=255)
         gAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGt                                                                                                                                         >  6:348108/1‑149 (MQ=255)
         gAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCAGCTGCCGACTCag                                                                                                                                               <  3:82341/143‑2 (MQ=255)
         gAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCAGCTGCCGACTCag                                                                                                                                               >  4:82341/1‑142 (MQ=255)
            ccTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTg                                                                                                                                        >  6:209931/1‑147 (MQ=255)
            ccTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGaa                                                                                                                                      <  6:117443/149‑1 (MQ=255)
                cATGAACCCCACCTTCGCCTCCCTCAGCCTCGGCCACCTGCGCACCCTCGACCACGTGTTGAAGCTTAAAAACCTCAGCCTCGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTTATCCGGCAGCTGGCTCACCTGCGTGAAGCCt                                                                                                                                  <  3:120084/149‑1 (MQ=255)
                 aTGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCtt                                                                                                                                 <  2:320389/149‑1 (MQ=255)
                   gAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTc                                                                                                                                >  7:294665/1‑148 (MQ=255)
                     aCCCCACCTTCGCCTCCCTCAGCCTGGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCATCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTTAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCga                                                                                                                              <  6:126749/148‑1 (MQ=255)
                     aCCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCgat                                                                                                                             >  5:27063/1‑149 (MQ=255)
                     aCCCCACCTTCGCCTACCTCAGCCTGGCCCACCTGCCCAACCTCGACCAACTGTTGCATCTTAAAAACCTCAGACACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCga                                                                                                                              <  6:132595/148‑1 (MQ=255)
                      ccccACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCgat                                                                                                                             >  7:137413/1‑148 (MQ=255)
                        ccACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCgatgat                                                                                                                          <  4:330852/149‑1 (MQ=255)
                            cTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGAGGATCCGt                                                                                                                      >  6:163225/1‑149 (MQ=255)
                             ttCGCCTCCCTCAGCCTAGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCGGCGTGAAGCCTTCGATGATCCGtt                                                                                                                     >  3:54033/1‑149 (MQ=255)
                                gCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCt                                                                                                                  >  3:368497/1‑149 (MQ=255)
                                gCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCAATCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCt                                                                                                                  >  3:368492/1‑149 (MQ=255)
                                     cctcAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTgc                                                                                                             >  3:365126/1‑149 (MQ=255)
                                      ctcAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAg                                                                                                                                     <  1:212334/124‑1 (MQ=255)
                                      ctcAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAg                                                                                                                                     >  2:212334/1‑124 (MQ=255)
                                      ctcAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTgcg                                                                                                            <  5:209931/149‑1 (MQ=255)
                                          gCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCaa                                                                                                        <  8:68345/149‑1 (MQ=255)
                                           ccTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAg                                                                                                       >  3:176446/1‑149 (MQ=255)
                                             tCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAggg                                                                                                     >  6:160236/1‑149 (MQ=255)
                                                cccACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGt                                                                                                  <  4:274311/149‑1 (MQ=255)
                                                 ccgcgcgcgcCCCCGCGACCACCTGGTGCCGCTTACACCCCTCAGCCCCGCGGCCGAGCGCCTGGGGGGCAGCCAATCGGCCCGGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTg                                                                                                 <  3:223133/143‑1 (MQ=255)
                                                 ccACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTg                                                                                                 <  8:294665/149‑1 (MQ=255)
                                                   acttGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGaa                                                                                                                                      >  4:87889/4‑110 (MQ=255)
                                                   acttGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGaa                                                                                                                                      <  3:87889/107‑1 (MQ=255)
                                                         gcACCCTCGACCCCCTGTTGCAGCTTAAAACCCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGCGCCGGCAGCTGGCTCACCTGCGGGAAGCAGTCGAGGATCCTGTGCGGGTACGCCAAGGGCGTg                                                                                                 >  4:223133/1‑141 (MQ=255)
                                                           aCCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCt                                                                                       >  7:135535/1‑149 (MQ=255)
                                                                   ccaccTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAAc                                                                               <  8:137413/149‑1 (MQ=255)
                                                                      ccTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTg                                                                                                 <  1:260295/128‑1 (MQ=255)
                                                                      ccTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTg                                                                                                 >  2:260295/1‑128 (MQ=255)
                                                                      ccTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCAGCGTGAAGACTTCGATGATGCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATg                                                                            >  4:310119/1‑149 (MQ=255)
                                                                       cTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAATCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGc                                                                           <  6:27063/149‑1 (MQ=255)
                                                                            gCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAAc                                                                               >  4:103597/1‑140 (MQ=255)
                                                                            gCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAAc                                                                               <  3:103597/140‑1 (MQ=255)
                                                                            gCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAg                                                                      >  8:44511/1‑149 (MQ=255)
                                                                                cTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCt                                                                  <  8:135535/149‑1 (MQ=255)
                                                                                   aaaaaCCTCAGCCTCGCTGCCGAGCTAATGGGTTTCAGCCAATCGGCCCTGAGCCGGAAGCTGGATCACCTGCGTGAAGCATTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTTGTTATGTGCTTAGCGAACATGCCGAAGCCCTgg                                                                <  4:368492/148‑1 (MQ=255)
                                                                                   aaaaaCCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTgg                                                                <  4:368497/148‑1 (MQ=255)
                                                                                         cTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCAGGTCGAGcc                                                         >  5:364004/1‑149 (MQ=255)
                                                                                            aGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGc                                                                                                                                           >  4:113391/1‑64 (MQ=255)
                                                                                            aGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGc                                                                                                                                           <  3:113391/64‑1 (MQ=255)
                                                                                             gCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTg                                                     <  5:348108/149‑1 (MQ=255)
                                                                                              ccACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGcc                                                                          >  4:5774/1‑127 (MQ=255)
                                                                                              ccACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCGGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTgc                                                    <  5:24428/149‑1 (MQ=255)
                                                                                                   gggCCGAGCGCCTGTGTGTCAGCCATTCGGCCCTGAGACGGAAGCTGGCTCACCTGAGTGAAGCCTTCGCTGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGTTCGAGCCTTTGCGCCAg                                               <  4:54033/148‑1 (MQ=255)
                                                                                                           cgcCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCga                                       <  4:176446/149‑1 (MQ=255)
                                                                                                            tccTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCgaa                                      <  4:365126/148‑1 (MQ=255)
                                                                                                              cTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGc                                                                                   >  6:328636/1‑102 (MQ=255)
                                                                                                              cTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGc                                                                                   <  5:328636/102‑1 (MQ=255)
                                                                                                                gggTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCg                                        <  3:270680/143‑1 (MQ=255)
                                                                                                                gggTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGACGAAGCCCTGGTCGAGCGTTTGAGCCAGGTGCTCgaatatc                                  >  4:270680/1‑145 (MQ=255)
                                                                                                                 ggTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGc                                                                           >  7:204321/1‑107 (MQ=255)
                                                                                                                 ggTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGc                                                                           <  8:204321/107‑1 (MQ=255)
                                                                                                                 ggTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCAGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAAc                                  <  6:364004/148‑1 (MQ=255)
                                                                                                                     tCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCAc                             >  2:378849/1‑149 (MQ=255)
                                                                                                                        gCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAg                                                                                                       <  6:7363/72‑1 (MQ=255)
                                                                                                                        gCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAg                                                                                                       >  5:7363/1‑72 (MQ=255)
                                                                                                                        gCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGt                                                               >  6:91887/1‑112 (MQ=255)
                                                                                                                        gCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGt                                                               <  5:91887/112‑1 (MQ=255)
                                                                                                                                 aggtGAGCCAGCTGCCGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTg                                                     <  5:284920/110‑1 (MQ=255)
                                                                                                                                 aggtGAGACAGCTGCCGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTg                                                     >  6:284920/4‑113 (MQ=255)
                                                                                                                                    tGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGACCTGCACGCGCTGCGCCAGccc              >  4:217548/1‑149 (MQ=255)
                                                                                                                                        ccGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCg             <  7:194525/146‑1 (MQ=255)
                                                                                                                                        ccGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCg             >  8:194525/1‑146 (MQ=255)
                                                                                                                                           gCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCGCCATCt       >  3:36633/1‑149 (MQ=255)
                                                                                                                                                tGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCGCCATCTTCGAc  >  6:220553/1‑149 (MQ=255)
                                                                                                                                                tGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCGCCATCTTCGAc  >  2:250218/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCGCCATCTTCGAC  >  NC_002947/2896665‑2896957

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