Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 45,978 C→T L58L (CTG→TTG)  yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091345,9780CT93.2% 159.9 / 2.5 60L58L (CTG→TTG) yaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (25/30);  total (28/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 = 3.03e-01

AGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAG  >  NC_000913/45734‑46223
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ggTTCCGTTACGGCGGATGTTGGTTTGATACGTAACGCCGCACTGCCTCTCATTGCAAAAAACAGGAATACCCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGgcg                                                                                                                                                                                                                                                 <  1:396493/250‑1 (MQ=255)
   tCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGcgcg                                                                                                                                                                                                                                               <  1:73772/250‑1 (MQ=255)
      gTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGa                                                                                                                                                                                                                                            >  1:29988/1‑250 (MQ=255)
               gATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGCTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCg                                                                                                                                                                                                                                  >  1:256221/1‑251 (MQ=255)
                 ggTTGGTTTGATACGTATCCCTGCACTGACTCTCATTGCAAAAACCAGAAATAACCATGCAACCTTCCAGAACCTTTGACTATCTTAATTTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGTTAATGATTGTCGTGGCGCGGGAGCAACTGACGCCTGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCggg                                                                                                                                                                                                                                <  1:346234/250‑1 (MQ=255)
                  gTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGTTGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGc                                                                                                                                                                                                                               >  1:310118/1‑251 (MQ=255)
                    tGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCt                                                                                                                                                                                                                              <  1:107192/250‑1 (MQ=255)
                     ggTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCACCCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGt                                                                                                                                                                                                                            <  1:338774/251‑1 (MQ=255)
                                tAACGCCGCACTGACTCTCACTGCAAAAAACAGGAATAAACATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGcac                                                                                                                                                                                                                  >  1:189842/1‑250 (MQ=255)
                                   cgccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCg                                                                                                                                                                                                              >  1:460990/1‑251 (MQ=255)
                                     ccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCACCCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTCTGGACTGGTAATGATTGGCGTGGAGCGGGAGCTACTGACACCTGCGCTAAAACTGGACGCTGACGGGCTTGGATTGGTGGGAGCGGGGACGCCCGCCGGGATGTTCGGTggacacttggt                                                                                                                                                                                                            >  2:311097/1‑242 (MQ=255)
                                      cgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTg                                                                                                                                                                                                           >  2:135750/1‑251 (MQ=255)
                                        cACTGACTCTAATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATACACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAATACTGGACGCTGCCTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCAGTTCGTTGGCACAACGCTGtt                                                                                                                                                                                                         >  1:422398/1‑251 (MQ=255)
                                             aCTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTggct                                                                                                                                                                                                    <  1:436255/251‑3 (MQ=255)
                                                  cATTGCAAAAAACAGGAATAACCATGCACCCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCtata                                                                                                                                                                                                 <  1:265722/249‑1 (MQ=255)
                                                        aaaaaaCAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGAt                                                                                                                                                                                         >  1:214971/1‑251 (MQ=255)
                                                         aaaaaCAGGAAAACCAGGGCATCCGTGCAGATACTTGGTTGATCTCAAATTCTCCTTGATTCACCGCCGGATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGGTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACGGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCTCATCGCTGTTTGGCTATATTTCCGATa                                                                                                                                                                                        <  2:188395/251‑1 (MQ=255)
                                                            aaCAGGAATACCCATGCACCCGGCCAATAACTTTGGCGATCTCAAATTCCCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAg                                                                                                                                                                                     <  2:335891/251‑1 (MQ=255)
                                                                    aaaCCATGCAACCTTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCGCCGCCGTATTTTGCTGTGGGGAAGCGGTGGGCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGTGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACgg                                                                                                                                                                             <  1:479600/250‑1 (MQ=255)
                                                                     aaCCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGtcggtcggg                                                                                                                                                                            >  2:443473/1‑250 (MQ=255)
                                                                        cATGCAACCGTCCAGAAACTTTGACGATCTCACATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGTATGGTTATGTACTGGTACTGGTTGGCGTGGCTCTGGAGCAACTGACTCCGGCTCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGTACGCACACCGGGCTGTTCTTTGACCCAGCGCTGTTTGTCTATATTTCCGGAAAAGACGGAccgccac                                                                                                                                                                         >  2:1836/1‑245 (MQ=255)
                                                                                ccGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGATCGTTGGCACATCGCTGTTTGGCTATATTTCCGCTAAAGTCGGACGGCGACAAATGtt                                                                                                                                                                  >  2:124481/1‑251 (MQ=255)
                                                                                    ccAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGGTCCTc                                                                                                                                                              >  2:286203/1‑251 (MQ=255)
                                                                                        aaaCTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAAc                                                                                                                                                                                                                                          <  2:155128/171‑1 (MQ=255)
                                                                                        aaaCTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAAc                                                                                                                                                                                                                                          >  1:155128/1‑171 (MQ=255)
                                                                                        aaaCTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCCCAACGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTg                                                                                                                                                          >  2:248896/1‑251 (MQ=255)
                                                                                          aCTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGat                                                                                                                                                        <  1:3822/251‑1 (MQ=255)
                                                                                                gcGATCACAAATTCTCCTTTAGACCCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGATATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGGTCGTTGGCACATCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATatcatc                                                                                                                                                  <  2:108180/250‑1 (MQ=255)
                                                                                                aCGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATatcatc                                                                                                                                                  >  1:477456/1‑251 (MQ=255)
                                                                                                   aTCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGcc                                                                                                                                               >  1:218949/1‑251 (MQ=255)
                                                                                                               cctcTATTCGCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGata                                                                                                                                   <  1:370361/251‑1 (MQ=255)
                                                                                                                 tcgaTTCACCGCCGCCTTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTAATGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCGGGCGCTAAAACTGGACGCTGAGTGGCTTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCGCATCGCTGTTTGGCTGTATTTTCGATAAAGACGGAGGTCGCAGAATGTTCCTCATTGATATCATCGCCATCTGCGTGgtat                                                                                                                                  >  2:259050/4‑250 (MQ=255)
                                                                                                                  ctATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGatat                                                                                                                                  <  1:135750/249‑1 (MQ=255)
                                                                                                                      tCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTgg                                                                                                                            >  2:225618/1‑251 (MQ=255)
                                                                                                                       cACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCAGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAACCTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGTTATCATCGCCATCGGCGTGATATCGGTGGc                                                                                                                           >  2:467328/1‑251 (MQ=255)
                                                                                                                           gccgcATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGcgacg                                                                                                                       >  1:443181/1‑251 (MQ=255)
                                                                                                                                     gCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCTCATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGAtgtttgtt                                                                                                              <  1:38384/250‑1 (MQ=255)
                                                                                                                                             gAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATccccc                                                                                                     >  1:55760/1‑251 (MQ=255)
                                                                                                                                                      gtCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTg                                                                                            <  1:481860/251‑1 (MQ=255)
                                                                                                                                                         cGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGGGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGGCGCTGTCTGGATTGGCTTGTTGGTCGCGGGAACGCTCGTCGGGTTGTTCGTTGGCACATCGCTGTGTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCCTCGTCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTg                                                                                         <  1:198042/251‑1 (MQ=255)
                                                                                                                                                              cTGGATGGTTATGTACTTGTAATGATTGTCTTGGCGCTGGAGCAACTGACGTCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTTTTGTGCGCGGGAACGCTCGCCGGGCTTGTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGTCGATGTTTGTTTCATCCCCCGTCg                                                                                                 <  2:94003/238‑1 (MQ=255)
                                                                                                                                                              cTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTAGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCg                                                                                                 >  1:94003/1‑238 (MQ=255)
                                                                                                                                                                ggATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGa                                                                                     <  1:443473/248‑1 (MQ=255)
                                                                                                                                                                ggATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACAACGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATg                                                                                   <  1:248896/250‑1 (MQ=255)
                                                                                                                                                                ggATGGTTATGTACTGGTAATGATTGGCGTGGCCCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATg                                                                                   <  1:232692/250‑1 (MQ=255)
                                                                                                                                                                        aTGTTCTGGTAATGATTGGGGGGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTTGAGTGGCTTGTTGGGCGCGGGACCGCTCGCCGGGCTGTTCGTTGGCACATCGCTTTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACtt                                                                          <  2:474904/251‑1 (MQ=255)
                                                                                                                                                                            aCTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGTATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCg                                                                      <  2:421628/251‑1 (MQ=255)
                                                                                                                                                                            aCTGGTAATGAGTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATc                                                                       <  2:113755/250‑1 (MQ=255)
                                                                                                                                                                                      aTTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCAt                                                            <  1:124481/251‑1 (MQ=255)
                                                                                                                                                                                           cgtggcgCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGGGATGCGGGTACTTATCGGCATTGTCATCGGTg                                                       >  1:357365/1‑251 (MQ=255)
                                                                                                                                                                                               gcgcTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGc                                                      >  1:267522/1‑248 (MQ=255)
                                                                                                                                                                                                       gCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCa                                           <  1:225618/251‑1 (MQ=255)
                                                                                                                                                                                                         aaCTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATc                                         <  1:377164/251‑1 (MQ=255)
                                                                                                                                                                                                                                gACGCTGACTGGATTTGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGTACGGCGCAAAATGTTCCTCATTTATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTc                                                                                                  <  2:28417/171‑1 (MQ=255)
                                                                                                                                                                                                                                gACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTc                                                                                                  >  1:28417/1‑171 (MQ=255)
                                                                                                                                                                                                                                gACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAAt                              <  1:43444/239‑1 (MQ=255)
                                                                                                                                                                                                                                gACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGGACTTATCGGCATTGTCATCGGGGCAGATTATCCCATCGCCACCTCAAt                              >  2:43444/1‑239 (MQ=255)
                                                                                                                                                                                                                                      gACTGGATTGGCTTGTTGGGCGCGGGAACGCTAGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCAGATAAAGCCGGACGGCGCCAAATGTTCCTCATTGCTATCATCGCCATCGGCGTGCTATCGGTGGCGACGATGTTTGTTTCATCACCCGTCGAACTGTTGGCGCTGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCCATGATCACCGAGTTCTcc              >  1:66930/1‑249 (MQ=255)
                                                                                                                                                                                                                                           gATTGGCTTGTTGGGCGCGGGAACGCTCGCCGTGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTAccc       <  1:16195/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTGTGGCTATATTTCCGATAAAGTCGGGCGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCGCAATGATCACCGAGTTCTCCAGTACCCGTCAg  <  2:443181/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
AGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAG  >  NC_000913/45734‑46223

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