Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,521,227 A→G 5.5% S171P (TCG→CCG)  icd ← NADP(+)‑specific isocitrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,521,2270AG 5.5% 218.2 / 2.0 73S171P (TCG→CCG) icdNADP(+)‑specific isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base G (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

CCACGTATTGCAGGGCCTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGCATACATACAGGTCCAGTTGCTG  >  NC_002947/4521080‑4521374
                                                                                                                                                    |                                                                                                                                                     
ccacGTATTGCAGGGCCTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGAc                                                                                                                                                      >  4:6773/1‑149 (MQ=255)
ccacGTATTGCAGGGCCTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGAc                                                                                                                                                      >  1:46508/1‑149 (MQ=255)
       ttGCAGGGCCTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGcc                                                                                                                                                <  2:181709/148‑1 (MQ=255)
          cAGGGCCTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGCTGACCCCCATTTCCTCCTTGAGAATCTTGATCACCTTGTTCGGCTCGGGCGACCCGGCCTTc                                                                                                                                             <  6:276306/148‑1 (MQ=255)
               ccTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGGCCCGGCCTTCCACTc                                                                                                                                        >  2:242448/1‑148 (MQ=255)
               ccTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGGCCCGGCCTTCCACTc                                                                                                                                        >  4:138593/1‑148 (MQ=255)
                cTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCg                                                                                                                                       >  6:357413/1‑148 (MQ=255)
                cTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCg                                                                                                                                       <  5:357413/148‑1 (MQ=255)
                cTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGa                                                                                                                                      <  3:342002/149‑1 (MQ=255)
                      cACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCgg                                                                                                                                >  5:125994/1‑149 (MQ=255)
                        ccAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCa                                                                                                                              >  1:203801/1‑149 (MQ=255)
                         cAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCa                                                                                                                              >  5:23524/1‑148 (MQ=255)
                         cAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCa                                                                                                                              >  5:175108/1‑148 (MQ=255)
                                   gTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTc                                                                                                                 <  6:23524/149‑1 (MQ=255)
                                     aCCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTGGGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCACGGGCGACCCGCCCTTCCACTCGATGCCGGCATAGATGTCTTCg                                                                                                                >  1:268950/1‑149 (MQ=255)
                                      cccTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGa                                                                                                              <  3:36386/149‑1 (MQ=255)
                                       cctcctTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAg                                                                                                             >  4:272853/1‑149 (MQ=255)
                                        ctcctTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGt                                                                                                            >  1:101970/1‑149 (MQ=255)
                                            ttCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTtctc                                                                                                        >  4:265668/1‑149 (MQ=255)
                                               gACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCGTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTcgcg                                                                                                     >  5:328175/1‑149 (MQ=255)
                                               gACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTcgc                                                                                                      <  2:268950/148‑1 (MQ=255)
                                               gACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTcgc                                                                                                      <  3:6773/148‑1 (MQ=255)
                                               gACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTcgc                                                                                                      <  3:272853/148‑1 (MQ=255)
                                                acacCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTcgcg                                                                                                     >  3:190027/1‑148 (MQ=255)
                                                   ccGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCGCTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAg                                                                                                 <  5:238026/149‑1 (MQ=255)
                                                           aCGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGa                                                                                                   >  7:77445/1‑139 (MQ=255)
                                                           aCGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGa                                                                                                   <  8:77445/139‑1 (MQ=255)
                                                             gCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCTCCATGt                                                                                       >  8:209225/1‑149 (MQ=255)
                                                                  tACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACg                                                                                  <  4:190027/149‑1 (MQ=255)
                                                                                        ggATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGGCCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGTAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCggg                                                             <  3:138593/149‑1 (MQ=255)
                                                                                        ggATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGGCCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCACCCGGCTTTTTCACCggg                                                             <  1:242448/149‑1 (MQ=255)
                                                                                               ggTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCttttt                                                                    >  7:69593/1‑135 (MQ=255)
                                                                                               ggTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCttttt                                                                    <  8:69593/135‑1 (MQ=255)
                                                                                               ggTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTgg                                                        <  2:101970/147‑1 (MQ=255)
                                                                                               ggTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGca                                                      >  1:142659/1‑149 (MQ=255)
                                                                                                gTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGcac                                                     <  1:21054/149‑1 (MQ=255)
                                                                                                  gACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTg                                                         >  2:70480/1‑143 (MQ=255)
                                                                                                  gACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTg                                                         <  1:70480/143‑1 (MQ=255)
                                                                                                     ccccATTTCCTCCTTGATAAACTTGATCACCTTCTTCGCCTCGGGCGACCCGGCCTTCCACTATATACCGGCATAGATGTCTTCGGAGATCTCGCGGAAGATCACCCTGTCGACGTCGCACGGCTTTTTCCCCGGGCATGGCACACCCt                                                >  6:316116/1‑149 (MQ=255)
                                                                                                     ccccATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCTGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTc                                                                                      <  1:331674/111‑1 (MQ=255)
                                                                                                     ccccATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCTGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTc                                                                                      >  2:331674/1‑111 (MQ=255)
                                                                                                     ccccATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCt                                                <  2:203801/149‑1 (MQ=255)
                                                                                                      cccATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTg                                               <  3:265668/149‑1 (MQ=255)
                                                                                                        cATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATc                                                                                              >  5:368715/1‑100 (MQ=255)
                                                                                                        cATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATc                                                                                              <  6:368715/100‑1 (MQ=255)
                                                                                                        cATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGt                                                                                       <  4:365692/107‑1 (MQ=255)
                                                                                                        cATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGt                                                                                       >  3:365692/1‑107 (MQ=255)
                                                                                                        cATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGa                                             <  1:366797/149‑1 (MQ=255)
                                                                                                              ctcctTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACcaca                                       <  6:129961/149‑1 (MQ=255)
                                                                                                              ctcctTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACcac                                        <  2:46508/148‑1 (MQ=255)
                                                                                                               tcctTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGAATGTGCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACca                                         >  1:307645/1‑148 (MQ=255)
                                                                                                               tcctTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAg                                      >  3:125722/1‑149 (MQ=255)
                                                                                                               tcctTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACAGGGCTTGGCACACCCTGGAACCACAg                                      >  4:118982/1‑149 (MQ=255)
                                                                                                                  ttGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCAc                                   >  5:245999/1‑149 (MQ=255)
                                                                                                                   tGAGAAACTTGATCACCTTGTTTGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGcccc                                  >  7:241269/1‑149 (MQ=255)
                                                                                                                    gagaAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCg                                 >  1:5929/1‑149 (MQ=255)
                                                                                                                    gagaAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCg                                 >  4:107596/1‑149 (MQ=255)
                                                                                                                      gaAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGc                               <  4:125722/149‑1 (MQ=255)
                                                                                                                       aaaCTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGcc                              >  2:55388/1‑149 (MQ=255)
                                                                                                                              aTCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGc                        <  6:328175/148‑1 (MQ=255)
                                                                                                                                cACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGcata                     >  8:101590/1‑149 (MQ=255)
                                                                                                                                 aCCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGcata                     <  6:245999/148‑1 (MQ=255)
                                                                                                                                   cTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGc                                                                         <  6:167367/94‑1 (MQ=255)
                                                                                                                                   cTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGc                                                                         >  5:167367/1‑94 (MQ=255)
                                                                                                                                        tCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGca                       >  2:164979/1‑139 (MQ=255)
                                                                                                                                        tCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGca                       <  1:164979/139‑1 (MQ=255)
                                                                                                                                            cTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGCATACATACAGGTCCa         <  6:175108/149‑1 (MQ=255)
                                                                                                                                               gggCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGCATACATACAGGTCCAg        <  3:118982/147‑1 (MQ=255)
                                                                                                                                                 gCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGCATACATACAGGTCCAGTTGc    >  3:52493/1‑149 (MQ=255)
                                                                                                                                                   gACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGCATACATACAGGTCCAGTTGCTg  <  2:142659/149‑1 (MQ=255)
                                                                                                                                                   gACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGCATACATACAGGTCCAGTTGCTg  >  3:43888/1‑149 (MQ=255)
                                                                                                                                                   gACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGCATACATACAGGTCCAGTTGCTg  <  8:241269/149‑1 (MQ=255)
                                                                                                                                                   gACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGCATACATACAGGTCCAGTTGCTg  <  3:107596/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCACGTATTGCAGGGCCTTGCGCACCAGACGCTTGGTACCCTCCTTCGACACCGGCTTGACGCCGATACCGCAGTCCTGGTCGAAGCGGATCTTGGTGACCCCCATTTCCTCCTTGAGAAACTTGATCACCTTGTTCGCCTCGGGCGACCCGGCCTTCCACTCGATGCCGGCATAGATGTCTTCGGAGTTCTCGCGGAAGATCACCATGTCGACGTCGCCCGGCTTTTTCACCGGGCTTGGCACACCCTGGAACCACAGCACCGGCCGCAGGCATACATACAGGTCCAGTTGCTG  >  NC_002947/4521080‑4521374

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