Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,837,288 C→A 7.3% P432Q (CCA→CAA)  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,2880CA 7.3% 252.0 / 8.4 82P432Q (CCA→CAA) ccoN‑Icbb3‑type cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (40/36);  new base A (4/2);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

GCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAA  >  NC_002947/4837142‑4837436
                                                                                                                                                     |                                                                                                                                                    
gcACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGGCGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAg                                                                                                                                                    >  7:109945/1‑149 (MQ=255)
      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                                                                                                                                             <  1:29873/103‑1 (MQ=255)
                                                     tCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgt                                                                                                                                             >  2:29873/1‑103 (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                                                                                    <  5:235803/144‑1 (MQ=255)
                                                                      gggCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                                    >  6:235803/1‑144 (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                                                                   <  6:164338/149‑1 (MQ=255)
                                                                                  gcgcgcGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGc                                                                   <  7:136921/149‑1 (MQ=255)
                                                                                  gcgcgcGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGc                                                                   <  1:361613/149‑1 (MQ=255)
                                                                                     cgcgGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGtt                                                                <  4:302462/149‑1 (MQ=255)
                                                                                     cgcgGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGtt                                                                <  5:299634/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                                >  5:348315/1‑148 (MQ=255)
                                                                                                                        gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGcc                                >  1:47795/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                         >  5:361977/1‑149 (MQ=255)
                                                                                                                              ggAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg                         >  7:43224/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)
                                                                                                                                                     |                                                                                                                                                    
GCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAA  >  NC_002947/4837142‑4837436

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: