Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,577,852 C→T L497L (CTG→CTA yddA ← ABC transporter family protein YddA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,577,8520CT75.5% 84.4 / 27.0 50L497L (CTG→CTAyddAABC transporter family protein YddA
Reads supporting (aligned to +/- strand):  ref base C (5/7);  new base T (18/19);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

CCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACACATACCAACTGTCAGCAGGAGAAGAAATATCACCTTTAAACCACGGCCAGCAGTGGGATAATGTTTTAAGCA  >  NC_000913/1577571‑1578117
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                         
ccACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCg                                                                                                                                                                                                                                                                     <  1:425427/288‑1 (MQ=255)
ccACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATAGCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCg                                                                                                                                                                                                                                                                     <  1:234706/288‑1 (MQ=255)
   caTCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTGTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGGGATa                                                                                                                                                                                                                                                                  >  1:224982/1‑288 (MQ=255)
    aTCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCCCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGGTAATGTCACAAATATCATAGGCCAGGTTCCAGACAACGGGTTGTTGTGTAACCATATTAACGCGGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCGCCAGATTAGCGGTAAGTTCGTACAGAAAAATGCCTTACGGACGTCGTAAAATTAAATGTGGGAGGGNGNTAc                                                                                                                                                                                                                                                                 >  2:175785/1‑288 (MQ=255)
              caTAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAGATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGttttt                                                                                                                                                                                                                                                        >  1:418584/1‑287 (MQ=255)
                                aTAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAata                                                                                                                                                                                                                                       <  1:175785/286‑1 (MQ=255)
                                         cGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGCCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCg                                                                                                                                                                                                                            >  1:224269/1‑288 (MQ=255)
                                               atttatttATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTTCTAGGGCGATACGTTGTTTTTc                                                                                                                                                                                                                                                       >  1:84803/1‑255 (MQ=255)
                                                     ttatATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGa                                                                                                                                                                                                                                              <  1:328901/258‑1 (MQ=255)
                                                               gTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCa                                                                                                                                                                                                      <  1:363669/288‑1 (MQ=255)
                                                                          tttttttTGTTTTTAAAAACCCCGTTAAATGTCACAAAGATCAGCGGCCAGTTTCCAGAAACCGGGTTGATGGGTAACCATAATGACGCCGCTTGTGGGTAGTTTGTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCTGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTccc                                                                                                                                                                                           <  2:224269/287‑1 (MQ=255)
                                                                          cTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATATGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTccc                                                                                                                                                                                           >  1:371790/1‑288 (MQ=255)
                                                                                     tttnaaaccncaagcTAGTGTCACCGATAGAATCGGCCAGTTTCCTGCAACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATAcgcg                                                                                                                                                                                                        <  2:418584/251‑1 (MQ=255)
                                                                                                                                tCCGGACACCTGGCAGATGTTTAACCATAATGACGCCGCTTGTGGGTAGTGTTTCACGAACTAAAGGCAGTAAGCGGATAGCCTATTGGTCCTCGAGATGAGAGGTAGTTTCGTCAAGAAATTTCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCaa                                                                                                                                                                                                 <  2:272200/229‑1 (MQ=255)
                                                                                                                                         cGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGCGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAg                                                                                                                                                                                        >  1:272200/1‑229 (MQ=255)
                                                                                                                                                       aCCATAATGACGCCGCTTGTGGGTAGTTTTTCAAGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGAt                                                                                                                                                                                                                                                                   >  1:449716/1‑140 (MQ=255)
                                                                                                                                                           tAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGGGAATTGCTTTGCAAATAATCTCTTTCa                                                                                                           >  1:189458/1‑288 (MQ=255)
                                                                                                                                                             aTGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCg                                                                                                                                                                                                                        >  2:213410/1‑177 (MQ=255)
                                                                                                                                                             aTGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCg                                                                                                                                                                                                                        <  1:213410/177‑1 (MQ=255)
                                                                                                                                                             aTGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGt                                                                                                                                                                        >  1:695098/1‑225 (MQ=255)
                                                                                                                                                                  gccgcTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCg                                                                                                                                                                                                                                                                     >  1:565645/1‑127 (MQ=255)
                                                                                                                                                                  gccgcTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCg                                                                                                                                                                                                                                                                     <  2:565645/127‑1 (MQ=255)
                                                                                                                                                                   ccgcTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACgg                                                                                                                                            >  1:347306/1‑247 (MQ=255)
                                                                                                                                                                   ccgcTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGaa                                                                                                                                                                                                                                         <  1:165000/154‑1 (MQ=255)
                                                                                                                                                                   ccgcTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGaa                                                                                                                                                                                                                                         >  2:165000/1‑154 (MQ=255)
                                                                                                                                                                                gTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAAGACGCGCAGCCAATTTCCCGAGACCAACCTGATGCAGTACATCGCCCAACGATTTATCGTCTACGGGCAGGTGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACg                                                                                      <  1:171526/288‑1 (MQ=255)
                                                                                                                                                                                            cTAAACGCAGTAAGCGGATGGCCTCTTGTTCCTCAAGAGGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAGGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGGCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACg                                                                                      <  1:699188/276‑1 (MQ=255)
                                                                                                                                                                                                    aGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGc                                                                                                                                                                                                                             >  2:238024/1‑133 (MQ=255)
                                                                                                                                                                                                    aGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGc                                                                                                                                                                                                                             <  1:238024/133‑1 (MQ=255)
                                                                                                                                                                                                                  ctctTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACACATACCAACt                                                                <  1:4025/276‑1 (MQ=255)
                                                                                                                                                                                                                  ctctTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAAGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCAATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTAGCTCAACGATTTATCGTCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACACATACCAACTGTCAGCAGgaga                                                    <  1:211178/288‑1 (MQ=255)
                                                                                                                                                                                                                    ctTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAg                                                                                                          <  1:497425/232‑1 (MQ=255)
                                                                                                                                                                                                                      tGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAAcc                                                                                                                                                                                 <  1:400205/159‑1 (MQ=255)
                                                                                                                                                                                                                      tGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAAcc                                                                                                                                                                                 >  2:400205/1‑159 (MQ=255)
                                                                                                                                                                                                                          ccTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACACATACCAACTGTCAGCAGGAGAAGAAatat                                            <  1:81091/288‑1 (MQ=255)
                                                                                                                                                                                                                             cAAGTTGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGt                                                                                                                                                                        >  1:364362/1‑161 (MQ=255)
                                                                                                                                                                                                                                        gTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACTGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACACATACCAACTGTCAGCAGGAGAAGAAATATCACCTTTAAACCAc                              >  1:665154/1‑288 (MQ=255)
                                                                                                                                                                                                                                        gTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACACATACCAACTGTCAGCAGGAGAAGAAatat                                            <  1:161589/274‑1 (MQ=255)
                                                                                                                                                                                                                                           gTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTcgccg                                                                                                                                                                                                                                                   >  1:51014/1‑72 (MQ=255)
                                                                                                                                                                                                                                           gTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTcgccg                                                                                                                                                                                                                                                   <  2:51014/72‑1 (MQ=255)
                                                                                                                                                                                                                                                gTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTa                                                                                                                                               <  2:308803/167‑1 (MQ=255)
                                                                                                                                                                                                                                                gTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTa                                                                                                                                               >  1:308803/1‑167 (MQ=255)
                                                                                                                                                                                                                                                 tCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACCCATACAAACTGTCACAAGGAGAAAAAAAATAAACTTTtaaaaatgccatca                       >  2:161589/1‑275 (MQ=255)
                                                                                                                                                                                                                                                 tCAAAAAATATCCCTTTCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACACATACCAACTGTCAGCAGGAGAAGAAATATCACCTTTAAACCACGGCcagcag                     <  1:572135/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                ttCGGACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCCGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATAGTCTACGGGCAGGGGAAGTGCTGTACAAATAATCCCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACCCATACCAACTGTCAGCAGGAGAAGAAATATCACGTTTAAACCACGGCCAGCAGTGGGATAATGTTTTa      >  1:645145/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                  cGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACACATACCAACTGTCAGCAGGAGAAGAAATATCACCTTTAAACCACGGCCAGCAGTGGGATAATGTTTTAAg    <  1:364452/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gACGTCGTAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACACATACCAACTGTCAGCAGGAGAAGAAATAACACCTTTAAACCACGGCCAGCAGTGGGATAATGTTTTAAGCa  >  1:302661/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                          gtAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTcgccg                                                                                                                                                                                                                                                   <  1:27964/41‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gtAAAATTAATCGTGCTAGGGCGATACGTTGTTTTTcgccg                                                                                                                                                                                                                                                   >  2:27964/1‑41 (MQ=255)
                                                                                                                                                                                                                                                                                    tCGTGCTAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGTGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACACATACCAACTGTCAGCAGGAGAAGAAATATCACCTTTAAACCACGGcc                          <  1:271610/248‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                         
CCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGTAGTTTTTCACGCACTAAACGCAGTAAGCGGATAGCCTCTTGTTCCTCAAGATGAGAGGTAGTTTCGTCAAGAAATATCCATTTCGGACGTCGTAAAATTAATCGTGCCAGGGCGATACGTTGTTTTTCGCCGCTGGAAAGAATATCTCCCCAGCGATCGTGGTCATGAATACGCGCAGCCAATTTCCCAAGACCAACCTGATGCAGTACTTCGCTCAACGATTTATCGTCTACGGGCAGGGGAAGTGCTTTACAAATAATCTCTTTCAGTAAGCCGGTTTTGATTAACGGTGTTTGTGACACATACCAACTGTCAGCAGGAGAAGAAATATCACCTTTAAACCACGGCCAGCAGTGGGATAATGTTTTAAGCA  >  NC_000913/1577571‑1578117

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