Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,613,529 G→C 8.0% E160D (GAG→GAC)  PP_4933 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,613,5290GC 8.0% 228.8 / 7.4 73E160D (GAG→GAC) PP_4933hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/32);  new base C (2/4);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GCCGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCTTGTTC  >  NC_002947/5613386‑5613676
                                                                                                                                                |                                                                                                                                                    
gCCGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCAGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCg                                                                                                                                                  <  1:54972/149‑1 (MQ=255)
    tACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGcc                                                                                                                                             <  6:26928/149‑1 (MQ=255)
      cAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCACCGGACTCCGGAGccag                                                                                                                                           <  5:255938/149‑3 (MQ=255)
      cAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCACCGGACTCCGGAGcca                                                                                                                                            <  4:282630/148‑2 (MQ=255)
       aGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCa                                                                                                                                         >  1:94057/1‑149 (MQ=255)
         gAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCACCGGACTCCGGAGccag                                                                                                                                           <  4:123768/146‑3 (MQ=255)
         gAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCACCGGACTCCGGAGccag                                                                                                                                           >  3:123768/1‑144 (MQ=255)
           aGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCATTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGcc                                                                                                                                     <  6:419561/149‑1 (MQ=255)
           aGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGC‑GGCACCGGACTCCGGAGCCAGCACgaa                                                                                                                                     <  2:294184/149‑3 (MQ=255)
             tGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCg                                                                                                                                   <  1:143312/149‑1 (MQ=255)
             tGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCg                                                                                                                                   >  6:6204/1‑149 (MQ=255)
               gCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGt                                                                                                                                 >  5:365106/1‑149 (MQ=255)
                     gagTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTg                                                                                                                           <  4:61056/149‑1 (MQ=255)
                        tGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCgg                                                                                                                        >  2:360904/1‑149 (MQ=255)
                              ctgGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACt                                                                                                                  <  6:140114/149‑1 (MQ=255)
                              ctcgcacccccGCAAGGTTGGCTGGTCGCTCCCCTCCACGACACCTCTCTGCCCATCCCCTCCCCCGCCGCCCTGCACTACCGC‑GACACACCCCAGATCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCaaa                                                                                                                    <  4:304947/140‑1 (MQ=255)
                               tgGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGTAAACt                                                                                                                  >  3:297247/1‑148 (MQ=255)
                                 gAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGAt                                                                                                               >  5:238423/1‑149 (MQ=255)
                                 gAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGAt                                                                                                               <  1:360904/149‑1 (MQ=255)
                                  aaaCCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATg                                                                                                              <  7:239573/149‑1 (MQ=255)
                                     ccccGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTa                                                                                                           >  6:135457/1‑149 (MQ=255)
                                      cccGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTAc                                                                                                          <  6:365106/149‑1 (MQ=255)
                                        cGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTg                                                                                                        <  5:6204/149‑1 (MQ=255)
                                                     gggcgggcGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACggg                                                                                           >  6:264068/1‑149 (MQ=255)
                                                       gcgggcgCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGcc                                                                                         <  6:238423/149‑1 (MQ=255)
                                                        cgggcgCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCt                                                                                        <  1:43133/149‑1 (MQ=255)
                                                             gcTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTAtgtg                                                                                   <  5:135457/149‑1 (MQ=255)
                                                               ttCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGAACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCa                                                                                 >  1:281790/1‑149 (MQ=255)
                                                                  aaCGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCCCTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCt                                                                              >  6:21607/1‑149 (MQ=255)
                                                                   aCGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTg                                                                             >  2:37170/1‑149 (MQ=255)
                                                                   aCGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGCCCTATGTGCACCTg                                                                             >  8:115106/1‑149 (MQ=255)
                                                                    cGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTAt                                                                                      <  1:271494/139‑1 (MQ=255)
                                                                    cGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTAt                                                                                      >  2:271494/1‑139 (MQ=255)
                                                                          cTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCg                                                                                                                         >  7:151394/1‑98 (MQ=255)
                                                                          cTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCg                                                                                                                         <  8:151394/98‑1 (MQ=255)
                                                                             tcTGCGCATCCCCTACGCCGACGCCCTGCACTACCGA‑GAAAGACCCCACTTCGTGATTGCTCAGGAGTCCGGTGCC‑GCACTCCCGGTGCATGGCGGCAACCTGATGCTACTGGCGCTCCCACGGGCCTATGTGCACCTGCCCACCGCtt                                                                   <  1:151154/148‑1 (MQ=255)
                                                                             gCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCGCACGGGCCTATGTGCACCTGCCCACCGCtt                                                                   <  1:37170/149‑1 (MQ=255)
                                                                                  ccATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCg                                                                 <  1:20278/146‑1 (MQ=255)
                                                                                  ccATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCg                                                                 >  2:20278/1‑146 (MQ=255)
                                                                                  ccATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTa                                                              <  7:318909/149‑1 (MQ=255)
                                                                                    aTCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTAtgt                                                                                    <  6:296052/125‑1 (MQ=255)
                                                                                    aTCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTAtgt                                                                                    >  5:296052/1‑125 (MQ=255)
                                                                                       cccTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAAATGATCATACTGAAGCGCCCACGGGCAAATGTGCACCTGCCCACCCCTTCGCAACAgat                                                         >  8:176918/1‑147 (MQ=255)
                                                                                        ccTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCt                                                         >  7:169688/1‑148 (MQ=255)
                                                                                           aCGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCa                                                                                               >  4:359489/1‑107 (MQ=255)
                                                                                           aCGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCa                                                                                               <  3:359489/107‑1 (MQ=255)
                                                                                           aCGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATc                                                     >  2:212947/1‑149 (MQ=255)
                                                                                             gccgccgccCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGACAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTa                                                   >  5:372898/1‑149 (MQ=255)
                                                                                               cgccgccCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATg                                                 <  2:281790/149‑1 (MQ=255)
                                                                                                     ccTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGc                                                                                                                       <  4:276899/73‑1 (MQ=255)
                                                                                                     ccTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGc                                                                                                                       >  3:276899/1‑73 (MQ=255)
                                                                                                     ccTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATgc                                           >  7:190269/1‑149 (MQ=255)
                                                                                                     ccTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATgc                                           >  2:58822/1‑149 (MQ=255)
                                                                                                       tGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATgcgc                                         >  3:64291/1‑149 (MQ=255)
                                                                                                         cACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCt                                       >  5:327048/1‑149 (MQ=255)
                                                                                                           cTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTg                                      >  2:176403/1‑148 (MQ=255)
                                                                                                            tACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGc                                     >  8:112027/1‑148 (MQ=255)
                                                                                                              ccGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCACACGGGCTTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCaa                                   <  5:264068/148‑1 (MQ=255)
                                                                                                               cgc‑gAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGt                                 >  3:188961/1‑149 (MQ=255)
                                                                                                                     aaCACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGcccc                            <  4:64291/149‑1 (MQ=255)
                                                                                                                         ccgcACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGa                        <  7:112027/146‑1 (MQ=255)
                                                                                                                           ccACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCgg                          <  3:55923/145‑1 (MQ=255)
                                                                                                                           ccACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCgg                          >  4:55923/1‑145 (MQ=255)
                                                                                                                            cACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCATCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACg                      >  2:256254/1‑148 (MQ=255)
                                                                                                                                 cGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGTAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCaa                  <  4:297247/147‑1 (MQ=255)
                                                                                                                                   tGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGaaa              <  1:176403/149‑1 (MQ=255)
                                                                                                                                     gcggCACCGGACTCCGGAGCCAGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAAc             >  2:274201/3‑149 (MQ=255)
                                                                                                                                         cTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGACAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGc        <  6:372898/149‑1 (MQ=255)
                                                                                                                                          tCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAACCTCAGCt       >  3:322314/1‑149 (MQ=255)
                                                                                                                                          tCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCt       >  5:78535/1‑149 (MQ=255)
                                                                                                                                            cGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCTTg     >  6:120423/1‑149 (MQ=255)
                                                                                                                                            cGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCTTg     >  1:160280/1‑149 (MQ=255)
                                                                                                                                               aGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCTTGTTc  <  2:94057/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCCGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCTTGTTC  >  NC_002947/5613386‑5613676

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