Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 645,711 C→T L290L (CTG→CTA citT ← citrate:succinate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913645,7110CT93.7% 189.2 / 3.1 64L290L (CTG→CTAcitTcitrate:succinate antiporter
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (29/30);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

GGTAATGTCTTTCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATCAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAGCAAGATA  >  NC_000913/645582‑645846
                                                                                                                                 |                                                                                                                                       
ggTAATGTCTTTCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAGTACAAGGCCAATTAGTGTCCAct                                                                                                                               >  2:673507/1‑140 (MQ=255)
ggTAATGTCTTTCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCAct                                                                                                                               >  2:175784/1‑140 (MQ=255)
  tAATGTCTTTCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCActct                                                                                                                             >  1:993381/1‑140 (MQ=255)
  tAATGTCTTTCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCActct                                                                                                                             >  1:653522/1‑140 (MQ=255)
  tAATGTCTTTCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCActct                                                                                                                             >  2:664245/1‑140 (MQ=255)
  tAATGTCTTTCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCActct                                                                                                                             >  1:576188/1‑140 (MQ=255)
  tAATGTCTTTCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCActct                                                                                                                             >  1:540227/1‑140 (MQ=255)
  tAATGTCTTTCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCActcc                                                                                                                             >  2:998702/1‑139 (MQ=255)
          ttCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCgtgt                                                                                                                     <  2:993381/140‑1 (MQ=255)
           tCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTc                                                                                                                    >  1:912510/1‑140 (MQ=255)
                ggCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAgcgc                                                                                                               >  1:946072/1‑140 (MQ=255)
                ggCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAgcgc                                                                                                               <  1:536636/140‑1 (MQ=255)
                     aaCGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCa                                                                                                          >  2:585813/1‑140 (MQ=255)
                     aaCGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATCAGTGTCCACTCTCTGCGTGTCAGCGCACCCa                                                                                                          >  1:750906/1‑140 (MQ=255)
                       cGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATg                                                                                                        <  2:847299/140‑1 (MQ=255)
                       cGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATg                                                                                                        <  2:875395/140‑1 (MQ=255)
                           cAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATAAGTGTCCTCTCCCTACGTGTCGGCGCACCCATGGttt                                                                                                    >  1:451027/1‑140 (MQ=255)
                            aGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATCAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGtttt                                                                                                   <  2:750906/140‑1 (MQ=255)
                                    cATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGTGCACCCATGGTTTTCAGTtcat                                                                                           >  2:246060/1‑140 (MQ=255)
                                     aTTAGCGAAACTGCCCGCAGACCAACCGCAGTAGCATCAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTtcatc                                                                                          <  1:497278/140‑1 (MQ=255)
                                      ttAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTtcatca                                                                                         <  2:912510/140‑1 (MQ=255)
                                         gCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCAccc                                                                                      >  2:172010/1‑140 (MQ=255)
                                            aaaCTGCCAGCAGACCAACCGCTGTAGCATTTATGCCTTGACTGCCAATTACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGcc                                                                                   <  2:375717/140‑1 (MQ=255)
                                                    agcagACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGc                                                                           <  2:832557/140‑1 (MQ=255)
                                                      cagACCAACCGCAGTCGCATTACTGTCTTTACTGCCACATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCa                                                                         <  1:310670/140‑1 (MQ=255)
                                                      cagACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGTGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCa                                                                         <  1:246060/140‑1 (MQ=255)
                                                      cagACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCa                                                                         <  1:175784/140‑1 (MQ=255)
                                                      cagACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCa                                                                         <  1:157739/140‑1 (MQ=255)
                                                            aaCCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCActtctt                                                                   >  2:200658/1‑140 (MQ=255)
                                                             aCCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGtt                                                                             <  2:955067/129‑1 (MQ=255)
                                                             aCCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGtt                                                                             >  1:955067/1‑129 (MQ=255)
                                                              ccGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCa                                                                 <  2:576188/140‑1 (MQ=255)
                                                              ccGCAGTAGCATTAACGACTTCACTGCCAAATACCCAGAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCa                                                                 <  1:357330/140‑1 (MQ=255)
                                                                gCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGGACCCATGGGTTTCAGTACATCACCCTCCCAGGTTGCCACTTCTTCACt                                                               >  2:357330/1‑140 (MQ=255)
                                                                gCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCATCCGCCCAGGTTGCCACTTCTTCACt                                                               <  1:227655/140‑1 (MQ=255)
                                                                 cAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACtg                                                              <  1:998702/140‑1 (MQ=255)
                                                                      gCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACtgtgtgt                                                         <  1:528029/140‑1 (MQ=255)
                                                                       cATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACtgtgtgtg                                                        >  2:799711/1‑140 (MQ=255)
                                                                         ttAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGAt                                                      >  1:92565/1‑140 (MQ=255)
                                                                           aaTGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGAttt                                                    <  1:381479/140‑1 (MQ=255)
                                                                              gACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCg                                                 >  2:1013555/1‑140 (MQ=255)
                                                                              gACTTCACTGCCAAAAACCCACTAACCTAGGCTGAGCTATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCg                                                 <  1:356725/140‑1 (MQ=255)
                                                                                 ttCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGtt                                              >  1:970761/1‑140 (MQ=255)
                                                                                 ttCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGACACATCCTCACTGTGTGTGATTTCCGGtt                                              >  1:418568/1‑140 (MQ=255)
                                                                                     cTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTa                                          <  2:2722/140‑1 (MQ=255)
                                                                                       gCCAAATACCCACAAGCTTAAGTTGAGCAATGCAAGGCCCATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACa                                        <  2:308738/140‑1 (MQ=255)
                                                                                       gCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACa                                        <  1:1013555/140‑1 (MQ=255)
                                                                                          aaaTACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCTACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCa                                     <  2:621252/140‑1 (MQ=255)
                                                                                          aaaTACCCACAAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCa                                     <  1:799711/140‑1 (MQ=255)
                                                                                                  acaAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGaa                             >  1:621009/1‑140 (MQ=255)
                                                                                                  acaAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGaa                             >  1:905829/1‑140 (MQ=255)
                                                                                                  acaAACCTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTATGATTTCCGGTTTCTACAGCACGTAGGaa                             >  1:393734/1‑140 (MQ=255)
                                                                                                       ccTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCa                        >  1:788052/1‑140 (MQ=255)
                                                                                                       ccTAAGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCa                        <  1:587404/140‑1 (MQ=255)
                                                                                                          aaGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACgg                     >  1:946937/1‑140 (MQ=255)
                                                                                                          aaGCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACgg                     >  1:665494/1‑140 (MQ=255)
                                                                                                            gCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGgcg                   <  1:587259/140‑1 (MQ=255)
                                                                                                            gCTGAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGgcg                   <  2:653522/140‑1 (MQ=255)
                                                                                                               gAGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCaa                <  1:540000/140‑1 (MQ=255)
                                                                                                                aGCAATACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAAt               >  2:885066/1‑140 (MQ=255)
                                                                                                                 gCAATACAAGGCCAATCAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATg              >  2:770617/1‑140 (MQ=255)
                                                                                                                    aTACAAGGCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATa           <  2:621009/140‑1 (MQ=255)
                                                                                                                      aCAAGGCCAATCAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAg         <  1:770617/140‑1 (MQ=255)
                                                                                                                           gCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAGCAAGa    <  2:699905/140‑1 (MQ=255)
                                                                                                                           gCCAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCGTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGTAAGCCACGGCGCAATGATAAGCAAGa    <  2:905829/140‑1 (MQ=255)
                                                                                                                             cAATTAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAGCAAGATa  >  1:1009170/1‑140 (MQ=255)
                                                                                                                                 |                                                                                                                                       
GGTAATGTCTTTCCAAGGCACAACGTGCAGGGCCAGCATTAGCGAAACTGCCAGCAGACCAACCGCAGTAGCATTAATGACTTCACTGCCAAATACCCACAAACCTAAGCTGAGCAATACAAGGCCAATCAGTGTCCACTCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAGCAAGATA  >  NC_000913/645582‑645846

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