Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,982,585 C→A 6.8% R319L (CGT→CTT)  araH ← L‑arabinose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,982,5850CA 6.8% 193.1 / 3.5 59R319L (CGT→CTT) araHL‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base A (2/2);  total (27/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 = 8.78e-01

GACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTT  >  NC_000913/1982440‑1982733
                                                                                                                                                  |                                                                                                                                                    
gagaTGTAAGTATTGACAGTGTGACGGGGGGCGGCATTTGCAACGAAGAACAGCCAACGGTTATGGCGGTTGGCTATGGTGGGAAAAAAGGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGt                                                                                                                                                   <  6:104887/146‑1 (MQ=255)
  cATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGTAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAATTTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCg                                                                                                                                                 <  8:275742/149‑1 (MQ=255)
     gTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCGAAATTGTTGCAG‑CAAAAAGCATCAGCCAGTGCGTTTCGCTTTTTGCTTGTAACGGGCg                                                                                                                                                 >  4:219899/1‑146 (MQ=255)
      tAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCg                                                                                                                                                 >  7:189156/1‑145 (MQ=255)
          tataGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAc                                                                                                                                         <  1:248137/149‑1 (MQ=255)
            taGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAcc                                                                                                                                        <  7:38908/148‑1 (MQ=255)
            taGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCg                                                                                                                                       >  7:278764/1‑149 (MQ=255)
                 aGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGc                                                                                                                                   >  3:306583/1‑148 (MQ=255)
                  gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGc                                                                                                                                   >  8:66822/1‑147 (MQ=255)
                       aCGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCcagcag                                                                                                                             >  7:220518/1‑148 (MQ=255)
                         gggggTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGAt                                                                                                                          <  7:113745/149‑1 (MQ=255)
                         gggggTCGCCATTTGAAAATAAAAAAAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAGAAACGCGATATTGGTGCAG‑CAAACAGCATCGGCCGGTGCGTTTGGCTGTTTGCTTGTATCGGTCGCAGATCAGCGCTGCCAGCAGGAt                                                                                                                          >  6:182321/1‑149 (MQ=255)
                           gggTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGAtt                                                                                                                         >  4:186016/1‑148 (MQ=255)
                                gCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTa                                                                                                                        >  7:48196/1‑144 (MQ=255)
                                gCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTa                                                                                                                        <  8:48196/144‑1 (MQ=255)
                                       gCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCcagca                                                                                                                              >  6:11252/1‑131 (MQ=255)
                                       gCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCcagca                                                                                                                              <  5:11252/131‑1 (MQ=255)
                                        cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGGGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGAGTCAGCCGCGAACCAc                                                                                                           >  6:77462/1‑149 (MQ=255)
                                         aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCcagc                                                                                                                               <  7:62177/128‑1 (MQ=255)
                                         aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCcagc                                                                                                                               >  8:62177/1‑128 (MQ=255)
                                         aaCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACg                                                                                                          >  2:41775/1‑149 (MQ=255)
                                                     ccAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGTCGCGAACCACGTACTGCGCGaaa                                                                                              <  1:284796/149‑1 (MQ=255)
                                                     ccAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaaa                                                                                              <  5:131290/149‑1 (MQ=255)
                                                      cAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAg                                                                                             <  3:186016/149‑1 (MQ=255)
                                                       aaCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAg                                                                                             >  2:304951/1‑148 (MQ=255)
                                                       aaCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAg                                                                                             >  3:250542/1‑148 (MQ=255)
                                                            ttATGGCGGGTGGCTATGGTGGGAAAAAACGCGAATTTGTTGCAG‑AAAAAAGCTTCAGACAGTGCGTTTCGCTTTTTCCTTGTAACGGACGAAGATAACCCCTACCACAAGGATTAAGCCGCGAACAACGAACTGCGCGAAAgaagaa                                                                                        >  1:239647/1‑148 (MQ=255)
                                                                 gCGGTTGGCTATTGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATAtt                                                                                   <  4:238955/149‑1 (MQ=255)
                                                                       ggCTATGGTTGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTTAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAg                                                                             <  6:277774/148‑1 (MQ=255)
                                                                          tATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGtt                                                                          >  6:116976/1‑148 (MQ=255)
                                                                          tATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAAGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGtt                                                                          >  2:27493/1‑148 (MQ=255)
                                                                          tATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAAGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGtt                                                                          >  4:287921/1‑148 (MQ=255)
                                                                            tggtggGAAAATACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACGACGTACTGCGGGAAAGGAGAAATATTAAGCAGGTTCa                                                                        <  2:239647/148‑1 (MQ=255)
                                                                              gtggGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATg                                                                      <  8:220518/148‑1 (MQ=255)
                                                                                 ggAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGt                                                                  >  7:8438/1‑149 (MQ=255)
                                                                                   aaaaaaCGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGttt                                                                >  5:51514/1‑149 (MQ=255)
                                                                                   aaaaaaCGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACATACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGttt                                                                <  4:306583/149‑1 (MQ=255)
                                                                                     aaaaCGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTc                                                              <  4:16675/149‑1 (MQ=255)
                                                                                     aaaaCGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTc                                                              <  6:200418/149‑1 (MQ=255)
                                                                                       aaCGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCTCCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCa                                                            <  1:6894/149‑1 (MQ=255)
                                                                                          gCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACgg                                                         >  1:140042/1‑149 (MQ=255)
                                                                                              aaTTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGc                                                      <  8:278764/148‑1 (MQ=255)
                                                                                                     gCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAAt                                               >  7:256193/1‑148 (MQ=255)
                                                                                                      cAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAAtt                                              >  6:281450/1‑148 (MQ=255)
                                                                                                       aG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAAtt                                              >  6:110300/1‑147 (MQ=255)
                                                                                                       aG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTaa                                            <  3:219899/149‑1 (MQ=255)
                                                                                                       aG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTaa                                            <  2:140042/149‑1 (MQ=255)
                                                                                                       aG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAAGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTaa                                            <  3:287921/149‑1 (MQ=255)
                                                                                                       aG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAAGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTaa                                            <  1:27493/149‑1 (MQ=255)
                                                                                                        g‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAg                                           <  6:51514/149‑1 (MQ=255)
                                                                                                            aaaaGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATa                                        <  5:182321/149‑1 (MQ=255)
                                                                                                               aGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATAccc                                     <  8:270760/149‑1 (MQ=255)
                                                                                                                      gACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGccacca                              <  1:41775/149‑1 (MQ=255)
                                                                                                                      gACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTGCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGccacca                              <  4:250542/149‑1 (MQ=255)
                                                                                                                         aGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAg                                                                             >  4:189971/1‑99 (MQ=255)
                                                                                                                         aGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAg                                                                             <  3:189971/99‑1 (MQ=255)
                                                                                                                           tGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGccac                                <  4:269634/142‑1 (MQ=255)
                                                                                                                           tGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGccac                                >  3:269634/1‑142 (MQ=255)
                                                                                                                                       tttGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGCCACCACATATGAGATTTTTc               <  7:66822/147‑1 (MQ=255)
                                                                                                                                                  cGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCTCCACCACATATGAGATTTTTCCGATGCCACCttt  <  1:304951/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAG‑AAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCCTAAAATTAAGATACCCGCCACCACATATGAGATTTTTCCGATGCCACCTTT  >  NC_000913/1982440‑1982733

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑