Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 109,430 C→T 8.9% V207V (GTC→GTT)  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,4300CT 8.9% 190.8 / 6.4 67V207V (GTC→GTT) PP_0103cytochrome c oxidase subunit 2
Reads supporting (aligned to +/- strand):  ref base C (34/27);  new base T (3/3);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

TGAGCACTACCTGCTGGAAGTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAGCGCAAGG  >  NC_002947/109286‑109575
                                                                                                                                                 |                                                                                                                                                  
tGAGCACTACCTGCTGGAAGTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACg                                                                                                                                                <  4:27891/149‑1 (MQ=255)
         ccagcgggAAGTCGACCACCCGCTGGTGCGGCCGGTGGGCGCGAAAGTGCGCTGCCTGGTCACAGCTGCCGAGGTGAGCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCACAC‑GCGACGCCATCCCAGGCTTCGTTGACAAAGCCTgg                                                                                                                                        <  6:355974/142‑1 (MQ=255)
             ctgGAAGTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCg                                                                                                                                  >  1:66616/1‑149 (MQ=255)
                gAAGTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGAt                                                                                                                               >  2:72631/1‑149 (MQ=255)
                   gTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCg                                                                                                                             >  8:170591/1‑148 (MQ=255)
                      gACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGaaaa                                                                                                                         <  7:170591/149‑1 (MQ=255)
                      gACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTCGTTGACAAAGCCTg                                                                                                                                         >  3:191419/1‑134 (MQ=255)
                      gACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTCGTTGACAAAGCCTg                                                                                                                                         <  4:191419/134‑1 (MQ=255)
                        ccAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGaaaa                                                                                                                         <  8:223761/147‑1 (MQ=255)
                        ccAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGc                                                                                                                       <  2:35251/149‑1 (MQ=255)
                          aGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGc                                                                                                                       >  4:315987/1‑147 (MQ=255)
                           gCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAg                                                                                                                    <  4:24455/149‑1 (MQ=255)
                                  gTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTa                                                                                                             >  5:210772/1‑149 (MQ=255)
                                    gctgcCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTCGTTGACAAAGCCTg                                                                                                                                         <  4:650/120‑1 (MQ=255)
                                    gctgcCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTCGTTGACAAAGCCTg                                                                                                                                         >  3:650/1‑120 (MQ=255)
                                      tgcCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACcgc                                                                                                         >  6:330577/1‑149 (MQ=255)
                                           gTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCa                                                                                                    <  6:210772/149‑1 (MQ=255)
                                             gggCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGc                                                                                                      <  5:289765/145‑1 (MQ=255)
                                             gggCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGc                                                                                                      >  6:289765/1‑145 (MQ=255)
                                               gcgcGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGCCGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGAGCGAAAAGCCAGGCATCTACCGCGGCCAGTGc                                                                                                >  8:115736/1‑149 (MQ=255)
                                                   gAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCAcc                                                                                             >  1:152210/1‑148 (MQ=255)
                                                        tGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGACCCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCt                                                                                        <  1:72631/149‑1 (MQ=255)
                                                         cggcTTCCTGGGCCCTGCTGCGTATGTTTTCAATTCATGGGGGGTGCGGGCCTGCGGGTGAAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCtgt                                                                                      <  2:66616/147‑1 (MQ=255)
                                                          ctcTTGCTGGTCACTGCTGCGTATTTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGGCATCCCAGGCTTTGTCAACGTAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCtgt                                                                                      <  2:152210/146‑1 (MQ=255)
                                                             ttCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCgg                                                                                  >  6:49898/1‑149 (MQ=255)
                                                                 tGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAg                                                                              >  1:272258/1‑149 (MQ=255)
                                                                   gTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGa                                                                            <  5:293496/149‑1 (MQ=255)
                                                                       ctgctgCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTAAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGACAGTGCACCGAGCTGTGCGGCAAGGAcaaa                                                                        >  8:135049/1‑146 (MQ=255)
                                                                        ttatgcCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCACATGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACg                                                                       <  7:135049/146‑1 (MQ=255)
                                                                        tgctgcCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACg                                                                       >  5:235801/1‑149 (MQ=255)
                                                                           tgcCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCt                                                                    <  6:235801/149‑1 (MQ=255)
                                                                               gATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCAt                                                                >  4:164614/1‑149 (MQ=255)
                                                                                     aTCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCggt                                                          >  6:269941/1‑149 (MQ=255)
                                                                                      tCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCggtg                                                         <  3:164614/149‑1 (MQ=255)
                                                                                        gatctcTGGTGGGGTCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGTCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCggtgg                                                        <  8:331160/143‑1 (MQ=255)
                                                                                         aCTCCTGGTGGGTGCCGGCCTTCGCGGTCAAACGGCGACGCCATCCAAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCggtggt                                                       >  1:222977/1‑149 (MQ=255)
                                                                                           tCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCggtggtggt                                                    >  1:135124/1‑149 (MQ=255)
                                                                                              tggtggGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGTA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGATTCATGCCcgtgggagatc                                                  >  7:331160/1‑142 (MQ=255)
                                                                                              tggtggGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTc                                                   >  6:302894/1‑147 (MQ=255)
                                                                                              tggtggGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCACTGCACCGCGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTTGTc                                                   >  8:79425/1‑147 (MQ=255)
                                                                                                   ggtgCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCAcc                                                                                             <  8:254087/100‑1 (MQ=255)
                                                                                                   ggtgCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCAcc                                                                                             >  7:254087/1‑100 (MQ=255)
                                                                                                    gtgCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCaa                                           <  2:222977/149‑1 (MQ=255)
                                                                                                         ggCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCa                                      <  5:269941/149‑1 (MQ=255)
                                                                                                             ttCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCaa                                                                               >  1:197312/1‑104 (MQ=255)
                                                                                                             ttCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCaa                                                                               <  2:197312/104‑1 (MQ=255)
                                                                                                             ttCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGc                                  >  6:21087/1‑149 (MQ=255)
                                                                                                             ttCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGc                                  >  6:300122/1‑149 (MQ=255)
                                                                                                             ttCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGc                                  >  4:150480/1‑149 (MQ=255)
                                                                                                               cgcgGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCATCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGc                                  <  3:315987/147‑1 (MQ=255)
                                                                                                               cgcgGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCg                                 >  1:56899/1‑148 (MQ=255)
                                                                                                                gcgGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGa                               <  3:150480/149‑1 (MQ=255)
                                                                                                                  ggTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGAtt                             <  2:56899/149‑1 (MQ=255)
                                                                                                                      aaaC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACg                          >  2:265581/1‑148 (MQ=255)
                                                                                                                      aaaC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGa                         >  5:117003/1‑149 (MQ=255)
                                                                                                                        aC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGaca                       >  3:156075/1‑149 (MQ=255)
                                                                                                                        aC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGaca                       >  8:179021/1‑149 (MQ=255)
                                                                                                                           gcgACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACAcc                     >  4:287722/1‑149 (MQ=255)
                                                                                                                                 ccATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCtt               <  5:330577/149‑1 (MQ=255)
                                                                                                                                   ggtccAGGCTTCGTTGACAAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTgg             >  3:182857/4‑149 (MQ=255)
                                                                                                                                      ccAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGc            >  4:175490/1‑147 (MQ=255)
                                                                                                                                        aGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAgc        >  6:152728/1‑149 (MQ=255)
                                                                                                                                          gCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGAGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAgcg       >  5:237692/1‑148 (MQ=255)
                                                                                                                                              tGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAGCGCaa    <  5:300122/147‑1 (MQ=255)
                                                                                                                                              tGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAGCGCaa    <  5:21087/147‑1 (MQ=255)
                                                                                                                                               gTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAGCGCAAgg  <  1:265581/148‑1 (MQ=255)
                                                                                                                                               gTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAGCGCAAgg  <  3:287722/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGAGCACTACCTGCTGGAAGTCGACCAGCCGCTGGTGCTGCCGGTGGGCGCGAAAGTGCGCTTCCTGGTCACTGCTGCCGATGTGATCCACTCCTGGTGGGTGCCGGCCTTCGCGGTCAAAC‑GCGACGCCATCCCAGGCTTTGTCAACGAAGCCTGGA‑CCCGGATCGAAAAGCCAGGCATCTACCGCGGCCAGTGCACCGAGCTGTGCGGCAAGGACCACGGCTTCATGCCGGTGGTGGTCGAGGTCAAGTCCAAGGCGGATTACGACACCTGGCTTGGCGAGCGCAAGG  >  NC_002947/109286‑109575

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: ‑