Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,560,900 T→A 7.6% E231D (GAA→GAT)  PP_4048 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,560,9000TA 7.6% 248.8 / 9.3 81E231D (GAA→GAT) PP_4048hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/35);  new base A (3/3);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

GCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGTCAAAAACAG  >  NC_002947/4560753‑4561049
                                                                                                                                                   |                                                                                                                                                     
gCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTc                                                                                                                                                     <  2:1255098/150‑1 (MQ=255)
    gagGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGc                                                                                                                                                 <  1:1124533/150‑1 (MQ=255)
      ggTCACCCGCCAAGCACTCCTGCCCATGGGCCCTGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGGGGGGAGGGCGCTGGCTTGCTCCAGCACCTTGGCCCGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCaa                                                                                                                                               <  2:643541/150‑1 (MQ=255)
        tCACCCGCCGAGTACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCATGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAgg                                                                                                                                             <  2:1202301/150‑1 (MQ=255)
            ccGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCTTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                         >  1:1356665/1‑150 (MQ=255)
            ccGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                         >  1:106077/1‑150 (MQ=255)
            ccGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                         >  1:1220212/1‑150 (MQ=255)
            ccGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATg                                                                                                                                         >  2:1535688/1‑150 (MQ=255)
             cGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCTTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGc                                                                                                                                        <  2:1356665/150‑1 (MQ=255)
                  aGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTg                                                                                                                                   >  2:395575/1‑150 (MQ=255)
                         aTGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCg                                                                                                                            >  1:1176854/1‑150 (MQ=255)
                          tGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGt                                                                                                                           >  1:171573/1‑150 (MQ=255)
                           ggCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTa                                                                                                                          <  1:1535688/150‑1 (MQ=255)
                           ggCCAGGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTa                                                                                                                          <  2:614338/150‑1 (MQ=255)
                             ccATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCGACCTGCGCGGGCAGGGCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTGGa                                                                                                                        >  2:454742/1‑150 (MQ=255)
                              cATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGaa                                                                                                                       <  2:106077/150‑1 (MQ=255)
                               aTGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTACCGCAAGGCATGCACCTGCACGGCGTAGAAc                                                                                                                      >  2:864626/1‑150 (MQ=255)
                                 gggCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTg                                                                                                                    <  1:671659/150‑1 (MQ=255)
                                     tCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCt                                                                                                                >  1:1648187/1‑150 (MQ=255)
                                     tCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCt                                                                                                                >  1:625674/1‑150 (MQ=255)
                                       cGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCg                                                                                                                            <  2:296710/136‑1 (MQ=255)
                                       cGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCg                                                                                                                            >  1:296710/1‑136 (MQ=255)
                                               gcgcATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTTGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTACCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCt                                                                                                      <  1:834835/150‑1 (MQ=255)
                                                 gcATGCACCCCGGCCTGCGCGACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGc                                                                                                    >  2:1019312/1‑150 (MQ=255)
                                                 gcATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACt                                                                                                                     >  2:891954/1‑133 (MQ=255)
                                                 gcATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACt                                                                                                                     <  1:891954/133‑1 (MQ=255)
                                                  cATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCa                                                                                                   >  1:870293/1‑150 (MQ=255)
                                                   aTGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCaa                                                                                                  >  2:562820/1‑150 (MQ=255)
                                                       aCCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCATGCTCCAGCACCTTGTCCAGCCTGGCATCTAGCTGCATTATGGTGATCTTTCTCGCCTCCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGc                                                                                              >  1:779447/1‑150 (MQ=255)
                                                        ccccGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGAAATCCCGCatgcg                                                                                                                                            >  1:1576580/1‑99 (MQ=255)
                                                        ccccGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGAAATCCCGCatgcg                                                                                                                                            <  2:1576580/103‑5 (MQ=255)
                                                         cccGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGt                                                                                                                           <  2:168394/119‑1 (MQ=255)
                                                         cccGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGt                                                                                                                           >  1:168394/1‑119 (MQ=255)
                                                             gCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTccc                                                                                                                                                   <  1:1407499/91‑1 (MQ=255)
                                                             gCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTccc                                                                                                                                                   >  2:1407499/1‑91 (MQ=255)
                                                                 gcgcAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGt                                                                                    >  1:140205/1‑150 (MQ=255)
                                                                 cagcgaCCTGCCCGCGCAGGTTGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGt                                                                                    >  1:213379/6‑150 (MQ=255)
                                                                                gCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGt                                                                     <  2:171573/150‑1 (MQ=255)
                                                                                gCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCAGCATGGAATGCACCTGCACAGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGt                                                                     <  1:1019312/150‑1 (MQ=255)
                                                                                  aGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTg                                                                   <  2:1220212/150‑1 (MQ=255)
                                                                                  aGGTCGCTAGCCAACTTTATCATCTTTGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTg                                                                   >  2:1636593/1‑150 (MQ=255)
                                                                                    gTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGaa                                                                 >  1:809515/1‑150 (MQ=255)
                                                                                     tCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAAc                                                                <  2:1648187/150‑1 (MQ=255)
                                                                                     tCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAAc                                                                <  2:737611/150‑1 (MQ=255)
                                                                                         tGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCAc                                                                                            >  2:37035/1‑118 (MQ=255)
                                                                                         tGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCAc                                                                                            <  1:37035/118‑1 (MQ=255)
                                                                                         tGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGcc                                                            >  1:802429/1‑150 (MQ=255)
                                                                                           gCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAg                                                          <  2:140205/150‑1 (MQ=255)
                                                                                                 tCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGACCTGCAACCGCACCAACAGGTGCGAGGCCGCCTg                                                                       >  2:501587/1‑131 (MQ=255)
                                                                                                 tCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGACCTGCAACCGCACCAACAGGTGCGAGGCCGCCTg                                                                       <  1:501587/131‑1 (MQ=255)
                                                                                                   cAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGc                                                                                                                                        >  2:768897/1‑64 (MQ=255)
                                                                                                   cAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGc                                                                                                                                        <  1:768897/64‑1 (MQ=255)
                                                                                                   cAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCAGGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAg                                                  >  2:1177125/1‑150 (MQ=255)
                                                                                                    aGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCaag                                                 <  2:802429/150‑3 (MQ=255)
                                                                                                     gCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGt                                                >  1:368468/1‑150 (MQ=255)
                                                                                                       aCCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAg                                              <  2:213379/150‑1 (MQ=255)
                                                                                                        ccTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAg                                                  <  2:1084392/145‑1 (MQ=255)
                                                                                                        ccTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAg                                                  >  1:1084392/1‑145 (MQ=255)
                                                                                                          ttGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCAc                                                                                            >  2:1093924/1‑101 (MQ=255)
                                                                                                          ttGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCAc                                                                                            <  1:1093924/101‑1 (MQ=255)
                                                                                                              ccAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAg                                                          <  2:1444100/131‑1 (MQ=255)
                                                                                                              ccAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAg                                                          >  1:1444100/1‑131 (MQ=255)
                                                                                                              ccAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGag                                       <  2:625674/150‑2 (MQ=255)
                                                                                                                  ccTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCtt                                   <  2:809515/150‑1 (MQ=255)
                                                                                                                        aTCGAGCTGCATCAGCGCATGCGGGATTTCACACAAAGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGACTCCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgccgt                             >  1:1207115/1‑150 (MQ=255)
                                                                                                                         tCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgccgtt                            >  1:1388479/1‑150 (MQ=255)
                                                                                                                         tCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGGGTAATGCCGCCGGGTTTAACGGCCATTGGCCCATGTATCTGCGAAGCTttgccgtt                            >  2:904078/1‑150 (MQ=255)
                                                                                                                              cTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGaa                                       >  2:1288885/1‑134 (MQ=255)
                                                                                                                              cTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGaa                                       <  1:1288885/134‑1 (MQ=255)
                                                                                                                                       cgcATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCtt              >  2:233584/1‑150 (MQ=255)
                                                                                                                                        gcATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCCTGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGCACCTTg             <  2:1388479/150‑1 (MQ=255)
                                                                                                                                        gcATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTg             >  2:1129058/1‑150 (MQ=255)
                                                                                                                                         cttGCGGGAAATCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTg             <  2:34734/147‑1 (MQ=255)
                                                                                                                                         cttGCGGGAAATCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTg             >  1:34734/3‑149 (MQ=255)
                                                                                                                                         cttGCGGGAAATCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTACATGGCCTGCAACCGCACCAACAGGTGCGATGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgcc                               <  1:1607262/129‑1 (MQ=255)
                                                                                                                                         cttGCGGGAAATCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTACATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTttgcc                               >  2:1607262/3‑131 (MQ=255)
                                                                                                                                         cATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTgg            <  2:1473533/150‑1 (MQ=255)
                                                                                                                                          aTGCGGGATTTACCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGt           <  1:864626/150‑1 (MQ=255)
                                                                                                                                            gCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGTCa         >  2:1196409/1‑150 (MQ=255)
                                                                                                                                               ggATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTg                                           <  2:494212/113‑1 (MQ=255)
                                                                                                                                               ggATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTg                                           >  1:494212/1‑113 (MQ=255)
                                                                                                                                                   ttCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGTCAAAAACAg  <  2:368468/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCCAGAGGTCACCCGCCCAGCACTCATGGCCATGGGCTCCGGCTTTTGCGCATGCACCCCGGCCTGCGCAACCTGCGCGGGCAGGTCGCTGGCCTGCTCCAGCACCTTGGCCAGCCTGGCATCGAGCTGCATCAGCGCATGCGGGATTTCCCGCAAGGCATGCACCTGCACGGCGTAGAACTGCGCTTCCATGGCCTGCAACCGCACCAACAGGTGCGAGGCCGCCTGGTTGAACGGCCAGTTGCCCAGGTAGCTGCGAAGCTTTGCCGTTGCCGCGAGAACCTTGGTCAAAAACAG  >  NC_002947/4560753‑4561049

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