Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,103,973 A→T 5.4% Q269L (CAA→CTA) ‡ PP_5353 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,103,9730AT 5.4% 285.9 / 4.6 92Q269L (CAA→CTA) ‡PP_5353hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (41/46);  new base T (3/2);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

ATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGT  >  NC_002947/6103828‑6104120
                                                                                                                                                   |                                                                                                                                                   
aTCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCaatg                                                                                                                                                  <  2:250075/149‑3 (MQ=255)
aTCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAgg                                                                                                                                                  >  5:249035/1‑149 (MQ=255)
 tCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATAACGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAgg                                                                                                                                                  <  1:298302/149‑1 (MQ=255)
   gATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCga                                                                                                                                               <  4:59204/149‑1 (MQ=255)
   gATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGCTCGCGCCCAAGGCga                                                                                                                                               <  7:4900/149‑1 (MQ=255)
        ccTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAAc                                                                                                                                           >  4:59887/1‑148 (MQ=255)
         cTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAAc                                                                                                                                           >  5:243608/1‑147 (MQ=255)
                tGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCTCGCCTTGGGCGCGAACCatgtgcca                                                                                                                                  <  8:336343/149‑8 (MQ=255)
                tGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTgg                                                                                                                                  <  4:241805/149‑1 (MQ=255)
                 gCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCTCGCCTTGGGCGCGAACCatgtgcca                                                                                                                                  <  4:307558/148‑8 (MQ=255)
                  cAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGc                                                                                                                                 <  4:354260/148‑1 (MQ=255)
                  cAAATGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTTTCGAAGATTTCGCGCTTCGACGAGGGCGAAGTGTTCAGCGAAAGCTGGCACATTGTTCGCGGCCAAGGCGAGAACCAGCCTTGGc                                                                                                                                 <  3:122498/148‑1 (MQ=255)
                   aaCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTGCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTg                                                                                                                               <  3:34012/149‑1 (MQ=255)
                    aCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              <  6:249035/149‑1 (MQ=255)
                    aCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              >  5:247633/1‑149 (MQ=255)
                          cGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCTCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCGCAAGTCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        <  8:149030/149‑1 (MQ=255)
                          cGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        >  3:181409/1‑149 (MQ=255)
                           ggTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCgg                                                                                                                       <  5:251697/149‑1 (MQ=255)
                               tCGATGACCGCGCCGAGCGGACCGACGCTACCCTGACCTTCCGGGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCTAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCTCCCAAGGCGAGAACCAGCGTTGGCTGGTAGCCGGTAt                                                                                                                    <  4:324506/148‑1 (MQ=255)
                                 gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              <  2:301422/136‑1 (MQ=255)
                                 gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              >  1:301422/1‑136 (MQ=255)
                                 gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCg                                                                                                                 >  4:185334/1‑149 (MQ=255)
                                      ccGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCACGGCGAAGTGTTCAGCGAAACCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCAGTCaa                                                                                                             >  6:174235/1‑148 (MQ=255)
                                          gcCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAg                                                                                                        <  3:132764/149‑1 (MQ=255)
                                          gcCCACCGCACCTAACGTACCCTGTCCTTACGTGGGGTTTCGAATAATTGGGGGTTTGCTCATGGGGAATGTCTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACa                                                                                                         <  1:117811/149‑1 (MQ=255)
                                                  cACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGc                                                                                                >  7:127943/1‑149 (MQ=255)
                                                  cACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGc                                                                                                >  3:296732/1‑149 (MQ=255)
                                                     cGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                             >  5:328896/1‑149 (MQ=255)
                                                     cGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                             >  5:133575/1‑149 (MQ=255)
                                                        cGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTg                                                                                          >  5:70200/1‑149 (MQ=255)
                                                             cccTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaa                                                                                      <  4:287399/148‑1 (MQ=255)
                                                             cccTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaa                                                                                      <  2:346478/148‑1 (MQ=255)
                                                              ccTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaaa                                                                                     >  7:345296/1‑148 (MQ=255)
                                                              ccTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaaa                                                                                     >  2:166267/1‑148 (MQ=255)
                                                                tGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCGTTGTCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                             <  3:351565/138‑1 (MQ=255)
                                                                tGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACAAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTCGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCgacaacgg                                                                                     >  4:351565/1‑139 (MQ=255)
                                                                   ccTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAg                                                                                >  5:342965/1‑148 (MQ=255)
                                                                   ccTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAg                                                                                >  7:346079/1‑148 (MQ=255)
                                                                     ttCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAc                                                                             >  3:63542/1‑149 (MQ=255)
                                                                     ttCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAc                                                                             >  4:276051/1‑149 (MQ=255)
                                                                        cGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAcccc                                                                          <  6:247633/149‑1 (MQ=255)
                                                                            gCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAATCAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                       <  6:243608/148‑1 (MQ=255)
                                                                            gCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                       <  2:254878/148‑1 (MQ=255)
                                                                            gCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                       <  8:132991/148‑1 (MQ=255)
                                                                            gCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTCGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCACGGGc                                                                      >  6:264426/1‑149 (MQ=255)
                                                                                 tCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        >  3:95013/1‑94 (MQ=255)
                                                                                 tCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        <  4:95013/94‑1 (MQ=255)
                                                                                 tCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGcc                                                                 >  1:259659/1‑149 (MQ=255)
                                                                                   gaagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCg                                                               <  4:102833/149‑1 (MQ=255)
                                                                                    aagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                              >  6:190046/1‑149 (MQ=255)
                                                                                    aagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                              >  6:219508/1‑149 (MQ=255)
                                                                                    aagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCAGGTATCCGTCAAAACAGCTAATCGCAGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                              >  6:346944/1‑149 (MQ=255)
                                                                                    aagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCACCGAAAGCTGGCACATGGTAAGCGCCCAAGGGGAGAACCAGCCTTGGCTGGTAGCCGCTATCCGTCAAAACAGCTAATCGCTGCCTGCTAAACACAGACCCACGGGCTTGCCCgg                                                              >  4:31648/1‑149 (MQ=255)
                                                                                     agaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                             <  8:346079/149‑1 (MQ=255)
                                                                                     agaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                             >  1:144031/1‑149 (MQ=255)
                                                                                     agaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                             <  6:342965/149‑1 (MQ=255)
                                                                                     agaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCAAGCCCggg                                                             <  4:296732/149‑1 (MQ=255)
                                                                                      gaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgggg                                                            >  2:95508/1‑149 (MQ=255)
                                                                                       aaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGt                                                           <  6:70200/149‑1 (MQ=255)
                                                                                             gcgcTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCGGGGttctgt                                                      <  1:145511/149‑5 (MQ=255)
                                                                                                cTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgctt                                                   >  5:161083/1‑149 (MQ=255)
                                                                                                   cGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgcttt                                                  <  5:190046/147‑1 (MQ=255)
                                                                                                    gACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgcttttg                                                <  4:181409/148‑1 (MQ=255)
                                                                                                            cGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTATTCTTTTGCCAAAGTgg                                       >  7:7162/1‑149 (MQ=255)
                                                                                                              aaGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAac                                                                                    >  5:94509/1‑102 (MQ=255)
                                                                                                              aaGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAac                                                                                    <  6:94509/102‑1 (MQ=255)
                                                                                                                   gtTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                                  >  7:17414/1‑148 (MQ=255)
                                                                                                                   gtTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                                  >  5:213418/1‑148 (MQ=255)
                                                                                                                    ttCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                                  >  6:14685/1‑147 (MQ=255)
                                                                                                                       aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAggg                             <  3:276051/149‑1 (MQ=255)
                                                                                                                       aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  6:180329/148‑1 (MQ=255)
                                                                                                                       aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  8:345296/148‑1 (MQ=255)
                                                                                                                       aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  2:259659/148‑1 (MQ=255)
                                                                                                                       aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  5:219508/148‑1 (MQ=255)
                                                                                                                           aaaGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              >  2:34280/1‑144 (MQ=255)
                                                                                                                           aaaGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGCTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  1:34280/144‑1 (MQ=255)
                                                                                                                             aGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAg                       <  5:310025/149‑1 (MQ=255)
                                                                                                                              gCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAgg                      <  3:237052/149‑1 (MQ=255)
                                                                                                                               cTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgc                     >  8:160733/1‑149 (MQ=255)
                                                                                                                                  aaggctGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgcg                    >  2:372894/6‑147 (MQ=255)
                                                                                                                                   aggctGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACACCGGGCTTGGCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAg                 >  5:242158/5‑149 (MQ=255)
                                                                                                                                    acaTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGc                >  2:127263/1‑149 (MQ=255)
                                                                                                                                    acaTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGc                <  8:17414/149‑1 (MQ=255)
                                                                                                                                       tGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGtt             >  8:99709/1‑149 (MQ=255)
                                                                                                                                         gTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt           <  5:224881/149‑1 (MQ=255)
                                                                                                                                         gTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt           >  6:71724/1‑149 (MQ=255)
                                                                                                                                          ttCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt           >  7:128062/1‑148 (MQ=255)
                                                                                                                                           tCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCTGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTc          <  6:328896/148‑1 (MQ=255)
                                                                                                                                             gcgcCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATc       <  6:213418/149‑1 (MQ=255)
                                                                                                                                                cccAATGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAg    <  3:59887/149‑1 (MQ=255)
                                                                                                                                                 ccAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATACGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCAGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAgg   >  5:298397/1‑149 (MQ=255)
                                                                                                                                                  cAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGt  <  8:128062/149‑1 (MQ=255)
                                                                                                                                                  cAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGt  <  1:166267/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGT  >  NC_002947/6103828‑6104120

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