Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,461 G→T 7.1% P67Q (CCG→CAG)  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,4610GT 7.1% 298.8 / 10.5 85P67Q (CCG→CAG) cyoCcytochrome o ubiquinol oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base G (32/47);  new base T (2/4);  total (34/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

CCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTG  >  NC_000913/448315‑448606
                                                                                                                                                  |                                                                                                                                                 
ccATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGc                                                                                                                                                 <  3:233717/149‑1 (MQ=255)
ccATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGc                                                                                                                                                 <  1:178679/149‑1 (MQ=255)
   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:324158/149‑1 (MQ=255)
              ccGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTc                                                                                                                                   <  1:348742/149‑1 (MQ=255)
              ccGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTc                                                                                                                                   <  2:513334/149‑1 (MQ=255)
              ccGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTc                                                                                                                                   <  4:500229/149‑1 (MQ=255)
                 gCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCtt                                                                                                                                >  2:446985/1‑149 (MQ=255)
                   aCCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTAc                                                                                                                              <  1:446985/149‑1 (MQ=255)
                      aaaCAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaa                                                                                                                                         >  4:395269/1‑135 (MQ=255)
                      aaaCAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaa                                                                                                                                         <  3:395269/135‑1 (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                                                                                                                                         <  2:485368/121‑3 (MQ=255)
                                    aCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACTGCCGTTCgtt                                                                                                                                         >  3:486305/1‑119 (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                                                                                                  <  1:675630/149‑1 (MQ=255)
                                               gagaTAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCAc                                                                                                  <  2:234283/149‑1 (MQ=255)
                                                            ttttGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCg                                                                                                             >  4:215912/1‑125 (MQ=255)
                                                            ttttGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCg                                                                                                             <  3:215912/125‑1 (MQ=255)
                                                            ttttGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAg                                                                                     >  2:345498/1‑149 (MQ=255)
                                                             tttGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCg                                                                                                             >  1:655161/1‑124 (MQ=255)
                                                             tttGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCg                                                                                                             <  2:655161/124‑1 (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                                                                                        <  2:187134/93‑1 (MQ=255)
                                                                                                                 aacaacAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCa                                                                                        >  1:187134/1‑93 (MQ=255)
                                                                                                                    aacaGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCg                                                    >  4:326317/1‑126 (MQ=255)
                                                                                                                    aacaGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCg                                                    <  3:326317/126‑1 (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            >  3:243528/1‑136 (MQ=255)
                                                                                                                                                  ggCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGc            <  4:243528/136‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTG  >  NC_000913/448315‑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: