Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,837,281 A→G 8.1% S430G (AGC→GGC)  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,2810AG 8.1% 211.1 / 6.5 74S430G (AGC→GGC) ccoN‑Icbb3‑type cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base A (36/32);  new base G (4/2);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

ATCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCC  >  NC_002947/4837134‑4837429
                                                                                                                                                      |                                                                                                                                                    
aTCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGgccag                                                                                                                                                     >  7:72834/1‑149 (MQ=255)
        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                                                                          <  8:255086/146‑1 (MQ=255)
                                                                               ggCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAacac                                                                          >  7:255086/1‑146 (MQ=255)
                                                                                          gcgcgcGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGc                                                            <  7:136921/149‑1 (MQ=255)
                                                                                          gcgcgcGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGc                                                            <  1:361613/149‑1 (MQ=255)
                                                                                          gcgcgcGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGc                                                            <  7:281836/149‑1 (MQ=255)
                                                                                          gcgcgcGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGc                                                            <  6:164338/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                                                                     <  4:56321/125‑1 (MQ=255)
                                                                                                         cGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                                     >  3:56321/1‑125 (MQ=255)
                                                                                                                        aCCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGc                                >  4:20968/1‑149 (MQ=255)
                                                                                                                             cTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACg                           >  7:323973/1‑149 (MQ=255)
                                                                                                                                gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGcc                         >  5:348315/1‑148 (MQ=255)
                                                                                                                                gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGcc                         >  1:47795/1‑148 (MQ=255)
                                                                                                                                gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGcc                         >  8:118487/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                  >  1:171060/1‑149 (MQ=255)
                                                                                                                                      ggAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg                  >  5:361977/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               <  7:118487/149‑1 (MQ=255)
                                                                                                                                         aaCCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt               >  1:199897/1‑149 (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                >  3:61898/1‑140 (MQ=255)
                                                                                                                                                 tggCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCt                <  4:61898/140‑1 (MQ=255)
                                                                                                                                                   gccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGtt                    <  5:152810/134‑1 (MQ=255)
                                                                                                                                                   gccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGtt                    >  6:152810/1‑134 (MQ=255)
                                                                                                                                                   gccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAAc     >  8:205308/1‑149 (MQ=255)
                                                                                                                                                     cagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGc   <  3:357669/149‑1 (MQ=255)
                                                                                                                                                      agccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGAAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGcc  >  4:367650/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ATCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCA‑‑CGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCC  >  NC_002947/4837134‑4837429

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: