Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,967,204 T→A 5.6% intergenic (+31/+88) PP_1762 → / ← yciK hypothetical protein/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,967,2040TA 5.6% 226.2 / 4.2 72intergenic (+31/+88)PP_1762/yciKhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base A (2/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

GAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCAT  >  NC_002947/1967057‑1967342
                                                                                                                                                    |                                                                                                                                          
ccaggcctaccccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCtt                                                                                                                                           >  1:1126583‑M1/133‑150 (MQ=255)
       taccccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCt                                                                                                                                    >  2:733665‑M1/126‑150 (MQ=255)
       gcAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTAGCGCGCTGGCCCGCgc                                                                                                                                   >  2:659709‑M1/1‑148 (MQ=255)
        accccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttgaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                   <  1:737201‑M1/26‑1 (MQ=255)
          cccggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTccc                                                                                                                                 <  2:878717‑M1/28‑1 (MQ=255)
            cggcgagaacccgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCcac                                                                                                                               <  2:464774‑M1/30‑1 (MQ=255)
               aTCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCC                                                                                                                                       >  1:330065‑M1/1‑138 (MQ=255)
               aTCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCC                                                                                                                                       <  2:330065‑M1/138‑1 (MQ=255)
                 cTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATc                                                                                                                         <  1:1398271/150‑1 (MQ=255)
                    cAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggc                                                                                                                      >  2:1345635‑M1/1‑144 (MQ=255)
                               ccTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGTTCAGCATGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                           <  1:798743/149‑1 (MQ=38)
                               aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                           <  2:375802/150‑1 (MQ=255)
                               aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                           <  2:425710/150‑1 (MQ=255)
                                  gCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGa                                                                                                        <  1:1067296/150‑1 (MQ=255)
                                    gctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGAt                                                                                                       <  2:1247556‑M1/54‑1 (MQ=255)
                                      gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCgaagag                                                                                                                               <  1:933910‑M1/123‑7 (MQ=255)
                                      gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCgaagag                                                                                                                               >  2:933910‑M1/1‑117 (MQ=255)
                                         cctgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGcc                                                                                                 <  1:956136/150‑1 (MQ=255)
                                           tgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCa                                                                                                                        >  2:521788/1‑125 (MQ=255)
                                           tgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCa                                                                                                                        <  1:521788/125‑1 (MQ=255)
                                            gttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGGTCACCGCAACCCt                                                                                                              >  1:91116/1‑134 (MQ=255)
                                            gttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCt                                                                                                              <  2:91116/134‑1 (MQ=255)
                                            gttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCgaaga                                                                                                                                <  1:1120314‑M1/116‑6 (MQ=255)
                                            gttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCgaaga                                                                                                                                >  2:1120314‑M1/1‑111 (MQ=255)
                                             ttCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGc                                                                                             >  2:1310949‑M2/1‑149 (MQ=255)
                                             ttCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGc                                                                                             >  1:566885‑M2/1‑149 (MQ=255)
                                                  ggAGGGCTGCTTGCGGCTGGAAGGGCGCCAGAAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaacccttgagggttgcggtgatcctgtg                                                                                        <  1:1045926‑M1/150‑37 (MQ=255)
                                                  ggAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtg                                                                                        >  1:1245693‑M2/1‑144 (MQ=255)
                                                   gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtgg                                                                                       <  1:1345635‑M1/150‑38 (MQ=255)
                                                   gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgg                                                                                       >  1:141177‑M2/1‑143 (MQ=255)
                                                   gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGCATCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtgg                                                                                       <  1:1083164‑M1/150‑38 (MQ=255)
                                                      gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                     >  1:464774‑M1/79‑150 (MQ=255)
                                                      cgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCC                                                                                                                                        >  1:30483‑M1/1‑98 (MQ=255)
                                                      cgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCC                                                                                                                                        <  2:30483‑M1/98‑1 (MQ=255)
                                                       ccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAg                                                                                    >  1:977503‑M1/78‑150 (MQ=255)
                                                         tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtgggagcgg                                                                                 >  1:464706‑M1/1‑107 (MQ=255)
                                                         tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtgggagcgg                                                                                 <  2:464706‑M1/150‑44 (MQ=255)
                                                         tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtg                                                                                        >  1:72148‑M1/1‑107 (MQ=255)
                                                         tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtg                                                                                        <  2:72148‑M1/143‑37 (MQ=255)
                                                         tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                 >  1:1337093‑M2/1‑137 (MQ=255)
                                                          gtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                 >  1:1247556‑M1/75‑150 (MQ=255)
                                                          gggtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                 >  2:300621‑M1/75‑150 (MQ=255)
                                                          atctacctgtaccttatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                 >  1:590879‑M1/75‑150 (MQ=255)
                                                               cctgtacgtgatggggccggacagcacaggggtggacgggcaggcgttcagcgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCAGGGCAACCCTCACGGATTGGGCCCTGCCTGTGGgagcgggagag                                                                            <  2:590879‑M1/81‑5 (MQ=255)
                                                                   tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtgg                                                                                       <  2:1241496‑M1/134‑38 (MQ=255)
                                                                   tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                 <  2:1197883‑M2/140‑14 (MQ=255)
                                                                   tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                 >  1:1197883‑M2/1‑127 (MQ=255)
                                                                   tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTAGCTGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtgg                                                                                       >  1:1241496‑M1/1‑97 (MQ=255)
                                                                         ggCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagag                                                                 >  1:659709‑M2/1‑121 (MQ=255)
                                                                           cgcCAGCAAGGCCGTATGCGGCCGGCGCTGTCTGCGGGTTGACCGCTTCCCTGTACTGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCGCAACCATCAAGGATTGGGCCCTGCCTGTGGGAGCGGCCAAGCCCGCGAAGAGGc                                                               <  2:713316/150‑1 (MQ=37)
                                                                              gccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                         <  1:283335‑M1/68‑1 (MQ=255)
                                                                              gccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                         >  2:283335‑M1/55‑122 (MQ=255)
                                                                                gCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                 >  1:586647‑M2/1‑114 (MQ=255)
                                                                                gCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                 <  2:586647‑M2/127‑14 (MQ=255)
                                                                                         caggggtgaacgggcaggcgttcaacgcccaataaacctgcacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                       <  2:1201420‑M1/70‑1 (MQ=255)
                                                                                         caggggtgaacgggcaggcgttcaacgcccaataaacctgcacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                       >  1:1201420‑M1/44‑113 (MQ=255)
                                                                                              gtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                     >  1:816858‑M1/39‑110 (MQ=255)
                                                                                              gtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                     >  1:816926‑M1/39‑110 (MQ=255)
                                                                                              gtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                     <  2:816926‑M1/72‑1 (MQ=255)
                                                                                              gtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                     <  2:816858‑M1/72‑1 (MQ=255)
                                                                                                                   gcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                     >  1:1126485‑M1/18‑89 (MQ=255)
                                                                                                                   gcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                     <  2:1126485‑M1/72‑1 (MQ=255)
                                                                                                                         cTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                 >  2:40373‑M2/1‑74 (MQ=255)
                                                                                                                         cTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                 <  1:40373‑M2/87‑14 (MQ=255)
                                                                                                                          tttCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                                                         <  1:269533‑M2/78‑6 (MQ=255)
                                                                                                                          tttCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                                                         >  2:269533‑M2/1‑73 (MQ=255)
                                                                                                                              ctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCcgtg                                                                                        >  1:201985‑M1/7‑75 (MQ=255)
                                                                                                                              ctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGc                                                                                   >  1:283002‑M1/7‑80 (MQ=255)
                                                                                                                              ctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGc                                                                                   <  2:283002‑M1/74‑1 (MQ=255)
                                                                                                                                        ctctTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcc   >  1:659709‑M1/1‑29 (MQ=255)
                                                                                                                                        CTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCa   <  1:751992‑M1/150‑1 (MQ=255)
                                                                                                                                         TCTTCGCGGGCTTGCCCGCTCCCCCAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCGGGGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCGGGTTTATTGGGCGTTGAAAGCCAGCCCGTTCACCCCTGTGCTGTCCGGCCCCAt  <  2:737201‑M1/150‑1 (MQ=255)
                                                                                                                                         TCTTCGCGGGCTTGCCCGCTCCCACAGGCTCACCACCACCCTCCAGGATTGCCCGCTGCCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCAt  <  1:249233‑M1/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCAT  >  NC_002947/1967057‑1967342

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: