Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,464 A→C 6.2% L66R (CTG→CGG)  cyoC ← cytochrome o ubiquinol oxidase subunit III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913448,4640AC 6.2% 311.2 / 10.6 88L66R (CTG→CGG) cyoCcytochrome o ubiquinol oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base A (34/48);  new base C (2/4);  total (36/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

TCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTG  >  NC_000913/448318‑448606
                                                                                                                                                  |                                                                                                                                              
tCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAATGTTTCAACCAGAACGAACGGCAGt                                                                                                                                              >  2:603360/1‑149 (MQ=255)
tCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGt                                                                                                                                              <  1:377031/149‑1 (MQ=255)
 cccGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGtt                                                                                                                                             <  2:409109/149‑1 (MQ=255)
  ccGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTc                                                                                                                                            >  1:408166/1‑149 (MQ=255)
   cGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCg                                                                                                                                           >  4:209950/1‑149 (MQ=255)
    gATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGa                                                                                                                                          <  4:512613/149‑1 (MQ=255)
     aTAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACTGCCGTTCgtt                                                                                                                                         <  1:294064/149‑3 (MQ=255)
     aTAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaa                                                                                                                                         >  2:589971/1‑149 (MQ=255)
      tAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACTGCCGTTCgttc                                                                                                                                        <  2:512576/149‑4 (MQ=255)
        aaTCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAAt                                                                                                                                      <  3:46907/149‑1 (MQ=255)
           ccGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTc                                                                                                                                   <  4:500229/149‑1 (MQ=255)
           ccGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTc                                                                                                                                   <  2:513334/149‑1 (MQ=255)
           ccGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTc                                                                                                                                   <  1:348742/149‑1 (MQ=255)
           ccGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTc                                                                                                                                   <  4:324158/149‑1 (MQ=255)
              gCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCtt                                                                                                                                >  2:446985/1‑149 (MQ=255)
                aCCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTAc                                                                                                                              <  1:446985/149‑1 (MQ=255)
                   aaaCAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaa                                                                                                                                         <  3:395269/135‑1 (MQ=255)
                   aaaCAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaa                                                                                                                                         >  4:395269/1‑135 (MQ=255)
                   aaaCAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTg                                                                                                                           >  2:550150/1‑149 (MQ=255)
                    aacaacCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGt                                                                                                                          <  4:220003/149‑1 (MQ=255)
                     acaacCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTc                                                                                                                         <  4:517843/149‑1 (MQ=255)
                       aacCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCgg                                                                                                                       <  2:601691/149‑1 (MQ=255)
                          cAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGgcc                                                                                                                    >  2:675630/1‑149 (MQ=255)
                            ggTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGgccgc                                                                                                                  <  1:205683/149‑1 (MQ=255)
                                aaCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGc                                                                                                              >  2:188129/1‑149 (MQ=255)
                                 aCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACTGCCGTTCgtt                                                                                                                                         >  1:485368/1‑119 (MQ=255)
                                 aCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACTGCCGTTCgtt                                                                                                                                         <  4:486305/121‑3 (MQ=255)
                                 aCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACTGCCGTTCgtt                                                                                                                                         >  3:486305/1‑119 (MQ=255)
                                 aCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACTGCCGTTCgtt                                                                                                                                         <  2:485368/121‑3 (MQ=255)
                                   gccagccaGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGt                                                                                                           <  3:100704/149‑1 (MQ=255)
                                          agGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTc                                                                                                    <  1:603360/149‑1 (MQ=255)
                                          agGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTc                                                                                                    <  1:550150/149‑1 (MQ=255)
                                           ggAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCa                                                                                                   <  2:442546/149‑1 (MQ=255)
                                            gagaTAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCAc                                                                                                  <  2:234283/149‑1 (MQ=255)
                                            gagaTAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCAc                                                                                                  <  1:675630/149‑1 (MQ=255)
                                                         ttttGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCg                                                                                                             <  3:215912/125‑1 (MQ=255)
                                                         ttttGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCg                                                                                                             >  4:215912/1‑125 (MQ=255)
                                                         ttttGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAg                                                                                     >  2:345498/1‑149 (MQ=255)
                                                          tttGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCg                                                                                                             <  2:655161/124‑1 (MQ=255)
                                                          tttGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCg                                                                                                             >  1:655161/1‑124 (MQ=255)
                                                          tttGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAgg                                                                                    <  3:383812/149‑1 (MQ=255)
                                                            tgttgtTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTa                                                                                  <  2:408166/149‑1 (MQ=255)
                                                              ttgttTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGc                                                                                >  1:370253/1‑149 (MQ=255)
                                                                 tttttGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCaaa                                                                             >  2:288453/1‑149 (MQ=255)
                                                                    ttGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACaa                                                                          <  1:167146/149‑1 (MQ=255)
                                                                    ttGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACaa                                                                          <  1:589971/149‑1 (MQ=255)
                                                                        aCATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATa                                                                      <  1:345498/149‑1 (MQ=255)
                                                                        aCATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGTGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATag                                                                     >  2:84799/1‑149 (MQ=255)
                                                                           tGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATagag                                                                   >  2:532134/1‑149 (MQ=255)
                                                                               gATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAAc                                                                <  4:680407/148‑1 (MQ=255)
                                                                               gATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAAc                                                                >  3:680407/1‑148 (MQ=255)
                                                                                aTAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAg                                                              >  3:175850/1‑149 (MQ=255)
                                                                                  aGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGaa                                                                                             >  2:676973/1‑116 (MQ=255)
                                                                                  aGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGaa                                                                                             <  1:676973/116‑1 (MQ=255)
                                                                                      gccATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCa                                                        <  1:532134/149‑1 (MQ=255)
                                                                                        cATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGt                                                      <  4:175850/149‑1 (MQ=255)
                                                                                        cATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGt                                                      <  3:519157/149‑1 (MQ=255)
                                                                                         aTGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTc                                                     <  1:140439/149‑1 (MQ=255)
                                                                                              gTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCa                                                                                                   >  2:368705/1‑98 (MQ=255)
                                                                                              gTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCa                                                                                                   <  1:368705/98‑1 (MQ=255)
                                                                                                aGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACGGAATGCAGTCGCtcatc                                              <  3:63687/149‑1 (MQ=255)
                                                                                                aGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTAGCtcatc                                              <  2:370253/149‑1 (MQ=255)
                                                                                                      tGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAg                                        >  2:57330/1‑149 (MQ=255)
                                                                                                      tGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAg                                        >  1:444563/1‑149 (MQ=255)
                                                                                                        gAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGAt                                      <  1:510257/149‑1 (MQ=255)
                                                                                                        gAGCTGAAAAAAAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGAt                                      <  1:188129/149‑1 (MQ=255)
                                                                                                          gCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATcc                                    >  4:398626/1‑149 (MQ=255)
                                                                                                              aacaacAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCa                                                                                        >  1:187134/1‑93 (MQ=255)
                                                                                                              aacaacAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCa                                                                                        <  2:187134/93‑1 (MQ=255)
                                                                                                                 aacaGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCg                                                    <  3:326317/126‑1 (MQ=255)
                                                                                                                 aacaGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCg                                                    >  4:326317/1‑126 (MQ=255)
                                                                                                                 aacaGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATcc                             >  2:68017/1‑149 (MQ=255)
                                                                                                                   caGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGa                           <  1:52905/149‑1 (MQ=255)
                                                                                                                        aGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCAt                                                                                       >  3:644095/1‑84 (MQ=255)
                                                                                                                        aGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCAt                                                                                       <  4:644095/84‑1 (MQ=255)
                                                                                                                        aGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGt                                                      <  4:332986/117‑1 (MQ=255)
                                                                                                                        aGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGt                                                      >  3:332986/1‑117 (MQ=255)
                                                                                                                        aGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGAtt                      >  3:681138/1‑149 (MQ=255)
                                                                                                                           aaGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTg                   >  4:360718/1‑149 (MQ=255)
                                                                                                                                tCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTcc              >  2:479426/1‑149 (MQ=255)
                                                                                                                                 cAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCg             <  1:479426/149‑1 (MQ=255)
                                                                                                                                    ccAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAg                                            >  1:388153/1‑115 (MQ=255)
                                                                                                                                    ccAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAg                                            <  2:388153/115‑1 (MQ=255)
                                                                                                                                           gaacgGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGc            >  3:6675/1‑140 (MQ=255)
                                                                                                                                           gaacgGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGc            <  4:6675/140‑1 (MQ=255)
                                                                                                                                           gaacgGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGTGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCgtg  <  1:84799/149‑1 (MQ=255)
                                                                                                                                               ggCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGc            <  4:243528/136‑1 (MQ=255)
                                                                                                                                               ggCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGc            >  3:243528/1‑136 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTG  >  NC_000913/448318‑448606

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: