Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,559,538 C→T 6.1% W337* (TGG→TAG) ‡ hisZ ← ATP phosphoribosyltransferase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,559,5380CT 6.1% 331.4 / 7.1 111W337* (TGG→TAG) ‡hisZATP phosphoribosyltransferase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base C (52/52);  new base T (3/4);  total (55/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 111 total aligned reads displayed.

CGTCTTGCTGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGG  >  NC_002947/5559392‑5559681
                                                                                                                                                   |                                                                                                                                                
cGTCTTGCTGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGa                                                                                                                                               <  1:1773742/150‑1 (MQ=255)
cGTCTTGCTGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGa                                                                                                                                               >  2:1294973/1‑150 (MQ=255)
 ttCTTGCTGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGAt                                                                                                                                              <  2:2200643/149‑1 (MQ=255)
 gTCTTGCTGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGAt                                                                                                                                              <  2:1640388/150‑1 (MQ=255)
 gTCTTGCTGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGAt                                                                                                                                              <  1:1370919/150‑1 (MQ=255)
  tcgtGCGGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATg                                                                                                                                             >  2:1844682/4‑150 (MQ=255)
       cTGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCTGGATgccgga                                                                                                                                       <  1:1454287/150‑3 (MQ=255)
          aaTCAATTGCCGATCACAATCCGCCTTAAGAGCAGCACTCAACGGCGGGCCAGGCAGAGCCTGAACCACACGCTGGCCCCCGTTGGCCAACTGGCAGCCTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGCTGCCTCCAgc                                                                                                                                    >  2:904269/1‑149 (MQ=255)
          aaTCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGt                                                                                                                                    >  1:917941/1‑150 (MQ=255)
              aaTTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGc                                                                                                                                >  1:648639/1‑150 (MQ=255)
              aaTTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCAACTGAACAGCATCCAGATGCCGCCAGTTGGc                                                                                                                                >  1:32/1‑150 (MQ=255)
              aaTTGCCCATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGc                                                                                                                                >  2:339963/1‑150 (MQ=255)
                 tGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAAt                                                                                                                             <  1:720147/150‑1 (MQ=255)
                     gATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTTCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGa                                                                                                                         >  2:18423/1‑150 (MQ=255)
                      aTCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGAc                                                                                                                        >  2:2006615/1‑150 (MQ=255)
                          caATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGGCGGCATCTGGATgccggacag                                                                                                                                    <  1:1759596/133‑6 (MQ=255)
                          caATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGGCGGCATCTGGATgccggacag                                                                                                                                    >  2:1759596/1‑128 (MQ=255)
                           aaTCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACATGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCg                                                                                                                   >  2:766615/1‑150 (MQ=255)
                            agccGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGc                                                                                                                  <  1:1844682/148‑1 (MQ=255)
                            aTCCGCCTCGAGAGCAGCACTCAACGGATGGTGGTGGTAGCGCTGAGCGAGACGTTGGCCCTCGTTGCGAAAAAGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGc                                                                                                                  <  1:540355/150‑1 (MQ=255)
                            aTCCGCCTCGAGACAAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACGGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGc                                                                                                                  <  2:6707/150‑1 (MQ=255)
                               cGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCCGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCg                                                                                                               >  1:1188368/1‑150 (MQ=255)
                                 ccTAGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCAACGCCATTCTCCGGCATCCACATGCCGCCAGTTGGCAATACGACCttcgttcggc                                                                                                             <  2:945107/150‑3 (MQ=255)
                                      agagCAGCACTCAATGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAgt                                                                                                        <  1:1821900/150‑1 (MQ=255)
                                         gcagcaCTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGtgtg                                                                                                     >  1:710736/1‑150 (MQ=255)
                                          cagcaCTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGGCTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACTGCCTCCGCTCGCCCCAGTGTGa                                                                                                    <  2:648639/150‑1 (MQ=255)
                                               cTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAgg                                                                                               <  2:710736/150‑1 (MQ=255)
                                               cTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAgg                                                                                               >  1:1949630/1‑150 (MQ=255)
                                                tCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGAc                                                                                                                        <  2:1222263/124‑1 (MQ=255)
                                                tCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGAc                                                                                                                        >  1:1222263/1‑124 (MQ=255)
                                                 cAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGt                                                                                             <  1:339963/150‑1 (MQ=255)
                                                     ggCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTc                                                                                         >  2:1064303/1‑150 (MQ=255)
                                                     ggCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTc                                                                                         <  2:2011346/150‑1 (MQ=255)
                                                       cTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGCCCGGCATCCAGATGCCGCCCGTTGGCAATAAGACCTCCGCTCGGACCGGTGTGCCCAGGGTCCTcga                                                                                       <  2:1900825/150‑3 (MQ=255)
                                                              aGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGt                                                                                 <  1:1895035/150‑1 (MQ=255)
                                                                gCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTgg                                                                               >  1:273012/1‑150 (MQ=255)
                                                                gCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTgg                                                                               >  1:778286/1‑150 (MQ=255)
                                                                gCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCGCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTgg                                                                               <  2:1006591/150‑1 (MQ=255)
                                                                  agagCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGgag                                                                             >  2:2085848/1‑150 (MQ=255)
                                                                     gCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGaaa                                                                          <  2:2031541/150‑1 (MQ=255)
                                                                     gCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGaaa                                                                          <  2:2031858/150‑1 (MQ=255)
                                                                      ccTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAAc                                                                         >  1:483021/1‑150 (MQ=255)
                                                                      ccTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAAc                                                                         >  1:2197340/1‑150 (MQ=255)
                                                                        tGACCCGCACTCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACAACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCa                                                                       <  2:917941/150‑1 (MQ=255)
                                                                        tGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGTCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCa                                                                       >  1:673468/1‑150 (MQ=255)
                                                                         gAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATTCGTGGAGAAACCAg                                                                      >  1:942882/1‑150 (MQ=255)
                                                                          aaCCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGt                                                                     <  2:290046/150‑1 (MQ=255)
                                                                              acacGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGcc                                                                 >  1:1462869/1‑150 (MQ=255)
                                                                              acacGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGcc                                                                 >  2:1634846/1‑150 (MQ=255)
                                                                                      gCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGgcgcgc                                                         >  1:1510549/1‑150 (MQ=255)
                                                                                        cgtCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGGCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGccc                                                       >  2:604860/3‑150 (MQ=255)
                                                                                         cTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCg                                                      >  1:109866/1‑150 (MQ=255)
                                                                                           cGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGc                                                    <  2:1949630/150‑1 (MQ=255)
                                                                                             ttGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGt                                                                                 >  2:2121709/1‑119 (MQ=255)
                                                                                             ttGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGt                                                                                 <  1:2121709/119‑1 (MQ=255)
                                                                                             ttGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCa                                                  >  1:1506184/1‑150 (MQ=255)
                                                                                             ttGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGACCCAGTGTGACCAGGGTCTTCAAATGCGTGGAGAAACAAGTGGCGGGGCCCGCCCGGCCa                                                  >  2:229549/1‑150 (MQ=255)
                                                                                               gcgcAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCact                                                <  2:1462869/150‑3 (MQ=255)
                                                                                               gcgcAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCaaa                                                <  1:2409970/150‑1 (MQ=255)
                                                                                                 gcAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGt                                              >  2:386584/1‑150 (MQ=255)
                                                                                                  cAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTgg                                                                                                                                 <  1:1492822/66‑1 (MQ=255)
                                                                                                  cAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTgg                                                                                                                                 >  2:1492822/1‑66 (MQ=255)
                                                                                                   ctcTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATAc                                                                                                                           <  1:2212093/69‑1 (MQ=255)
                                                                                                    aCTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCgg                                           <  1:1228843/150‑1 (MQ=255)
                                                                                                       ggCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGgcgc                                        <  2:483021/150‑1 (MQ=255)
                                                                                                          aGACTTGCTGCCAGAGGGTCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCGGGGCGCGCCCGGCCAAAGTCGGCGCCGa                                     >  2:437625/1‑150 (MQ=255)
                                                                                                          aGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGa                                     <  1:766615/150‑1 (MQ=255)
                                                                                                            acgtgctgtCAGAGTCCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGata                                   <  2:32/147‑1 (MQ=255)
                                                                                                             cTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGa                                     >  1:783958/1‑147 (MQ=255)
                                                                                                             cTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGa                                     <  2:783958/147‑1 (MQ=255)
                                                                                                              ttGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATAtc                                 <  1:1294973/150‑1 (MQ=255)
                                                                                                              ttGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATAtc                                 >  1:1568279/1‑150 (MQ=255)
                                                                                                                 ctgcCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAgg                                                                                               >  2:168754/1‑85 (MQ=255)
                                                                                                                 ctgcCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAgg                                                                                               <  1:168754/85‑1 (MQ=255)
                                                                                                                  tgcCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCa                                                                       <  1:1848387/108‑1 (MQ=255)
                                                                                                                  tgcCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCa                                                                       >  2:1848387/1‑108 (MQ=255)
                                                                                                                  tgcCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATAtcgtcg                             >  2:1781954/1‑150 (MQ=255)
                                                                                                                    ccAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTa                           >  2:1393356/1‑150 (MQ=255)
                                                                                                                    ccAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTa                           <  1:1781954/150‑1 (MQ=255)
                                                                                                                     cAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTTGCAAGACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAg                          >  2:1069592/1‑150 (MQ=255)
                                                                                                                         gggCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGAc                      <  2:1400817/150‑1 (MQ=255)
                                                                                                                            ccGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGc                   <  1:1278953/150‑1 (MQ=255)
                                                                                                                              gCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTc                                                                                         >  2:101978/1‑78 (MQ=255)
                                                                                                                              gCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTc                                                                                         <  1:101978/78‑1 (MQ=255)
                                                                                                                              gCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCgg                                                               >  1:478001/1‑104 (MQ=255)
                                                                                                                              gCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCgg                                                               <  2:478001/104‑1 (MQ=255)
                                                                                                                               cATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGAc                                                                                                   <  1:1543152/67‑1 (MQ=255)
                                                                                                                                 gcGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAAACCGGGGAAAACCCAATGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTgg              <  1:904269/149‑1 (MQ=255)
                                                                                                                                  cGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTggg             <  1:1393356/150‑1 (MQ=255)
                                                                                                                                   gCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGgcgc                                        <  1:955689/122‑1 (MQ=255)
                                                                                                                                   gCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGc            <  1:1064303/150‑1 (MQ=255)
                                                                                                                                     cACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCCCTTCGCTCATTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcga          <  2:1642593/150‑1 (MQ=255)
                                                                                                                                         tggcGGCATCTGGATGCCGGACAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgatc        <  1:601999/146‑1 (MQ=255)
                                                                                                                                         tggcGGCATCTGGATGCCGGCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGgc                                                             >  2:1443169/4‑95 (MQ=255)
                                                                                                                                         tggcGGCATCTGGATGCCGGCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGgc                                                             <  1:1443169/92‑1 (MQ=255)
                                                                                                                                          tCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTg                                                                                >  2:1242913/1‑75 (MQ=255)
                                                                                                                                          tCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTg                                                                                <  1:1242913/75‑1 (MQ=255)
                                                                                                                                           ccGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGGGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGAtt    <  2:2346563/150‑1 (MQ=255)
                                                                                                                                            cGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTg   >  2:1486577/1‑150 (MQ=255)
                                                                                                                                             ggCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCTTGGGCGATCGATTgg  >  1:936646/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGTCTTGCTGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATTGG  >  NC_002947/5559392‑5559681

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: