Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,381,315 T→C 6.3% intergenic (+32/+33) bamE → / ← PP_4732 outer membrane protein assembly factor BamE/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,381,3150TC 6.3% 279.2 / 6.5 95intergenic (+32/+33)bamE/PP_4732outer membrane protein assembly factor BamE/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (46/43);  new base C (4/2);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

CGCAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGACAG  >  NC_002947/5381168‑5381464
                                                                                                                                                   |                                                                                                                                                     
cgcAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTg                                                                                                                                                     <  2:1980/150‑1 (MQ=255)
cgcAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTg                                                                                                                                                     <  2:1182170/150‑1 (MQ=255)
cgcAACCGGAAGAAAAACCGGCCAAGCATGGTTAGACAGAAGAGTACATCCAGCGAGAGATCGAAACCACCCAAACTAACCCGGCCCCGACCACGGCACAGCTGGAAACCCGCTATGCGGGAACAAAAAAAAAAAAAAAGCCCGGACTTg                                                                                                                                                     <  2:685294/150‑1 (MQ=255)
 gcAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTgg                                                                                                                                                    >  1:412761/1‑150 (MQ=255)
 gcAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTgg                                                                                                                                                    <  1:1129133/150‑1 (MQ=255)
   aaCCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGTTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTc                                                                                                                                                  <  2:660865/150‑1 (MQ=255)
      cTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACCTCACGTGCAAAACAGCCCGCAGTTGGTCCgg                                                                                                                                               <  2:1256818/150‑1 (MQ=255)
       tGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTggctcttt                                                                                                                                              <  2:1161893/150‑7 (MQ=255)
       tGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCggg                                                                                                                                              <  2:1579094/150‑1 (MQ=255)
       tGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCggg                                                                                                                                              >  2:1908810/1‑150 (MQ=255)
        gaagaaAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGc                                                                                                                                             >  1:1125716/1‑150 (MQ=255)
           gaaAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCttt                                                                                                                                          >  2:833085/1‑150 (MQ=255)
           gaaAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCttt                                                                                                                                          >  1:526899/1‑150 (MQ=255)
             aaaaCCGGCCAAGTCTGGTTTGACCGAAGAGCCCATCCAGCGCGATATCGACACCATCGAGACTACCCCGGTCCCGCCCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCttttt                                                                                                                                        <  1:938250/150‑1 (MQ=255)
                 ccGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTtttgtt                                                                                                                                    <  2:418431/150‑1 (MQ=255)
                  cGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTtttgttt                                                                                                                                   <  2:412761/150‑1 (MQ=255)
                   ggCCCACCCTGGTTCGACCGAAGAGTCAATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATACAAAAAAGCCCGGACCAAGTCCGGGCTTTtttgcatc                                                                                                                                  <  2:1668154/150‑5 (MQ=38)
                                 cGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCa                                                                                                                    <  2:1168372/150‑1 (MQ=255)
                                   aCCGAAGAGTCCATCCAGCGCGATATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAgg                                                                                                                  <  2:1125716/150‑1 (MQ=255)
                                       acgagTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCAtt                                                                                                              >  1:249959/3‑150 (MQ=255)
                                       aaGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCttt                                                                                                                                          <  2:1316864/122‑1 (MQ=255)
                                       aaGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCttt                                                                                                                                          >  1:1316864/1‑122 (MQ=255)
                                       aaGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTgcgc                                                                                                                               <  2:1640349/133‑1 (MQ=255)
                                       aaGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTgcgc                                                                                                                               >  1:1640349/1‑133 (MQ=255)
                                       aaGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTgc                                                                                                                                 <  2:1078186/131‑1 (MQ=255)
                                       aaGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTgc                                                                                                                                 >  1:1078186/1‑131 (MQ=255)
                                                   aGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCTTTtttgcat                                                                                                                                   <  1:2202813/117‑4 (MQ=38)
                                                   aGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACCAAGTCCGGGCTTTtttgcat                                                                                                                                   >  2:2202813/1‑114 (MQ=38)
                                                   aGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACAGCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGc                                                                                                  >  1:1687773/1‑150 (MQ=255)
                                                       cgAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCCGAGCGGGGACATGAGATGTAAAAAAGCCCGGTCTTGGTCCGGGCTTTTTTTTTTTTGCCTTCTACGTTAGGCATTCAATTGCCTGGCTCGc                                                                                              >  2:862508/1‑150 (MQ=255)
                                                        gagaTCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTTCGCCTTATACGTCAGGCATTCAATTTCTTGGCACCCt                                                                                             >  2:235334/1‑150 (MQ=255)
                                                        gagaTCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCt                                                                                             >  1:533261/1‑150 (MQ=255)
                                                            tCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAgg                                                                                                                  >  2:1924017/1‑125 (MQ=255)
                                                            tCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAgg                                                                                                                  <  1:1924017/125‑1 (MQ=255)
                                                              gACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGTGGTTGCTGTATCATTCATTCAATTGCTTGGCACGCTGGCAGa                                                                                       >  2:1355176/1‑150 (MQ=255)
                                                                     tCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATc                                                                                >  2:2230949/1‑150 (MQ=255)
                                                                        agaCTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGAc                                                                             <  2:427261/150‑1 (MQ=255)
                                                                         gaCTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGAcc                                                                            >  1:2107964/1‑150 (MQ=255)
                                                                          aCTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCa                                                                           <  1:1908810/150‑1 (MQ=255)
                                                                            tACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGc                                                                         >  1:1290105/1‑150 (MQ=255)
                                                                            ccccccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGTGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGc                                                                         <  1:903761/148‑1 (MQ=255)
                                                                              ccccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTg                                                                                                                                                     >  2:1836088/1‑72 (MQ=255)
                                                                              ccccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTg                                                                                                                                                     <  1:1836088/72‑1 (MQ=255)
                                                                               cccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTg                                                                                                    <  1:1146959/120‑1 (MQ=255)
                                                                               cccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTg                                                                                                    >  2:1146959/1‑120 (MQ=255)
                                                                               cccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTgg                                                                                           <  2:628716/129‑1 (MQ=255)
                                                                               cccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTgg                                                                                           >  1:628716/1‑129 (MQ=255)
                                                                               cccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTa                                                                      >  1:1417302/1‑150 (MQ=255)
                                                                               cccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCAGGGCTGTTTTGTTTGCGCCTTATACGTCAGGGATTCAATTGCTTGGCACGCTGGCGGAACGCATAGACTAGCgag                                                                      >  2:903761/1‑148 (MQ=255)
                                                                                ccGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTAt                                                                     >  1:1099621/1‑150 (MQ=255)
                                                                                    tCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGc                                                                         <  2:2243641/142‑1 (MQ=255)
                                                                                    tCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGc                                                                         >  1:2243641/1‑142 (MQ=255)
                                                                                      ccGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTcgccg                                                               >  2:628125/1‑150 (MQ=255)
                                                                                        gACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTta                                                                                                                           <  1:1235680/88‑1 (MQ=255)
                                                                                        gACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTta                                                                                                                           >  2:1235680/1‑88 (MQ=255)
                                                                                           cccGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTc                                                                                                                                                  <  1:1867276/62‑1 (MQ=255)
                                                                                           cccGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTc                                                                                                                                                  >  2:1867276/1‑62 (MQ=255)
                                                                                               gCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGa                                                      >  1:741116/1‑150 (MQ=255)
                                                                                                cACCGCTGGAAACCCGCTAATCTGGAACATCAGATGACGAGTAGCCCGGACTTGGTCCGTGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGaa                                                     <  1:1192181/150‑1 (MQ=255)
                                                                                                 aCCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGa                                                                                       >  1:450817/1‑115 (MQ=255)
                                                                                                 aCCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGa                                                                                       <  2:450817/115‑1 (MQ=255)
                                                                                                       ggAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCggg                                              >  2:1648068/1‑150 (MQ=255)
                                                                                                        gAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCg                                                                                                                                                >  2:31727/1‑51 (MQ=255)
                                                                                                        gAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCg                                                                                                                                                <  1:31727/51‑1 (MQ=255)
                                                                                                        gAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTg                                                                                            >  1:1871218/1‑103 (MQ=255)
                                                                                                        gAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTg                                                                                            <  2:1871218/103‑1 (MQ=255)
                                                                                                             ccGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTGAATTGCTTGGCACGCTGGCAGAAc                                                                                     >  1:987175/1‑105 (MQ=255)
                                                                                                             ccGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAAc                                                                                     <  2:987175/105‑1 (MQ=255)
                                                                                                             ccGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAAc                                        >  2:730811/1‑150 (MQ=255)
                                                                                                                cTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGATCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTg                                     <  1:2230949/150‑1 (MQ=255)
                                                                                                                cTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTcgccgc                                                              >  2:444378/1‑125 (MQ=255)
                                                                                                                cTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTcgccgc                                                              <  1:444378/125‑1 (MQ=255)
                                                                                                                    gCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGaa                                                                                      >  2:661089/1‑97 (MQ=255)
                                                                                                                    gCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGaa                                                                                      <  1:661089/97‑1 (MQ=255)
                                                                                                                            aTCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATc                         <  1:833085/150‑1 (MQ=255)
                                                                                                                             tCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCg                        >  1:866133/1‑150 (MQ=255)
                                                                                                                              cAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGACTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCgg                       >  2:459626/1‑150 (MQ=255)
                                                                                                                                 aTGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGccc                    >  1:1679193/1‑150 (MQ=255)
                                                                                                                                   gCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAAt                          >  1:14981/1‑142 (MQ=255)
                                                                                                                                   gCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAAt                          <  2:14981/142‑1 (MQ=255)
                                                                                                                                    cAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGc                 <  2:526899/150‑1 (MQ=255)
                                                                                                                                    cAAAAAAGCCCGGACTTGGTCCGGGATTTTGTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGAACAGCGTATTCGCCGCCTGATTGAATAGCGGGCCAAACGTGGCGCGCACAATCGGCCCGGc                 >  1:654692/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTTAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATATCGGGCCCAACGCCGCGCGCACAACCGGCCCGGCg                >  2:1167025/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGccc                    <  1:1907518/146‑1 (MQ=255)
                                                                                                                                     aaaaaaGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGccc                    >  2:1907518/1‑146 (MQ=255)
                                                                                                                                      aaaaaGCCCGGACCAAGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGt               >  2:416407/1‑150 (MQ=255)
                                                                                                                                        aaaGCCCGGACCAAGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTaa             >  1:1127410/1‑150 (MQ=255)
                                                                                                                                         aaGCCCGGACCAAGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAAt            >  1:160960/1‑150 (MQ=255)
                                                                                                                                             ccGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTTATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGaa        >  2:758015/1‑150 (MQ=255)
                                                                                                                                             ccGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGaa        <  1:1204936/150‑1 (MQ=255)
                                                                                                                                             ccGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGaa        <  2:1290105/150‑1 (MQ=255)
                                                                                                                                              cGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAg       >  1:106340/1‑150 (MQ=255)
                                                                                                                                              cGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAg       <  2:1679193/150‑1 (MQ=255)
                                                                                                                                              cGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAACTGCTTGGCACGCTGGCAGAGCGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAg       <  1:862508/150‑1 (MQ=255)
                                                                                                                                               ggACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAgg      <  1:235334/150‑1 (MQ=255)
                                                                                                                                                   ttGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGACAg  <  2:249959/150‑1 (MQ=255)
                                                                                                                                                   ttGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGACAg  >  1:1016516/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCAACCTGAAGAAAAACCGGCCAAGCCTGGTTCGACCGAAGAGTCCATCCAGCGCGAGATCGACACCATCGAGACTACCCCGGTCCCGACCCCGGCACCGCTGGAAACCCGCTAAGCGGGAACATCAGATGCAAAAAAGCCCGGACTTGGTCCGGGCTTTTTTGTTTGCGCCTTATACGTCAGGCATTCAATTGCTTGGCACGCTGGCAGAACGCATCGACCAGCGTATTCGCCGCCTGATTGAATAGCGGGCCCAACGTGGCGCGCACAATCGGCCCGGCGTAATCGAAGGACAG  >  NC_002947/5381168‑5381464

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