Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,160,790 C→T 6.2% E105K (GAG→AAG)  PP_4543 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,160,7900CT 6.2% 234.9 / 3.9 80E105K (GAG→AAG) PP_4543hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (41/34);  new base T (2/3);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

CTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCACTGAGAAAAGCATGGCGCAGGAAG  >  NC_002947/5160641‑5160938
                                                                                                                                                     |                                                                                                                                                    
cTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTc                                                                                                                                                      <  1:178742/150‑1 (MQ=255)
         caacaaACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGa                                                                                                                                             <  2:1481785/150‑1 (MQ=255)
           acaaACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGt                                                                                                                                           <  1:1460953/150‑1 (MQ=255)
             aaaCCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGCACAACATAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGc                                                                                                                                         >  2:1141560/1‑150 (MQ=255)
                  gTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACtt                                                                                                                                    >  1:880330/1‑150 (MQ=255)
                          gATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCAGAGCGCCAAGAGGTTTGTGGTAGTAACATAAACCCTCCCTCGCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTAcg                                                                                                                            >  1:1691931/1‑150 (MQ=255)
                            tCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTAcgcg                                                                                                                          >  2:1789649/1‑150 (MQ=255)
                               gCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCgg                                                                                                                         >  1:1614072/1‑148 (MQ=255)
                               gCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCgg                                                                                                                         <  2:1614072/148‑1 (MQ=255)
                                 cGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCaa                                                                                                                     <  2:1029319/150‑1 (MQ=255)
                                   tCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACa                                                                                                                   <  2:752/150‑1 (MQ=255)
                                     agcagcCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAgc                                                                                                                 >  1:1478277/1‑150 (MQ=255)
                                      gcagcCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCgg                                                                                                                         >  2:1757010/1‑141 (MQ=255)
                                      gcagcCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCgg                                                                                                                         <  1:1757010/141‑1 (MQ=255)
                                        agcCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGTGACGCTGTCAACAGCTcc                                                                                                              <  2:119785/150‑1 (MQ=255)
                                             tGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTAc                                                                                                         <  1:993657/150‑1 (MQ=255)
                                                    cTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTTGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCg                                                                                                  >  2:844889/1‑150 (MQ=255)
                                                    cTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTTGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCg                                                                                                  >  2:844806/1‑150 (MQ=255)
                                                    cTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCg                                                                                                  <  1:1932652/150‑1 (MQ=255)
                                                      gACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTc                                                                                                <  2:1896042/150‑1 (MQ=255)
                                                       aCCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTcc                                                                                                   >  2:863510/1‑146 (MQ=255)
                                                       aCCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTcc                                                                                                   <  1:863510/146‑1 (MQ=255)
                                                              tGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGc                                                                                                                                                    >  2:766210/1‑90 (MQ=255)
                                                              tGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGc                                                                                                                                                    <  1:766210/90‑1 (MQ=255)
                                                              tGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCa                                                                                        >  2:1371969/1‑150 (MQ=255)
                                                                  cGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCAc                                                                                    >  1:969820/1‑150 (MQ=255)
                                                                     cacaGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAg                                                                                 <  1:231042/150‑1 (MQ=255)
                                                                     cacaGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAg                                                                                 <  2:814789/150‑1 (MQ=255)
                                                                     cacaGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAg                                                                                 >  2:233432/1‑150 (MQ=255)
                                                                      acaGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAgg                                                                                >  1:379616/1‑150 (MQ=255)
                                                                      acaGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCCCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAgg                                                                                >  2:553922/1‑150 (MQ=255)
                                                                        aGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCAAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTTCGGTACGGGCGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCGGGTg                                                                              >  2:377399/1‑150 (MQ=255)
                                                                         ccgcCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGc                                                                             <  1:1741005/149‑1 (MQ=255)
                                                                            ccAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTcg                                                                          >  1:1470769/1‑150 (MQ=255)
                                                                            ccAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTcg                                                                          >  1:8857/1‑150 (MQ=255)
                                                                                 cAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTTGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCggcggc                                                                     <  1:1789649/150‑1 (MQ=255)
                                                                                  aGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGcc                                                                    >  1:625705/1‑150 (MQ=255)
                                                                                    ttttGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCa                                                                  <  1:1909945/150‑1 (MQ=255)
                                                                                    ttttGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCa                                                                  >  1:621276/1‑150 (MQ=255)
                                                                                      ttGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGc                                                                <  2:871698/150‑1 (MQ=255)
                                                                                          gtagtaGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCgt                                                            >  1:1026094/1‑150 (MQ=255)
                                                                                           tagtagGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCgtg                                                           <  2:1078839/150‑1 (MQ=255)
                                                                                                ggCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTgcg                                                      <  1:517506/150‑1 (MQ=255)
                                                                                                 gCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTgcgc                                                     >  1:305719/1‑150 (MQ=255)
                                                                                                 gCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTgcgc                                                     <  2:621276/150‑1 (MQ=255)
                                                                                                    aTACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGc                                                                >  1:932204/1‑136 (MQ=255)
                                                                                                    aTACCCGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGc                                                                <  2:932204/136‑1 (MQ=255)
                                                                                                        cGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACGCATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGc                                              >  2:1119256/1‑150 (MQ=255)
                                                                                                         ggCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCt                                             >  1:1255159/1‑150 (MQ=255)
                                                                                                         ggCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCt                                             >  1:1890269/1‑150 (MQ=255)
                                                                                                          gCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTg                                            >  2:1865882/1‑150 (MQ=255)
                                                                                                          gCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTg                                            >  2:1016928/1‑150 (MQ=255)
                                                                                                              gCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGc                                        >  2:1291395/1‑150 (MQ=255)
                                                                                                              gCTCACAAGCATCGGCCCAACCATGACCGAAAACATAATCGCCCCGCGAGTGCGAATCGCGGTACTCGTGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCCGGTGCTCGGCGGCCCAGCCCCCGGTGCGCGCCACGCTGCCGc                                        >  1:242686/1‑150 (MQ=255)
                                                                                                                   caAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGtcc                                   <  2:8857/150‑1 (MQ=255)
                                                                                                                    aaGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGc                                                                             >  2:804726/1‑107 (MQ=255)
                                                                                                                    aaGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGc                                                                             <  1:804726/107‑1 (MQ=255)
                                                                                                                     aGCATCGGCCCAACCATGATCGAACACATACTTGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGtcctc                                 <  1:844889/150‑1 (MQ=255)
                                                                                                                     aGCATCGGCCCAACCATGATCGAACACATACTTGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGtcctc                                 <  1:844806/150‑1 (MQ=255)
                                                                                                                       cATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCa                               <  2:1691931/150‑1 (MQ=255)
                                                                                                                        aTCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGc                                                                                                                                         >  2:112142/1‑43 (MQ=255)
                                                                                                                        aTCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGc                                                                                                                                         <  1:112142/43‑1 (MQ=255)
                                                                                                                         tCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCAc                                                <  1:756698/131‑1 (MQ=255)
                                                                                                                         tCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCAc                                                >  2:756698/1‑131 (MQ=255)
                                                                                                                            gCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGc                                           <  2:934406/133‑1 (MQ=255)
                                                                                                                            gCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGc                                           >  1:934406/1‑133 (MQ=255)
                                                                                                                            gCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCa                          >  2:881707/1‑150 (MQ=255)
                                                                                                                              ccAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCACt                        <  1:1865882/150‑1 (MQ=255)
                                                                                                                                   cATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCACTGAGaa                   >  1:473874/1‑150 (MQ=255)
                                                                                                                                   cATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCACTGAGaa                   >  1:1688032/1‑150 (MQ=255)
                                                                                                                                   cATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGAAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTGGGCGGCCCAGCCCGTGTTGTGCGCCACGCTGCTGCTGTCCTCCATTGCACTGAGaa                   >  2:1741750/1‑150 (MQ=255)
                                                                                                                                    aTGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCACTGAGaaa                  <  2:379616/150‑1 (MQ=255)
                                                                                                                                     tGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCAATGAGaaaa                 <  1:1141560/150‑1 (MQ=255)
                                                                                                                                       aTCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCt                                             >  2:1376375/1‑120 (MQ=255)
                                                                                                                                       aTCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCt                                             <  1:1376375/120‑1 (MQ=255)
                                                                                                                                         cGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCACTGAGAAAAGCAt             <  2:1688032/150‑1 (MQ=255)
                                                                                                                                             cacaTACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCACTGAGAAAAGCATGgcg         <  1:684529/150‑1 (MQ=255)
                                                                                                                                               caTACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTAc                                                                                                         >  2:1063939/1‑52 (MQ=255)
                                                                                                                                               caTACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTAc                                                                                                         <  1:1063939/52‑1 (MQ=255)
                                                                                                                                                   cTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCACTGAGAAAAGCATGGCGCAGGaa   >  1:1023855/1‑150 (MQ=255)
                                                                                                                                                    tCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCACTGAGAAAAGCATGGCGCAGGAAg  >  1:1172866/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCGCCCGTACCTGTCCGTCACGCTCCAGCACCAGGTGCTCGGCGGCCCAGCCCGTGTTGCGCGCCACGCTGCCGCTGTCCTCCATTGCACTGAGAAAAGCATGGCGCAGGAAG  >  NC_002947/5160641‑5160938

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