Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,896,806 A→T 6.3% Q42L (CAG→CTG)  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,8060AT 6.3% 239.6 / 5.6 79Q42L (CAG→CTG) PP_2551LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (35/39);  new base T (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 = 7.49e-01

CACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCGCCATCT  >  NC_002947/2896665‑2896952
                                                                                                                                             |                                                                                                                                                  
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                                                                                                                                          >  4:82341/1‑142 (MQ=255)
         gAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCAGCTGCCGACTCag                                                                                                                                          <  3:82341/143‑2 (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                                                                                                                                 <  3:87889/107‑1 (MQ=255)
                                                   acttGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGaa                                                                                                                                 >  4:87889/4‑110 (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                                                                          <  3:103597/140‑1 (MQ=255)
                                                                            gCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAAc                                                                          >  4:103597/1‑140 (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                                                                                                                                      <  3:113391/64‑1 (MQ=255)
                                                                                            aGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGc                                                                                                                                      >  4:113391/1‑64 (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                                                                              <  5:328636/102‑1 (MQ=255)
                                                                                                              cTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGc                                                                              >  6:328636/1‑102 (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        >  8:194525/1‑146 (MQ=255)
                                                                                                                                        ccGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCg        <  7:194525/146‑1 (MQ=255)
                                                                                                                                           gCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCGCCATCt  >  3:36633/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CACCGCAGGGAGCCTGCATGAACCCCACCTTCGCCTCCCTCAGCCTCGCCCACCTGCGCACCCTCGACCACCTGTTGCAGCTTAAAAACCTCAGCCACGCGGCCGAGCGCCTGGGTGTCAGCCAATCGGCCCTGAGCCGGCAGCTGGCTCACCTGCGTGAAGCCTTCGATGATCCGTTGCTGGTGCGCCAAGGGCGTGGTTATGTGCTTAGCGAACATGCCGAAGCCCTGGTCGAGCCTTTGCGCCAGGTGCTCGAAGAACTGCACGCGCTGCGCCAGCCCGCCATCT  >  NC_002947/2896665‑2896952

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: