Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,840,005 A→T 6.9% Y124F (TAC→TTC)  ccoN‑II → cbb3‑type cytochrome c oxidase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,840,0050AT 6.9% 200.6 / 2.6 67Y124F (TAC→TTC) ccoN‑IIcbb3‑type cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base A (36/26);  new base T (2/3);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TTCGGCACCAGCTACTACGTGGTGCAGCGCACCTGCCAGACCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTG  >  NC_002947/4839857‑4840144
                                                                                                                                                      |                                                                                                                                             
ttCGGCACCAGCTACTACGTGGTGCAGCGCACCTTCCAGACCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATa                                                                                                                                               <  7:317333/149‑1 (MQ=255)
   gtcACCAGCTACTACGTGGTGCAGCGCACCTGCCAGTCCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATGGGCGGCGCGCTGATCACCCTGGCGATGGGCTACACCACCACCAAGGAATTCGc                                                                                                                                            <  4:82815/147‑1 (MQ=255)
   ggCACCAGCTACTACGTGGTGCAGCGCACCTGCCAGACCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGc                                                                                                                                            >  2:324843/1‑149 (MQ=255)
    gCACCAGCTACTACGTGGTGCAGCGCACCTGCCAGACCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGcc                                                                                                                                           >  5:105120/1‑149 (MQ=255)
     cACCAGCTACTACGTGGTGCAGCGCACCTGCCAGACCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGcc                                                                                                                                           >  2:328519/1‑148 (MQ=255)
     cACCAGATACTACGTGGTGCAGCGCACCTGCCAGACCCGGCTGATCTCCGACAGCATCGCCGCCTTCACCTTCTGCGGTTGGCAAGCCGTGATCGTCGGCGCGATGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGcc                                                                                                                                           >  5:378319/1‑148 (MQ=255)
      aCCAGCTACTACGTGGTGCAGCGCACCTGCCAGACCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGcc                                                                                                                                           >  7:53534/1‑147 (MQ=255)
                      tgCAGCGCACCTGCCAGACCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTg                                                                                                                        <  1:324843/149‑1 (MQ=255)
                         aGCGCACCTGCCAGACCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTAGGGCAAGCCGTGATCGTCGGCGCGCTGATAACCCTGCCGATGGGCTACCCCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGc                                                                                                                      >  1:299264/1‑149 (MQ=255)
                                  gCCAGACCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCctgctg                                                                                                            <  4:219500/149‑1 (MQ=255)
                                    cAGACCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGc                                                                                                          <  1:214383/149‑1 (MQ=255)
                                            gCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCt                                                                                                  >  7:205403/1‑149 (MQ=255)
                                            gCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCt                                                                                                  <  1:51446/149‑1 (MQ=255)
                                              tGATATCCGACAGCATGGCCGCCTTCACATTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATAAGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTgg                                                                                                >  3:293232/1‑149 (MQ=255)
                                                aTCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTggg                                                                                               <  1:328519/148‑1 (MQ=255)
                                                aTCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTgg                                                                                                <  8:205403/147‑1 (MQ=255)
                                                  ctcCGACAGAATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCCATCGGCTAAACCACAACCAAGGAATACGCCGAGCTGGCGTGGCCGCTGGCGCTCATGCTGGCCATCGTCTGGGtt                                                                                             >  7:196762/1‑148 (MQ=255)
                                                   tcCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTAc                                                                                          >  3:242970/1‑149 (MQ=255)
                                                       aCAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTAc                                                                                      >  3:73495/1‑149 (MQ=255)
                                                         aGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACgg                                                                                    >  7:86943/1‑149 (MQ=255)
                                                            aTGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTAc                                                                                          <  6:59633/140‑1 (MQ=255)
                                                            aTGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTAc                                                                                          >  5:59633/1‑140 (MQ=255)
                                                               gccgccTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGgt                                                                              >  7:361941/1‑149 (MQ=255)
                                                                       cACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGttcttc                                                                       <  2:299264/148‑1 (MQ=255)
                                                                       cACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGttcttc                                                                       <  4:242970/148‑1 (MQ=255)
                                                                        aCCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCg                                                                      <  6:105120/148‑1 (MQ=255)
                                                                           ttCTGGGGTTGGCAAGCCGTGATCGTGCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCa                                                                   >  1:196166/1‑149 (MQ=255)
                                                                             cTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCa                                                                   <  2:196166/146‑1 (MQ=255)
                                                                                gggTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTa                                                                <  4:253013/146‑1 (MQ=255)
                                                                                   ttGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTg                                                                                >  4:24180/1‑127 (MQ=255)
                                                                                      ggCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCaaa                                                        >  4:329649/1‑149 (MQ=255)
                                                                                      ggCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCaaa                                                        >  4:8624/1‑149 (MQ=255)
                                                                                       gCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTg                                                                                >  4:192642/1‑124 (MQ=255)
                                                                                       gCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTg                                                                                <  3:192642/124‑1 (MQ=255)
                                                                                       gCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAAc                                                       <  8:53534/149‑1 (MQ=255)
                                                                                          aGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATTCGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGc                                                     >  2:20657/1‑148 (MQ=255)
                                                                                          aGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGc                                                     >  5:167301/1‑148 (MQ=255)
                                                                                          aGCCGTGATCGTCGGCGCGCTGATCACACTGCCGATGGGCTAAACCACCACCAAGGAATTAGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATAGTCTGGGTTACCTACGGGATGGTGTTCTTCGGCACAATCGTCAAAccc                                                     >  2:264048/1‑146 (MQ=255)
                                                                                             cGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATTCGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGa                                                 <  1:20657/149‑1 (MQ=255)
                                                                                             cGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGa                                                 >  2:346503/1‑149 (MQ=255)
                                                                                              gTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATTCGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTAGCCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGa                                                 <  1:264048/149‑1 (MQ=255)
                                                                                              gTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGATGGTGTTCTTCGGCACTATCGTCAAACGCAAGa                                                 >  5:81807/1‑148 (MQ=255)
                                                                                                  tcgtcgGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAATGCAAGACCAAg                                            <  3:8624/149‑1 (MQ=255)
                                                                                                            cgcTGATGAGCCTGGCGATGGGCTACACCACCACCACGGAATACGCCGAGCTGGAGTGGCCGCTTGCTATCCTGCTGGCCATCTGCTGGGTTACCTACGGTCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTAt                                   <  8:196762/149‑1 (MQ=255)
                                                                                                             gcTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATg                                  <  8:86943/149‑1 (MQ=255)
                                                                                                                gATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCg                               >  3:31320/1‑149 (MQ=255)
                                                                                                                  tCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCgg                              >  2:140951/1‑148 (MQ=255)
                                                                                                                   cACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGc                             >  4:330996/1‑148 (MQ=255)
                                                                                                                     cccTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAAc                          >  3:95654/1‑149 (MQ=255)
                                                                                                                      ccTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGACGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAAc                          >  1:269594/1‑149 (MQ=255)
                                                                                                                      ccTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTc                                <  1:182791/142‑1 (MQ=255)
                                                                                                                      ccTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTc                                >  2:182791/1‑142 (MQ=255)
                                                                                                                       cTGCCGATGGGCTACACCACCACCAAGGAATAGGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTg                        <  4:73495/149‑1 (MQ=255)
                                                                                                                        tGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATg                                  >  6:338468/1‑138 (MQ=255)
                                                                                                                        tGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACAACTATGTCGGAAACTgg                       >  7:151681/1‑149 (MQ=255)
                                                                                                                           cGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGc                             <  4:146078/140‑1 (MQ=255)
                                                                                                                           cGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGc                             >  3:146078/1‑140 (MQ=255)
                                                                                                                                 gCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGt              >  7:79573/1‑149 (MQ=255)
                                                                                                                                  cTACACCACCACCAAGGAATACGCCGAGCTGGAGTGTCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCTGCACTATCGTCAAACGCAAGACCAAGAACATCTATGTc                                >  7:281658/1‑130 (MQ=255)
                                                                                                                                     caccaccacCAATGAATACGCCTAGCTGGAGTGGCCGAGGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGt          <  8:151681/149‑1 (MQ=255)
                                                                                                                                     caccaccacCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGt          <  4:293232/149‑1 (MQ=255)
                                                                                                                                     caccaccacCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGt          >  8:147096/1‑149 (MQ=255)
                                                                                                                                       ccaccaccaAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTc        <  1:140951/149‑1 (MQ=255)
                                                                                                                                        caccaccaAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCa       >  8:323402/1‑149 (MQ=255)
                                                                                                                                          ccaccaAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATc     <  8:79573/149‑1 (MQ=255)
                                                                                                                                             ccaAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTg  <  8:361941/149‑1 (MQ=255)
                                                                                                                                                  gAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATc                                                                                                      <  5:61323/45‑1 (MQ=255)
                                                                                                                                                  gAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATc                                                                                                      >  6:61323/1‑45 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
TTCGGCACCAGCTACTACGTGGTGCAGCGCACCTGCCAGACCCGGCTGATCTCCGACAGCATGGCCGCCTTCACCTTCTGGGGTTGGCAAGCCGTGATCGTCGGCGCGCTGATCACCCTGCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTG  >  NC_002947/4839857‑4840144

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: