Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,521,033 T→A 5.3% L83* (TTG→TAG) ‡ murD → UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,521,0330TA 5.3% 231.5 / 2.3 74L83* (TTG→TAG) ‡murDUDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base A (2/2);  total (36/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.86e-01

TGCCTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAAT  >  NC_002947/1520884‑1521180
                                                                                                                                                     |                                                                                                                                                   
tGCCTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCtt                                                                                                                                                     <  2:1071792/150‑1 (MQ=255)
            cGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCGTGGCCAAGGCCAGGCCggg                                                                                                                                         <  1:1363675/150‑2 (MQ=255)
            cGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGc                                                                                                                                         <  1:461942/150‑1 (MQ=255)
             gCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGc                                                                                                                                                  <  1:161291/140‑1 (MQ=255)
             gCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGc                                                                                                                                                  >  2:161291/1‑140 (MQ=255)
              ccGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCt                                                                                                                                       <  1:680627/150‑1 (MQ=255)
                 aTACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgc                                                                                                                                    <  1:683313/150‑1 (MQ=255)
                     ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGCACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCCCGCCGGCCTTgcagcc                                                                                                                                >  2:1449191/1‑149 (MQ=255)
                     ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTTGCCTTGGCCACGCCGGCCTTgcagca                                                                                                                                >  2:1279383/1‑150 (MQ=255)
                                     gAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtggc                                                                                                                <  1:1400941/150‑1 (MQ=255)
                                      aaCTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtggcg                                                                                                               >  2:1093991/1‑150 (MQ=255)
                                         tGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGcaag                                                                                                                       <  2:391797/139‑4 (MQ=255)
                                         tGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGcaag                                                                                                                       >  1:391797/1‑136 (MQ=255)
                                          ggAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGaa                                                                                                           >  1:120559/1‑150 (MQ=255)
                                                  tGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTg                                                                                                   <  2:361907/150‑1 (MQ=255)
                                                  tGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGCCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTg                                                                                                   <  2:1291490/150‑1 (MQ=255)
                                                      ccGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGa                                                                                               <  2:1479120/150‑1 (MQ=255)
                                                       cGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGAc                                                                                              <  1:1123636/150‑1 (MQ=255)
                                                            ttACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                         >  2:1356692/1‑150 (MQ=255)
                                                             tACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATagag                                                                                        >  2:1625150/1‑150 (MQ=255)
                                                               cccGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCt                                                                                      <  1:1073327/150‑1 (MQ=255)
                                                                ccGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTg                                                                                                                                                    <  1:60519/87‑1 (MQ=255)
                                                                ccGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTg                                                                                                                                                    >  2:60519/1‑87 (MQ=255)
                                                                ccGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTg                                                                                     >  1:51376/1‑150 (MQ=255)
                                                                  gCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGtt                                                                                   <  2:120559/150‑1 (MQ=255)
                                                                    aGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTg                                                                                 >  2:472782/1‑150 (MQ=255)
                                                                      gTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCt                                                                               >  1:509606/1‑150 (MQ=255)
                                                                       tGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTg                                                                                                                                     <  2:90041/95‑1 (MQ=255)
                                                                       tGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTg                                                                                                                                     >  1:90041/1‑95 (MQ=255)
                                                                           aGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCa                                                                          <  2:51376/150‑1 (MQ=255)
                                                                                         aaCTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGa                                                            >  1:184348/1‑150 (MQ=255)
                                                                                          aCTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCAt                                                                         <  1:873375/136‑1 (MQ=255)
                                                                                          aCTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCAt                                                                         >  2:873375/1‑136 (MQ=255)
                                                                                              gACGTGGACTTCCTGTGCCGCGCCAAGGAGCCGAACGTGAACCCCGGGACTGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTTCTCGCCATGCCAAGGCGCCGATCGtt                                                       >  1:205155/1‑150 (MQ=255)
                                                                                              gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTTGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGtt                                                       <  1:1279383/150‑1 (MQ=255)
                                                                                              gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGtt                                                       >  2:1536832/1‑150 (MQ=255)
                                                                                                cGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGc                                                     >  2:864085/1‑150 (MQ=255)
                                                                                                    gACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATc                                                 <  1:1477374/150‑1 (MQ=255)
                                                                                                      cTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAg                                               <  2:205155/150‑1 (MQ=255)
                                                                                                         ccTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTg                                                                                     >  1:1259057/1‑109 (MQ=255)
                                                                                                         ccTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTg                                                                                     <  2:1259057/109‑1 (MQ=255)
                                                                                                         ccTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTgg                                            <  1:1093991/150‑1 (MQ=255)
                                                                                                          cTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGt                                           <  2:509606/150‑1 (MQ=255)
                                                                                                           tgtgCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCa                                                <  2:1645650/144‑1 (MQ=255)
                                                                                                           tgtgCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGCTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGtt                                          <  1:1216612/150‑1 (MQ=255)
                                                                                                           tgggCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCa                                                >  1:1645650/4‑144 (MQ=255)
                                                                                                                cgcgCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCATTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAAc                                     <  2:184348/150‑1 (MQ=255)
                                                                                                                 gcgcCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACg                                    >  1:746700/1‑150 (MQ=255)
                                                                                                                     cAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAACTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCaa                                >  1:108120/1‑150 (MQ=255)
                                                                                                                      aaCGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAag                               >  2:630288/1‑150 (MQ=255)
                                                                                                                       aCGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAaga                              <  1:1625150/150‑1 (MQ=255)
                                                                                                                       aCGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAaga                              <  1:1536832/150‑1 (MQ=255)
                                                                                                                                  gTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCacg                   >  2:1159163/1‑150 (MQ=255)
                                                                                                                                   tGAGCCACGGCCTGGCCTTGGACACGCCGGCCGTGCAGCGGCCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGTGCTTTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGTTTCCAACGCCAAGAGCACCGTCacga                  <  2:623840/150‑1 (MQ=255)
                                                                                                                                    gAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGtt                                                       <  2:638770/112‑1 (MQ=255)
                                                                                                                                    gAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGtt                                                       >  1:638770/1‑112 (MQ=255)
                                                                                                                                    gAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCacgac                 >  1:1618763/1‑150 (MQ=255)
                                                                                                                                         ccGGCGTGGCCAAGGCCAGGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGcc                                                                           >  2:627080/1‑87 (MQ=38)
                                                                                                                                         ccGGCGTGGCCAAGGCCAGGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGcc                                                                           <  1:627080/87‑1 (MQ=38)
                                                                                                                                         ccGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGt                                                                                    <  1:142562/78‑1 (MQ=255)
                                                                                                                                         ccGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGt                                                                                    >  2:142562/1‑78 (MQ=255)
                                                                                                                                          cGGCGTGGCCAAGGCCAGGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGAAATAGAGCTGTTTGCTCGCCATGCCAAGCCGCCGATCGTTGCAATCAGTGGTTCCAACGCCAAGAGCACCGTCCCGAGGCTggc           >  2:279658/1‑149 (MQ=255)
                                                                                                                                           ggcctggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGt                                           >  1:1516225/1‑117 (MQ=255)
                                                                                                                                           ggcctggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGt                                           <  2:1516225/117‑1 (MQ=255)
                                                                                                                                              ctggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTg                                             >  1:642565/1‑112 (MQ=255)
                                                                                                                                              ctggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTg                                             <  2:642565/112‑1 (MQ=255)
                                                                                                                                              ctggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGc               >  1:732765/1‑142 (MQ=255)
                                                                                                                                              ctggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGc               <  2:732765/142‑1 (MQ=255)
                                                                                                                                               tggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCaa                                      >  2:27535/1‑118 (MQ=255)
                                                                                                                                               tggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCaa                                      <  1:27535/118‑1 (MQ=255)
                                                                                                                                                  cctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGttt                                                                                  <  1:106856/71‑1 (MQ=255)
                                                                                                                                                  cctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGttt                                                                                  >  2:106856/1‑71 (MQ=255)
                                                                                                                                                   ctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCa                           >  2:674448/1‑125 (MQ=255)
                                                                                                                                                   ctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCa                           <  1:674448/125‑1 (MQ=255)
                                                                                                                                                   ctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAAt  >  2:249286/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGCCTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAAT  >  NC_002947/1520884‑1521180

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