Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,875 G→A 6.0% S589S (AGC→AGT cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913448,8750GA 6.0% 219.6 / 3.0 67S589S (AGC→AGTcyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base A (2/2);  total (34/33)
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.77e-01

ACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCAC  >  NC_000913/448729‑449022
                                                                                                                                                  |                                                                                                                                                   
aCGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCt                                                                                                                                                   >  3:270755/1‑149 (MQ=255)
     gtaATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGtttt                                                                                                                                              <  1:401745/149‑1 (MQ=255)
       aaTCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTc                                                                                                                                            <  1:558040/149‑1 (MQ=255)
        aTCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCg                                                                                                                                           >  4:573274/1‑149 (MQ=255)
          ccACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGc                                                                                                                                         >  4:163742/1‑149 (MQ=255)
            aCGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCat                                                                                                                                       <  4:181230/149‑1 (MQ=255)
            aCGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCat                                                                                                                                       >  1:544911/1‑149 (MQ=255)
              gTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCatat                                                                                                                                     <  1:244575/149‑1 (MQ=255)
                 cTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCACTGTTTTTCGGCATATGaa                                                                                                                                  >  1:502859/1‑149 (MQ=255)
                  tcgtcgAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAAt                                                                                                                                 >  3:460839/1‑149 (MQ=255)
                  tcgtcgAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAAt                                                                                                                                 <  4:44114/149‑1 (MQ=255)
                  tcgtcgAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAAt                                                                                                                                 <  1:215346/149‑1 (MQ=255)
                    gtcgAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAAttt                                                                                                                               <  3:289442/149‑1 (MQ=255)
                     tcgAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCACTGTTTTTCGGCATATGAATttc                                                                                                                              <  2:502859/149‑1 (MQ=255)
                       gAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATttctt                                                                                                                            >  2:286488/1‑149 (MQ=255)
                            ttttCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAg                                                                                                                       >  1:446241/1‑149 (MQ=255)
                                cACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGt                                                                                                                   >  3:47672/1‑149 (MQ=255)
                                   gATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAg                                                                                                                <  3:22549/149‑1 (MQ=255)
                                      ccAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCt                                                                                                             <  4:270755/149‑1 (MQ=255)
                                        aGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCttt                                                                                                           >  1:444194/1‑149 (MQ=255)
                                        aGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCttt                                                                                                           <  1:286488/149‑1 (MQ=255)
                                        aGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCttt                                                                                                           <  3:96026/149‑1 (MQ=255)
                                        aGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCttt                                                                                                           <  2:444194/149‑1 (MQ=255)
                                        aGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCttt                                                                                                           >  2:501415/1‑149 (MQ=255)
                                            gatgatCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTg                                                                                                                     >  3:295119/1‑135 (MQ=255)
                                            gatgatCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTg                                                                                                                     <  4:295119/135‑1 (MQ=255)
                                             atgatCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTtata                                                                                                      <  4:214688/149‑1 (MQ=255)
                                                       cTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCt                                                                                            <  3:156558/149‑1 (MQ=255)
                                                                caacaaTCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTc                                                                                   >  2:584701/1‑149 (MQ=255)
                                                                 aacaaTCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCa                                                                                  >  2:342925/1‑149 (MQ=255)
                                                                 aacaaTCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCa                                                                                  >  1:112073/1‑149 (MQ=255)
                                                                 aacaaTCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCACTGTTTTTCGGCatat                                                                                                                                     <  4:312092/98‑1 (MQ=255)
                                                                 aacaaTCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCACTGTTTTTCGGCatat                                                                                                                                     >  3:312092/1‑98 (MQ=255)
                                                                      tCGCCAGCCACTAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCGTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTcc                                                                             >  4:53517/1‑149 (MQ=255)
                                                                      tCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTcc                                                                             >  3:546044/1‑149 (MQ=255)
                                                                           agccaCCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGaa                                                                        <  4:197775/149‑1 (MQ=255)
                                                                                       tGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTc                                                            <  3:157826/149‑1 (MQ=255)
                                                                                         ccAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGt                                                          <  2:581445/149‑1 (MQ=255)
                                                                                         ccAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGt                                                          >  1:303006/1‑149 (MQ=255)
                                                                                         ccAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGt                                                          <  3:547688/149‑1 (MQ=255)
                                                                                           aGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGa                                                        >  3:66707/1‑149 (MQ=255)
                                                                                                aTGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTg                                                   >  4:401082/1‑149 (MQ=255)
                                                                                                 tGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGc                                                  <  4:201845/149‑1 (MQ=255)
                                                                                                  ggCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCg                                                 <  4:460839/149‑1 (MQ=255)
                                                                                                  ggCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgat                                                 <  4:7633/149‑3 (MQ=255)
                                                                                                  ggCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTgat                                                 <  1:140990/149‑3 (MQ=255)
                                                                                                    cGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGc                                               >  3:373716/1‑149 (MQ=255)
                                                                                                     gAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGc                                                                     >  1:517879/1‑126 (MQ=255)
                                                                                                     gAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGc                                                                     <  2:517879/126‑1 (MQ=255)
                                                                                                           gAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACgg                                        >  3:256380/1‑149 (MQ=255)
                                                                                                               aTGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCaaa                                    <  1:501415/149‑1 (MQ=255)
                                                                                                                  gtggAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGtt                                 >  4:367107/1‑149 (MQ=255)
                                                                                                                  gtggAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGtt                                 <  1:584701/149‑1 (MQ=255)
                                                                                                                         aaGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAAc                          >  2:50654/1‑149 (MQ=255)
                                                                                                                                 aTGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCgg                  >  3:256414/1‑149 (MQ=255)
                                                                                                                                  tGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCggg                 <  3:53517/149‑1 (MQ=255)
                                                                                                                                   gACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATTACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCgggg                <  2:176931/149‑1 (MQ=255)
                                                                                                                                      gATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAag             <  2:501288/149‑1 (MQ=255)
                                                                                                                                      gATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAag             <  4:551665/149‑1 (MQ=255)
                                                                                                                                        tACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAagag           <  4:256380/149‑1 (MQ=255)
                                                                                                                                         aCCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAgg          >  3:305990/1‑149 (MQ=255)
                                                                                                                                           cTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGtt        >  2:281491/1‑149 (MQ=255)
                                                                                                                                           cTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGtt        >  1:369518/1‑149 (MQ=255)
                                                                                                                                           cTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGtt        >  1:300300/1‑149 (MQ=255)
                                                                                                                                            tGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTg       >  1:275404/1‑149 (MQ=255)
                                                                                                                                               aCCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGccc    >  3:305317/1‑149 (MQ=255)
                                                                                                                                                 cGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCAc  >  2:15664/1‑149 (MQ=255)
                                                                                                                                                  gCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTga                                                  <  4:171652/100‑2 (MQ=255)
                                                                                                                                                  gCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCACGAACGTga                                                  >  3:171652/1‑99 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATACCTGCACCGCTGTTTTTCGGCATATGAATTTCTTCATAGTGGTCAGGCTTTTTATACGCTTCGCCTTTCTCTTTCATTTCCCAGAATGCATCACGTTCGTGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCAC  >  NC_000913/448729‑449022

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: