Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,837,292 C→T 7.0% G433G (GGC→GGT ccoN‑I → 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,837,2920CT 7.0% 253.9 / 6.8 86G433G (GGC→GGTccoN‑Icbb3‑type cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (40/39);  new base T (4/2);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

CTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAA  >  NC_002947/4837148‑4837437
                                                                                                                                                   |                                                                                                                                                 
ctggctggCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtc                                                                                                                                                <  5:312098/148‑1 (MQ=255)
   gctggcCACCATCGGTACCGTGCTGTACATCGCCTCTATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgtc                                                                                                                                             <  5:141215/148‑1 (MQ=255)
   gctggcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgtcg                                                                                                                                            <  1:225464/149‑1 (MQ=255)
   gctggcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgtc                                                                                                                                             <  7:206393/148‑1 (MQ=255)
         caccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCt                                                                                                                                      >  5:197158/1‑149 (MQ=255)
          accaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCtt                                                                                                                                     >  8:281836/1‑149 (MQ=255)
           ccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTtc                                                                                                                                    <  8:142734/149‑1 (MQ=255)
           ccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTtc                                                                                                                                    <  2:120379/149‑1 (MQ=255)
           ccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTGGGTGGCTGGCCACCAAGGTTtccacgaac                                                                                                                                         >  1:325257/1‑137 (MQ=18)
           ccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTGGGTGGCTGGCCACCAAGGTTtccacgaac                                                                                                                                         <  2:325257/144‑8 (MQ=18)
                 gTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCACAGGGCATGATGTGGCGCGCGGTCAACAGCGCCGGCA‑‑CGATCACATACTCGTTCTTGGAAACATTGGTGGCCAGCCACCCCGGCTACGTCGTGCGGTTCGTggg                                                                                                                               >  1:77898/1‑149 (MQ=18)
                           gTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCt                                                                                                                    >  2:366978/1‑149 (MQ=255)
                              cATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGACTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTc                                                                                                                  <  3:171865/148‑1 (MQ=255)
                              cATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTCGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTcc                                                                                                                 >  6:299634/1‑149 (MQ=255)
                               aTCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCt                                                                                                                >  3:302462/1‑149 (MQ=255)
                                  gCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCa                                                                                                              <  3:347496/148‑1 (MQ=255)
                                   ccTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTGCTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGc                                                                                                            >  1:364129/1‑149 (MQ=255)
                                    cTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCg                                                                                                           >  5:164338/1‑149 (MQ=255)
                                     tCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCAAGGGATCACATACACGTTCTTGGAAACATTGGTGCCCAGCCACCCAGGATACGACGTCCGATTCGTCGGCGGCGCGAAATTCCTCAGc                                                                                                            >  1:106626/1‑149 (MQ=21)
                                     tCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCgg                                                                                                          >  6:198340/1‑149 (MQ=255)
                                       gATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGc                                                                                                         <  1:366978/148‑1 (MQ=255)
                                        aTGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCAt                                                                                                       >  7:348477/1‑149 (MQ=255)
                                           tgGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGt                                                                                                     <  6:50174/148‑1 (MQ=255)
                                               tCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCTTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCAATCTTCCTAAGCGGCATGTTCCt                                                                                                 >  4:254554/1‑148 (MQ=255)
                                               tCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgt                                                                                                                                              >  2:29873/1‑103 (MQ=255)
                                               tCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgt                                                                                                                                              <  1:29873/103‑1 (MQ=255)
                                                 aaCGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGGCGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGAt                                                                                              <  8:109945/149‑1 (MQ=255)
                                                  aCGGGATCTCGCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAAGCTTTGTGGCCAGCCAGCCAGGCTTCGTCGTTCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGAt                                                                                              <  3:151008/148‑1 (MQ=255)
                                                      cATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTTCTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCt                                                                                         <  2:364129/149‑1 (MQ=255)
                                                      cATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCt                                                                                         >  6:41085/1‑149 (MQ=255)
                                                       aTCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTg                                                                                        <  2:106626/149‑1 (MQ=255)
                                                          gCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTACTGATGGCCTGGac                                                                                     >  2:55687/1‑148 (MQ=255)
                                                            ccAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAaca                                                                                   >  3:177181/1‑149 (MQ=255)
                                                                gggCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                                     >  6:235803/1‑144 (MQ=255)
                                                                gggCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                                     <  5:235803/144‑1 (MQ=255)
                                                                 ggCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAacac                                                                                  >  7:255086/1‑146 (MQ=255)
                                                                 ggCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAacac                                                                                  <  8:255086/146‑1 (MQ=255)
                                                                            gcgcgcGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGc                                                                    <  7:281836/149‑1 (MQ=255)
                                                                            gcgcgcGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGc                                                                    <  1:361613/149‑1 (MQ=255)
                                                                            gcgcgcGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGc                                                                    <  6:164338/149‑1 (MQ=255)
                                                                            gcgcgcGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGc                                                                    <  7:136921/149‑1 (MQ=255)
                                                                               cgcgGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGtt                                                                 <  5:299634/149‑1 (MQ=255)
                                                                               cgcgGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGtt                                                                 <  4:302462/149‑1 (MQ=255)
                                                                                      aaCAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAgc                                                          >  2:129089/1‑149 (MQ=255)
                                                                                        cAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAgcgc                                                        >  7:33872/1‑149 (MQ=255)
                                                                                           cGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                                             >  3:56321/1‑125 (MQ=255)
                                                                                           cGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                                             <  4:56321/125‑1 (MQ=255)
                                                                                                          aCCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGc                                        >  4:20968/1‑149 (MQ=255)
                                                                                                               cTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACg                                   >  7:323973/1‑149 (MQ=255)
                                                                                                                  gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGcc                                 >  8:118487/1‑148 (MQ=255)
                                                                                                                  gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGcc                                 >  1:47795/1‑148 (MQ=255)
                                                                                                                  gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGcc                                 >  5:348315/1‑148 (MQ=255)
                                                                                                                  gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGc                                  <  2:77898/147‑1 (MQ=255)
                                                                                                                   tgcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGACGTGCGCTTCGTCGCCGGCGCGATCTTCCTCAGCCGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCa                               <  8:72834/147‑1 (MQ=255)
                                                                                                                      gtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGtt                            <  8:348477/149‑1 (MQ=255)
                                                                                                                      gacgAAGCCTGGGTGGCTGGCCACCAAGGTTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGtt                            >  4:297083/4‑149 (MQ=255)
                                                                                                                        ggAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg                          >  7:43224/1‑149 (MQ=255)
                                                                                                                        ggAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg                          >  5:361977/1‑149 (MQ=255)
                                                                                                                        ggAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg                          >  1:171060/1‑149 (MQ=255)
                                                                                                                         gAAGCCTGGGTGGCTGGCCACCAAGGTTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGc                         >  6:235388/1‑149 (MQ=255)
                                                                                                                         gAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGc                         <  4:177181/149‑1 (MQ=255)
                                                                                                                           aaCCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCTGTTGGCtt                       <  2:199897/149‑1 (MQ=255)
                                                                                                                           aaCCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt                       >  1:199897/1‑149 (MQ=255)
                                                                                                                           aaCCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt                       <  7:118487/149‑1 (MQ=255)
                                                                                                                            aCCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt                       >  1:358623/1‑148 (MQ=255)
                                                                                                                              ctgggtggCTGGCCACCAAGGTTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgag                    <  3:297083/146‑1 (MQ=255)
                                                                                                                              cTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgag                    >  8:95614/1‑149 (MQ=255)
                                                                                                                              cTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTCGCGCACCCTGACTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTagag                    >  4:246916/1‑149 (MQ=255)
                                                                                                                               tgggtggCTGGCCACCAAGGTTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGACCACCTGGCGCCCCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgag                    >  8:140261/3‑148 (MQ=255)
                                                                                                                                 ggtggCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCAATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgagg                   <  3:254554/147‑1 (MQ=255)
                                                                                                                                  gtggCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgaggag                 >  4:357669/1‑148 (MQ=255)
                                                                                                                                   tggCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCt                        <  4:61898/140‑1 (MQ=255)
                                                                                                                                   tggCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCt                        >  3:61898/1‑140 (MQ=255)
                                                                                                                                     gccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGtt                            >  6:152810/1‑134 (MQ=255)
                                                                                                                                     gccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGtt                            <  5:152810/134‑1 (MQ=255)
                                                                                                                                     gccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAAc             >  8:205308/1‑149 (MQ=255)
                                                                                                                                       cagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGc           <  3:357669/149‑1 (MQ=255)
                                                                                                                                        agccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGAAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGcc          >  4:367650/1‑149 (MQ=255)
                                                                                                                                         gccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCt         <  1:129089/149‑1 (MQ=255)
                                                                                                                                         gccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCAATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCt         >  4:19986/1‑149 (MQ=255)
                                                                                                                                            aCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCtgat      <  6:348315/149‑1 (MQ=255)
                                                                                                                                              ccAGGCTTCGTCGTGCGCTTCGCCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCtgatg     >  6:240859/1‑148 (MQ=255)
                                                                                                                                               cAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCAATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGaa   <  3:19986/149‑1 (MQ=255)
                                                                                                                                                aGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATTGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGaaa  <  3:246916/149‑1 (MQ=255)
                                                                                                                                                aGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGaaa  <  7:283473/149‑1 (MQ=255)
                                                                                                                                                aGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGaaa  >  8:283473/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAA  >  NC_002947/4837148‑4837437

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: