Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 109,427 T→C 9.5% F206F (TTT→TTC PP_0103 → cytochrome c oxidase subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947109,4270TC 9.5% 166.7 / 7.0 63F206F (TTT→TTCPP_0103cytochrome c oxidase subunit 2
Reads supporting (aligned to +/- strand):  ref base T (34/23);  new base C (3/3);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CAAAGGATGAGCACTACCTGCTGGAAGTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAGCGCAA  >  NC_002947/109279‑109573
                                                                                                                                                     |                                                                                                                                                   
cAAAGGATGAGCACTACCTGCTGGAAGTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCttt                                                                                                                                                     <  4:135064/149‑1 (MQ=255)
       tGAGCACTACCTGCTGGAAGTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACg                                                                                                                                              <  4:27891/149‑1 (MQ=255)
                ccagcgggAAGTCGACCACCCGCTGGTGCGGCCGGTGGGCGCGAAAGTGCGCTGCCTGGTCACAGCTGCCGAGGTGAGCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCACACGCGACGCCATCCCAGGCTTCGTTGACAAAGCCTgg                                                                                                                                      <  6:355974/142‑1 (MQ=255)
                    ctgGAAGTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCg                                                                                                                                >  1:66616/1‑149 (MQ=255)
                       gAAGTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGAt                                                                                                                             >  2:72631/1‑149 (MQ=255)
                          gTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCg                                                                                                                           >  8:170591/1‑148 (MQ=255)
                             gACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGaaaa                                                                                                                       <  7:170591/149‑1 (MQ=255)
                             gACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTCGTTGACAAAGCCTg                                                                                                                                       >  3:191419/1‑134 (MQ=255)
                             gACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTCGTTGACAAAGCCTg                                                                                                                                       <  4:191419/134‑1 (MQ=255)
                               ccAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGaaaa                                                                                                                       <  8:223761/147‑1 (MQ=255)
                               ccAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGc                                                                                                                     <  2:35251/149‑1 (MQ=255)
                                 aGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGc                                                                                                                     >  4:315987/1‑147 (MQ=255)
                                  gCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAg                                                                                                                  <  4:24455/149‑1 (MQ=255)
                                         gTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTa                                                                                                           >  5:210772/1‑149 (MQ=255)
                                           gctgcCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTCGTTGACAAAGCCTg                                                                                                                                       >  3:650/1‑120 (MQ=255)
                                           gctgcCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTCGTTGACAAAGCCTg                                                                                                                                       <  4:650/120‑1 (MQ=255)
                                             tgcCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACcgc                                                                                                       >  6:330577/1‑149 (MQ=255)
                                                  gTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCa                                                                                                  <  6:210772/149‑1 (MQ=255)
                                                    gggCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGc                                                                                                    >  6:289765/1‑145 (MQ=255)
                                                    gggCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGc                                                                                                    <  5:289765/145‑1 (MQ=255)
                                                      gcgcGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGCCGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGAGCGAAAAGCCAGGCATCTACCGCGGCCAGTGc                                                                                              >  8:115736/1‑149 (MQ=255)
                                                          gAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCAcc                                                                                           >  1:152210/1‑148 (MQ=255)
                                                               tGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCt                                                                                      <  1:72631/149‑1 (MQ=255)
                                                                cggcTTCCTGGGCCCTGCTGCGTATGTTTTCAATTCATGGGGGGTGCGGGCCTGCGGGTGAAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCtgt                                                                                    <  2:66616/147‑1 (MQ=255)
                                                                 ctcTTGCTGGTCACTGCTGCGTATTTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGGCATCCCAGGCTTTGTCAACGTAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCtgt                                                                                    <  2:152210/146‑1 (MQ=255)
                                                                    ttCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCgg                                                                                >  6:49898/1‑149 (MQ=255)
                                                                        tGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAg                                                                            >  1:272258/1‑149 (MQ=255)
                                                                          gTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGa                                                                          <  5:293496/149‑1 (MQ=255)
                                                                              ctgctgCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTAAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGACAGTGCACCGAGCTGTGCGGCAAGGAcaaa                                                                      >  8:135049/1‑146 (MQ=255)
                                                                               ttatgcCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCACATGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACg                                                                     <  7:135049/146‑1 (MQ=255)
                                                                               tgctgcCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACg                                                                     >  5:235801/1‑149 (MQ=255)
                                                                                  tgcCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCt                                                                  <  6:235801/149‑1 (MQ=255)
                                                                                      gATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCAt                                                              >  4:164614/1‑149 (MQ=255)
                                                                                            aTCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCggt                                                        >  6:269941/1‑149 (MQ=255)
                                                                                             tCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCggtg                                                       <  3:164614/149‑1 (MQ=255)
                                                                                               gatctcTGGTGGGGTCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGTCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCggtgg                                                      <  8:331160/143‑1 (MQ=255)
                                                                                                aCTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGGCGACGCCATCCAAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCggtggt                                                     >  1:222977/1‑149 (MQ=255)
                                                                                                  tCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCggtggtggt                                                  >  1:135124/1‑149 (MQ=255)
                                                                                                     tggtggGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGTACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGATTCATGCCcgtgggagatc                                                >  7:331160/1‑142 (MQ=255)
                                                                                                     tggtggGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTc                                                 >  6:302894/1‑147 (MQ=255)
                                                                                                     tggtggGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCACTGCACCGCGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTTGTc                                                 >  8:79425/1‑147 (MQ=255)
                                                                                                          ggtgCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCAcc                                                                                           <  8:254087/100‑1 (MQ=255)
                                                                                                          ggtgCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCAcc                                                                                           >  7:254087/1‑100 (MQ=255)
                                                                                                           gtgCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCaa                                         <  2:222977/149‑1 (MQ=255)
                                                                                                                ggCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCa                                    <  5:269941/149‑1 (MQ=255)
                                                                                                                    ttCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCaa                                                                             >  1:197312/1‑104 (MQ=255)
                                                                                                                    ttCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCaa                                                                             <  2:197312/104‑1 (MQ=255)
                                                                                                                    ttCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGc                                >  6:21087/1‑149 (MQ=255)
                                                                                                                    ttCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGc                                >  6:300122/1‑149 (MQ=255)
                                                                                                                    ttCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGc                                >  4:150480/1‑149 (MQ=255)
                                                                                                                      cgcgGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCATCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGc                                <  3:315987/147‑1 (MQ=255)
                                                                                                                      cgcgGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCg                               >  1:56899/1‑148 (MQ=255)
                                                                                                                       gcgGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGa                             <  3:150480/149‑1 (MQ=255)
                                                                                                                         ggTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGAtt                           <  2:56899/149‑1 (MQ=255)
                                                                                                                             aaaCGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACg                        >  2:265581/1‑148 (MQ=255)
                                                                                                                             aaaCGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGa                       >  5:117003/1‑149 (MQ=255)
                                                                                                                               aCGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGaca                     >  3:156075/1‑149 (MQ=255)
                                                                                                                               aCGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGaca                     >  8:179021/1‑149 (MQ=255)
                                                                                                                                  gcgACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACAcc                   >  4:287722/1‑149 (MQ=255)
                                                                                                                                        ccATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCtt             <  5:330577/149‑1 (MQ=255)
                                                                                                                                          ggtccAGGCTTCGTTGACAAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTgg           >  3:182857/4‑149 (MQ=255)
                                                                                                                                             ccAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGc          >  4:175490/1‑147 (MQ=255)
                                                                                                                                               aGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAgc      >  6:152728/1‑149 (MQ=255)
                                                                                                                                                 gCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGAGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAgcg     >  5:237692/1‑148 (MQ=255)
                                                                                                                                                     tGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAGCGCaa  <  5:300122/147‑1 (MQ=255)
                                                                                                                                                     tGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAGCGCaa  <  5:21087/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAAAGGATGAGCACTACCTGCTGGAAGTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAGCGCAA  >  NC_002947/109279‑109573

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: