Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,251,595 T→C 6.2% M83V (ATG→GTG)  PP_1094 ← electron transport complex protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,251,5950TC 6.2% 240.2 / 4.3 80M83V (ATG→GTG) PP_1094electron transport complex protein
Reads supporting (aligned to +/- strand):  ref base T (36/39);  new base C (2/3);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

GATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGCCGGC  >  NC_002947/1251448‑1251742
                                                                                                                                                    |                                                                                                                                                   
gATCGGCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCAGGATCGGCCCGGGATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCAGCCATa                                                                                                                                                    <  8:163981/149‑1 (MQ=255)
 aTCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAg                                                                                                                                                   <  2:368170/149‑1 (MQ=255)
  tCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGt                                                                                                                                                  <  6:376608/149‑1 (MQ=255)
                 cTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCGTGCCACTATGGCAGGTAAgccaa                                                                                                                                   <  2:331647/149‑3 (MQ=255)
                    tCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCa                                                                                                                                 <  1:299083/148‑1 (MQ=255)
                    tCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCa                                                                                                                                 <  2:280764/148‑1 (MQ=255)
                      cTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGa                                                                                                                              <  3:269505/149‑1 (MQ=255)
                       tgttgtTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGa                                                                                                                              >  2:321103/1‑148 (MQ=255)
                        gttgttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGAGAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATa                                                                                                                            >  6:272513/1‑149 (MQ=255)
                         ttgttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCGTGCCACTATGGCAGGTAAgccaa                                                                                                                                   <  5:338872/141‑3 (MQ=255)
                         ttgttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCGTGCCACTATGGCAGGTAAgccaa                                                                                                                                   >  6:338872/1‑139 (MQ=255)
                           gttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACg                                                                                                                          <  2:68896/148‑1 (MQ=255)
                            ttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGc                                                                                                                        >  3:109042/1‑149 (MQ=255)
                             tCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCGTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGc                                                                                                                        >  2:243139/1‑148 (MQ=255)
                                 cgaCCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTg                                                                                                                   >  6:206704/3‑149 (MQ=255)
                                    ccTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGc                                                                                                                >  4:93310/1‑149 (MQ=255)
                                    ccTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGc                                                                                                                <  4:358000/149‑1 (MQ=255)
                                    ccTGGCTGCCGATCAGCCGGCCCTCCGTGGTGATGGCGACGCGCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCACGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGc                                                                                                                <  2:137724/149‑1 (MQ=255)
                                      tGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCct                                                                                                              >  4:86904/1‑149 (MQ=255)
                                        gCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCctct                                                                                                            >  3:46318/1‑149 (MQ=255)
                                          tGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCtgt                                                                                                          >  2:43343/1‑149 (MQ=255)
                                            ccGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGAt                                                                                                                             >  6:363106/1‑128 (MQ=255)
                                            ccGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGAt                                                                                                                             <  5:363106/128‑1 (MQ=255)
                                               aTCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTgcggc                                                                                                     >  3:365324/1‑149 (MQ=255)
                                                   gccggccgTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGc                                                                                                 >  3:78880/1‑149 (MQ=255)
                                                       gccgTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTg                                                                                             >  5:207691/1‑149 (MQ=255)
                                                        ccgtccgtGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTCgg                                                                                            <  4:365324/149‑1 (MQ=255)
                                                             cgtGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCg                                                                                       >  5:236064/1‑149 (MQ=255)
                                                             cgtGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCg                                                                                       >  5:241469/1‑149 (MQ=255)
                                                              gtggtgATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCgg                                                                                      >  3:196135/1‑149 (MQ=255)
                                                              gtggtgATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCgg                                                                                      <  4:235921/149‑1 (MQ=255)
                                                                 gtgATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCt                                                                                   <  3:56298/149‑1 (MQ=255)
                                                                  tgATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCtt                                                                                  <  2:234970/149‑1 (MQ=255)
                                                                    aTGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTc                                                                                <  2:372409/149‑1 (MQ=255)
                                                                       gCGACGCTCAGCATGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCa                                                                              <  1:321103/148‑1 (MQ=255)
                                                                         gACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCCCATTGCATACCTGCCATAGTGGCACGGTAAGCCGCCAGGCTTCCGCGAATTGGCCGATGTGCGGGGGGCCTTGGCCGGGGGCTTTACAtgg                                                                           >  1:98566/1‑149 (MQ=255)
                                                                             cTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTg                                                                       <  7:42522/149‑1 (MQ=255)
                                                                                 gCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTg                                                                   <  4:78880/149‑1 (MQ=255)
                                                                                  cACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTg                                                                   >  5:137820/1‑148 (MQ=255)
                                                                                  cACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGt                                                                  <  4:196135/149‑1 (MQ=255)
                                                                                  cACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGt                                                                  >  8:347755/1‑149 (MQ=255)
                                                                                    cGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCa                                                                <  2:111059/149‑1 (MQ=255)
                                                                                    cGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCa                                                                <  4:46318/149‑1 (MQ=255)
                                                                                     gATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCAt                                                               <  2:170924/149‑1 (MQ=255)
                                                                                     gATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCAt                                                               <  6:137820/149‑1 (MQ=255)
                                                                                     gATCGGCCCGGCATACCCCTGGGCCTGTGTGAGCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCa                                                                >  5:118835/1‑148 (MQ=255)
                                                                                              ggCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGt                                                      >  3:266033/1‑149 (MQ=255)
                                                                                               gCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATa                                                              >  1:244662/1‑140 (MQ=255)
                                                                                               gCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATa                                                              <  2:244662/140‑1 (MQ=255)
                                                                                                 aTACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCgg                                                   <  4:109042/149‑1 (MQ=255)
                                                                                                  tACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCggg                                                  <  1:5934/149‑1 (MQ=255)
                                                                                                    ccccTGCGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGgca                                                <  2:254186/149‑1 (MQ=255)
                                                                                                     cccTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATa                                                              <  6:286151/134‑1 (MQ=255)
                                                                                                     cccTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATa                                                              >  5:286151/1‑134 (MQ=255)
                                                                                                     cccTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGgcag                                               >  6:193694/1‑149 (MQ=255)
                                                                                                       cTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCGTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGt                                                      <  4:252782/140‑1 (MQ=255)
                                                                                                       cTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGt                                                      >  3:252782/1‑140 (MQ=255)
                                                                                                            ccTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGcc                                        <  5:378699/149‑1 (MQ=255)
                                                                                                                   gATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCt                                 >  1:25996/1‑149 (MQ=255)
                                                                                                                    aTCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACCGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTg                                >  6:279138/1‑149 (MQ=255)
                                                                                                                      cAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGgc                                                 >  6:92007/1‑130 (MQ=255)
                                                                                                                      cAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTggc                              >  1:368971/1‑149 (MQ=255)
                                                                                                                       aGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTgtcg                             <  4:266033/149‑1 (MQ=255)
                                                                                                                          aTAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTACCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCt                          >  8:248217/1‑149 (MQ=255)
                                                                                                                           tAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCtt                         >  3:201785/1‑149 (MQ=255)
                                                                                                                             aTGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTg                       <  2:98566/149‑1 (MQ=255)
                                                                                                                                   ggcggCTTACCGTGCCACTATGGCAGGTAAGCCACCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAgg                                                                            >  6:295421/4‑90 (MQ=21)
                                                                                                                                   ggcggCTTACCGTGCCACTATGGCAGGTAAGCCACCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAgg                                                                            <  5:295421/87‑1 (MQ=21)
                                                                                                                                     tGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCTTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAgctgc                <  8:130715/148‑1 (MQ=255)
                                                                                                                                      ggCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCa                     >  6:11494/1‑142 (MQ=255)
                                                                                                                                      ggCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCa                     <  5:11494/142‑1 (MQ=255)
                                                                                                                                           aCCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCa          <  3:93310/148‑1 (MQ=255)
                                                                                                                                             cctctcATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAgc        <  7:344389/144‑1 (MQ=255)
                                                                                                                                               tGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAgcgc      >  4:1575/1‑149 (MQ=255)
                                                                                                                                                gcaaTAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTTCAGCTGCTTTTCAGCGcc     >  2:89911/4‑149 (MQ=255)
                                                                                                                                                gCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGcc     >  4:139890/1‑149 (MQ=255)
                                                                                                                                                gCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGcc     >  6:44089/1‑149 (MQ=255)
                                                                                                                                                  cATATTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGCCgg   <  7:152269/149‑1 (MQ=255)
                                                                                                                                                  cATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGCCgg   <  5:206704/149‑1 (MQ=255)
                                                                                                                                                   aTAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGCCGGc  <  3:139890/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGCCGGC  >  NC_002947/1251448‑1251742

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: