Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,983,083 G→A 7.2% A345A (GCG→GCA)  PP_1775 → metallo‑beta‑lactamase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,983,0830GA 7.2% 206.1 / 3.1 69A345A (GCG→GCA) PP_1775metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (2/3);  total (34/35)
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.53e-01

AGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAGGGGACGCCGGGGGCGGTAATCCAGGCGAGTGAAGGGG  >  NC_002947/1982938‑1983231
                                                                                                                                                 |                                                                                                                                                    
aGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATc                                                                                                                                                   <  8:16623/149‑1 (MQ=255)
aGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATc                                                                                                                                                   <  4:123196/149‑1 (MQ=255)
    gCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCAATcacga                                                                                                                                               >  5:302465/1‑145 (MQ=255)
     cACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGa                                                                                                                                               <  4:264818/148‑1 (MQ=255)
       cGAATATTTGATTGCAGATGCGAGGGCTGGTGTTGCAGCTGGGCGAGATCCACTCGGCTTCAGCCAATTGATCTGTGTGTACACCCATGCCCGCCAACAGCAGGTTGTCAATTACATCAAGAGCACGGGGCGACCGGCGATCGTTATTg                                                                                                                                            <  4:291281/149‑1 (MQ=255)
          atatTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCg                                                                                                                                         <  8:221464/149‑1 (MQ=255)
           tatTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCAATCACGATCGCCgg                                                                                                                                        <  6:302465/149‑1 (MQ=255)
            atTGGAATGCAGAAGCCAGGGCGCGGGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCAATCACGATCGCCggt                                                                                                                                       <  8:327483/149‑2 (MQ=255)
             ttGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCAATCACGATCGCCggt                                                                                                                                       <  8:145015/148‑2 (MQ=255)
               ggAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCaa                                                                                                                                     <  2:321961/148‑1 (MQ=255)
                gAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGa                                                                                                                                               <  4:138368/137‑1 (MQ=255)
                             aGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCggg                                                                                                                      >  1:213409/1‑149 (MQ=255)
                              gggCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGt                                                                                                                     <  5:337643/149‑1 (MQ=255)
                                  gcgTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCgg                                                                                                                       >  1:400884/1‑143 (MQ=255)
                                  gcgTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCgg                                                                                                                       <  2:400884/143‑1 (MQ=255)
                                   cgTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGgcg                                                                                                                >  5:378496/1‑149 (MQ=255)
                                    gtgtTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGgcg                                                                                                                >  3:191325/1‑148 (MQ=255)
                                      gtTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCAt                                                                                                             >  2:369613/1‑149 (MQ=255)
                                          cagcagGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTc                                                                                                         <  2:213409/149‑1 (MQ=255)
                                            gcagGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCaa                                                                                                       >  5:207871/1‑149 (MQ=255)
                                            gcagGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCaa                                                                                                       <  2:153808/149‑1 (MQ=255)
                                              agGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAAt                                                                                                      >  4:289252/1‑148 (MQ=255)
                                                     gATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCAt                                                                                                             <  3:254567/134‑1 (MQ=255)
                                                     gATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCAt                                                                                                             >  4:254567/1‑134 (MQ=255)
                                                       tCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAgg                                                                                            <  6:89674/149‑1 (MQ=255)
                                                         cACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGc                                                                                           >  2:320662/1‑148 (MQ=255)
                                                          aCTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAAt                                                                                                                                    >  6:309169/1‑106 (MQ=255)
                                                          aCTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAAt                                                                                                                                    <  5:309169/106‑1 (MQ=255)
                                                          aCTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCa                                                                                         >  3:262559/1‑149 (MQ=255)
                                                          aCTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCa                                                                                         >  2:197071/1‑149 (MQ=255)
                                                              ggCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCt                                                                                     >  5:195079/1‑149 (MQ=255)
                                                                 ttCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTggg                                                                                  >  8:342123/1‑149 (MQ=255)
                                                                      ccAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATc                                                                             >  3:270050/1‑149 (MQ=255)
                                                                      ccAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATc                                                                             >  5:353109/1‑149 (MQ=255)
                                                                         aTTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCcgc                                                                          <  1:369613/149‑1 (MQ=255)
                                                                               cTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCAt                                                                     >  8:76119/1‑148 (MQ=255)
                                                                               cTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATg                                                                    >  8:169468/1‑149 (MQ=255)
                                                                                  tgtgGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAAtt                                                                                                     >  1:175106/1‑113 (MQ=255)
                                                                                  tgtgGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAAtt                                                                                                     <  2:175106/113‑1 (MQ=255)
                                                                                  tgtgGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATg                                                                    <  1:197071/146‑1 (MQ=255)
                                                                                    tgGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGt                                                                <  7:342123/148‑1 (MQ=255)
                                                                                    tgGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGt                                                                <  1:165654/148‑1 (MQ=255)
                                                                                    tgGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGt                                                                >  2:165654/1‑148 (MQ=255)
                                                                                         acCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTgg                                                                                   >  4:222332/1‑124 (MQ=255)
                                                                                         acCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTgg                                                                                   <  3:222332/124‑1 (MQ=255)
                                                                                           ccATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGTCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCg                                                        >  4:397325/1‑149 (MQ=255)
                                                                                           ccATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCg                                                        >  7:85550/1‑149 (MQ=255)
                                                                                           ccATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCg                                                        >  5:139359/1‑149 (MQ=255)
                                                                                           ccATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCg                                                        >  4:277098/1‑149 (MQ=255)
                                                                                             aTGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTg                                                      <  4:191325/149‑1 (MQ=255)
                                                                                                 ccGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGTTCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGt                                                  <  6:353109/149‑1 (MQ=255)
                                                                                                   gCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTc                                                                                                         >  1:111473/1‑92 (MQ=255)
                                                                                                   gCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTc                                                                                                         <  2:111473/92‑1 (MQ=255)
                                                                                                        cagcagGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGc                                           <  1:320662/149‑1 (MQ=255)
                                                                                                           cagGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTaa                                        <  6:139359/149‑1 (MQ=255)
                                                                                                           cagGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTaa                                        <  6:207871/149‑1 (MQ=255)
                                                                                                             ggtggtCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAgg                                      <  4:270050/149‑1 (MQ=255)
                                                                                                              gtggtCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAggg                                     <  3:277098/149‑1 (MQ=255)
                                                                                                               tggtCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAgggg                                    >  7:274106/1‑149 (MQ=255)
                                                                                                                ggtCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAGGGGa                                   >  6:281184/1‑149 (MQ=255)
                                                                                                                   cAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAGGGGACGc                                >  5:265989/1‑149 (MQ=255)
                                                                                                                    aaTTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCAt                                                                                                                               >  5:286624/1‑53 (MQ=255)
                                                                                                                    aaTTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCAt                                                                                                                               <  6:286624/53‑1 (MQ=255)
                                                                                                                        aCCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAGGGGACGCCggg                            <  6:265989/148‑1 (MQ=255)
                                                                                                                                   aCCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAGGGGACGCCGGGGGCGGTAATc                  >  7:158939/1‑147 (MQ=255)
                                                                                                                                        gCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAGGGGACGCCGGGGGCGGTAATCCAGGCGa           <  6:195079/149‑1 (MQ=255)
                                                                                                                                         ttgccGGCAATCACGATCGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGAGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAGGGGACGCCGGGGGCGGTAATCCAGGCGAg          >  4:103416/4‑149 (MQ=255)
                                                                                                                                          gACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGTCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAGGGGACGCCGGGGGCGGTAATCCAGGCGAGt         <  3:397325/149‑1 (MQ=255)
                                                                                                                                              ggCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAGGGGACGCCGGGGGCGGTAATCCAGGCGAGTGAAg     >  2:67324/1‑149 (MQ=255)
                                                                                                                                                 gATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAGGGGACGCCGGGGGCGGTAATCCAGGCGAGTGAAgggg  >  2:42133/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGCTGCACGAATATTGGAATGCAGAAGCCAGGGCGCGTGTTGCAGCAGGGCGCGATCCACTCGGCTTCAGCCAATTGATCTGTGTGGACACCCATGCCCGCCACCAGCAGGTGGTCAATTACCTCAAGAGCACCGGGCGACCGGCGATCGTGATTGCCGGCAATGGCATGTGTTCGGGTGGGCGCATCGTCAATTACCTCAAGGCCATGCTGGGGGATCCGCGGCATGAGGTGATGTTCGTGGGGTACCAGGCTAAGGGGACGCCGGGGGCGGTAATCCAGGCGAGTGAAGGGG  >  NC_002947/1982938‑1983231

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