Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,246,232 C→G 5.8% S29S (TCG→TCC alr ← alanine racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,246,2320CG 5.8% 324.3 / 9.6 104S29S (TCG→TCCalralanine racemase
Reads supporting (aligned to +/- strand):  ref base C (47/50);  new base G (4/2);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

Only 100 of 109 total aligned reads displayed.

AGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGTTCGA  >  NC_002947/4246084‑4246381
                                                                                                                                                     |                                                                                                                                                     
aGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCg                                                                                                                                                      <  1:57489/150‑1 (MQ=255)
 ggCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGa                                                                                                                                                     <  1:738256/150‑1 (MQ=255)
    aTCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATc                                                                                                                                                       >  1:1569395/1‑145 (MQ=255)
    aTCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACaa                                                                                                                                                  >  1:907097/1‑150 (MQ=255)
     tCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAAc                                                                                                                                                 >  2:254576/1‑150 (MQ=255)
         ccTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCTTCCAccatttgg                                                                                                                                             <  2:221720/150‑8 (MQ=255)
          cTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACACAGGACTCGGTGTTTTGGATGGAAAGCTTCATGTTGCCTTGGTCCATCGATAACgctgg                                                                                                                                            >  2:1205172/1‑150 (MQ=255)
               gCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCg                                                                                                                                       <  1:872491/150‑1 (MQ=255)
                   ggCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTTACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATa                                                                                                                                   <  1:1855647/150‑1 (MQ=255)
                       cacaGCTTGGACTTGCCGGCCAGCTCGGCCCGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAgg                                                                                                                               <  1:805013/150‑1 (MQ=255)
                        acaGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggg                                                                                                                              >  1:1071174/1‑150 (MQ=255)
                          aGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggggg                                                                                                                            <  1:1571113/150‑1 (MQ=255)
                          aGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCAACAACGGTGGTGCCGCATACAggggg                                                                                                                            <  1:1639297/150‑1 (MQ=255)
                            cTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGc                                                                                                                          <  1:1101799/150‑1 (MQ=255)
                            cTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGc                                                                                                                          >  1:1545681/1‑150 (MQ=255)
                            cTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGc                                                                                                                          <  1:1676352/150‑1 (MQ=255)
                            cTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGc                                                                                                                          >  1:457142/1‑150 (MQ=255)
                            cTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGc                                                                                                                          >  1:451549/1‑150 (MQ=255)
                                 aCTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTc                                                                                                                     <  2:1071174/150‑1 (MQ=255)
                                     gCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGt                                                                                                                 <  2:734633/150‑1 (MQ=255)
                                      ccGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTg                                                                                                                <  2:437110/150‑1 (MQ=255)
                                        ggCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGAt                                                                                                              >  1:1230833/1‑150 (MQ=255)
                                        ggCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGAt                                                                                                              <  2:639613/150‑1 (MQ=255)
                                           cAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAg                                                                                                           >  1:1054525/1‑150 (MQ=255)
                                             gCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGaa                                                                                                         <  1:1384974/150‑1 (MQ=255)
                                             gCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGaa                                                                                                         >  2:2044395/1‑150 (MQ=255)
                                             gCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGaa                                                                                                         <  1:215381/150‑1 (MQ=255)
                                             cgatctccTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACATCGACCCAGGCATTGCTGTTTTGCACGGTCCGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCGTGTCCGGTGATCAGaa                                                                                                         <  2:1225239/144‑1 (MQ=255)
                                                   ccTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCa                                                                                                   <  2:907097/150‑1 (MQ=255)
                                                    cTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGGTGAGTTAACAGATGGCACTGCGGTTTTGCACGTTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGGCCGGTGATCAGAAgagatcg                                                                                                  >  1:1225239/1‑144 (MQ=255)
                                                    cTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAg                                                                                                  <  2:751486/150‑1 (MQ=255)
                                                    cTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAg                                                                                                  <  1:568969/150‑1 (MQ=255)
                                                            tGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGATGCTGTTGTCCATCCCCGACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAg                                                                                          <  1:2124197/150‑1 (MQ=255)
                                                               ggATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCa                                                                                       >  1:432165/1‑150 (MQ=255)
                                                               ggATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCa                                                                                       >  1:113268/1‑150 (MQ=255)
                                                                gATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCgg                                                                                                                  >  1:1985930/1‑122 (MQ=255)
                                                                gATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCgg                                                                                                                  <  2:1985930/122‑1 (MQ=255)
                                                                   gttgtGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAg                                                                                   >  2:2380600/1‑150 (MQ=255)
                                                                      gtgCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGt                                                                                >  1:441562/1‑150 (MQ=255)
                                                                       tgctgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCa                                                                                                   >  2:1248734/1‑130 (MQ=255)
                                                                       tgctgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCa                                                                                                   <  1:1248734/130‑1 (MQ=255)
                                                                         ctgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACg                                                                             <  1:2044395/150‑1 (MQ=255)
                                                                            cAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCGATGGACAACggc                                                                                                                                              >  2:710284/1‑81 (MQ=255)
                                                                            cAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCGATGGACAACggc                                                                                                                                              <  1:710284/82‑2 (MQ=255)
                                                                               ggcgctggcgctGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGaaa                                                                       <  1:983808/150‑1 (MQ=255)
                                                                                    tggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCa                                                                  >  1:621489/1‑150 (MQ=255)
                                                                                        gctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGt                                                              <  2:113268/150‑1 (MQ=255)
                                                                                         ctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggggg                                                                                                                            >  2:1856084/1‑87 (MQ=255)
                                                                                         ctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggggg                                                                                                                            <  1:1856084/87‑1 (MQ=255)
                                                                                         ctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTg                                                             >  1:270264/1‑150 (MQ=255)
                                                                                             tttCGACCCAGGCATTGCTTTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCAACCAGGGGGTCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTgcggg                                                         <  2:2300736/149‑5 (MQ=255)
                                                                                             cTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTa                                                         <  2:270264/150‑1 (MQ=255)
                                                                                              ttCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGa                                                            <  2:1766927/146‑1 (MQ=255)
                                                                                              ttCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGa                                                            <  2:1766725/146‑1 (MQ=255)
                                                                                              ttCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGa                                                            >  1:1766725/1‑146 (MQ=255)
                                                                                              ttCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTAc                                                        >  1:464993/1‑150 (MQ=255)
                                                                                                cGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTc                                                      <  2:1545681/150‑1 (MQ=255)
                                                                                                    ccAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggggg                                                                                                                            >  2:229936/1‑76 (MQ=255)
                                                                                                    ccAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAggggg                                                                                                                            <  1:229936/76‑1 (MQ=255)
                                                                                                     cAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGtt                                                 >  1:1454469/1‑150 (MQ=255)
                                                                                                      aGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCgg                                                                                                                  <  2:1151782/84‑1 (MQ=255)
                                                                                                      aGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCgg                                                                                                                  >  1:1151782/1‑84 (MQ=255)
                                                                                                      aGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAgg                                                                                                 <  2:1662646/101‑1 (MQ=255)
                                                                                                      aGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAgg                                                                                                 >  1:1662646/1‑101 (MQ=255)
                                                                                                         gCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCtt                                              >  1:2125471/1‑150 (MQ=255)
                                                                                                         gCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCtt                                              <  2:2407248/150‑1 (MQ=255)
                                                                                                         gCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCtt                                              <  1:921503/150‑1 (MQ=255)
                                                                                                          cATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTg                                             <  2:432165/150‑1 (MQ=255)
                                                                                                           aTTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGAt                                                                                              <  2:812130/100‑1 (MQ=255)
                                                                                                           aTTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGAt                                                                                              >  1:812130/1‑100 (MQ=255)
                                                                                                               cTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATTTTCAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTa                                        <  1:1205172/150‑1 (MQ=255)
                                                                                                               cTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTa                                        <  1:2172508/150‑1 (MQ=255)
                                                                                                                tGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAg                                       >  1:1473445/1‑150 (MQ=255)
                                                                                                                      gCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCa                                                                                                   <  2:752341/84‑1 (MQ=255)
                                                                                                                      gCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCa                                                                                                   >  1:752341/1‑84 (MQ=255)
                                                                                                                        aCGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGcc                               <  2:2179958/150‑1 (MQ=255)
                                                                                                                              aGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACa                         <  2:1938106/150‑1 (MQ=255)
                                                                                                                                   gtTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGc                    >  1:633223/1‑150 (MQ=255)
                                                                                                                                   gtTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCAAggt                                                                                                <  2:921643/74‑2 (MQ=255)
                                                                                                                                   gtTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCAAggt                                                                                                >  1:921643/1‑73 (MQ=255)
                                                                                                                                     tGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAg                  <  1:1099565/150‑1 (MQ=255)
                                                                                                                                       gTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATTGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgt                >  1:2366322/1‑150 (MQ=255)
                                                                                                                                        caccGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAg                  >  2:1986817/3‑147 (MQ=255)
                                                                                                                                        caccGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAg                  <  1:1986817/145‑1 (MQ=255)
                                                                                                                                        caccGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtg               >  1:903510/3‑150 (MQ=255)
                                                                                                                                         gCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGAGGGCCTGAGCGGTGATCATAAGTGCCATGTATTCAGCCAGAAGGGTACGGCGAACGGGCCTGCTGACTACTAGTGTTCCTGTCCTAGTGCGCGCCCCAACATACGCAGgtgg              <  2:2366322/150‑2 (MQ=255)
                                                                                                                                         accGTTGTCGATGGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAATGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgtgt              >  2:1775956/2‑150 (MQ=255)
                                                                                                                                           cGTTGTCCATCGACAACGGTGGTACCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCGAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgagtgt            >  2:1185497/1‑150 (MQ=255)
                                                                                                                                             ttGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTgc                                    <  2:1488026/124‑1 (MQ=255)
                                                                                                                                             ttGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTgc                                    >  1:1488026/1‑124 (MQ=255)
                                                                                                                                               gTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTGCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGa        >  2:604451/1‑150 (MQ=255)
                                                                                                                                                  cATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGtt     >  2:1740950/1‑150 (MQ=255)
                                                                                                                                                  cATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGtt     >  2:1059181/1‑150 (MQ=255)
                                                                                                                                                  cATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGtt     >  1:1324471/1‑150 (MQ=255)
                                                                                                                                                   aTCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGaa                                                                                    >  2:2027191/1‑70 (MQ=255)
                                                                                                                                                   aTCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGaa                                                                                    <  1:2027191/70‑1 (MQ=255)
                                                                                                                                                   aTCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGTTc    >  2:1324915/1‑150 (MQ=255)
                                                                                                                                                     cGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGTTCGa  <  1:1492046/150‑1 (MQ=255)
                                                                                                                                                     cGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGTTCGa  <  2:457142/150‑1 (MQ=255)
                                                                                                                                                     cGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGTTCGa  <  2:451549/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACGGTGGTGCCGCATACAGGGGGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGTGTTTGAGTTCGA  >  NC_002947/4246084‑4246381

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