Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 471,899 C→T Q422* (CAA→TAA)  mdlB → ABC transporter family protein MdlB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913471,8990CT98.5% 205.2 / ‑4.1 67Q422* (CAA→TAA) mdlBABC transporter family protein MdlB
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (35/31);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.17e-01

GAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGCAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCGTGCTGGTCGAGACGCCGCAAATCCTGATCCTTGATGAGGCA  >  NC_000913/471659‑472144
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        
gAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGAGGCGCTGGTCGCGCATACCGGCAGTGGCAAAAGCACCCTCGCAAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCggg                                                                                                                                                                                                                                                >  2:365115/1‑250 (MQ=255)
     cGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAACATTAATCTCTCGGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCTCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGgtggtg                                                                                                                                                                                                                                            <  1:439938/250‑1 (MQ=255)
       aTAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTg                                                                                                                                                                                                                                         >  1:435109/1‑251 (MQ=255)
       aTAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCGTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGAGTCGCCTTGATGGGCGGCCATTAAGGTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGGTCCGGTGGTTCTg                                                                                                                                                                                                                                         <  2:24550/251‑1 (MQ=255)
          tcatgaaattggcttctagcgATGAAAAACTGGTTATAAAAACCATTAATATCTGTGGGCCTTCGCGCAATATTGGGGCGATGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCAAGTTTATAGATGGGCTATTAACCGCTAACGGAAGGGGAGATACGCCTTGTTGGTCGACAGTTAAGTACGCTACGTCATAGCGCGCTGCGCCCGGGCGTGGAAATGGAGCAGTAAGGTCCGGTAGTGcagacg                                                                                                                                                                                                                                      >  2:276504/19‑246 (MQ=255)
             tgtCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGAt                                                                                                                                                                                                                                   >  1:424918/1‑251 (MQ=255)
                            gcgATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGTTTCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAAc                                                                                                                                                                                                                    <  1:402884/251‑1 (MQ=255)
                            gcgATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAAc                                                                                                                                                                                                                    >  1:49681/1‑251 (MQ=255)
                                 gACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCACTTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTCACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCACGCCAACgggcc                                                                                                                                                                                                               >  2:291309/1‑247 (MQ=255)
                                    aaTCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCt                                                                                                                                                                                                            >  1:183796/1‑251 (MQ=255)
                                     atatGGTGCTAACGTACATTAATCTCTCTGTGGCTTCGCTCAATTTTATGGCGCTGTTCGGGCATACCGGCAGTGGCAAATGCACCCTCGCCAGTTTATTGATGGGCTATTACCCTCTAACGGAAGGGGAGATCCGCCTTGATGGTCTTCCATTAATGTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGTAATGGTTCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACTTGACGCTg                                                                                                                                                                                                           <  1:289335/248‑1 (MQ=255)
                                        tGGTGCCAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTgggg                                                                                                                                                                                                        >  1:117607/1‑251 (MQ=255)
                                        gggTGCTAAAGACAATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCGGGTCGGGCATACCGGCAGGGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGAGCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTgggg                                                                                                                                                                                                        <  2:26173/250‑1 (MQ=255)
                                             cTAAAGAACATTAATCTCTCTGTGCCTTCGCGCACTTTTGTGGCGATGGTCGGACATACCGGCAGTGGCAAAAGCACCCTCGACAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGCTGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGAGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATCCCTTCCTCGACAACGTGACGCTGGGGCgggg                                                                                                                                                                                                   >  1:64208/1‑250 (MQ=255)
                                               aaaGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGa                                                                                                                                                                                                   <  1:99694/249‑1 (MQ=255)
                                               aaaGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACACTCGCCAGTTTATTGCTGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGCTACCTTCCTCGCCAACGGGGCGCTGGGGCGGGatg                                                                                                                                                                                                 >  2:115303/1‑250 (MQ=255)
                                                 aGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATAtc                                                                                                                                                                                               >  1:341117/1‑251 (MQ=255)
                                                     cATTAATCTCTCTGTGCCGTCGCGCAATTTTGGGGCGCGGGTCGGGCATACCGGCAGTGGCAAAAGCCCCCTCGCCAGTTTATTGATGGGCTATGACCCTCGAACTGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTTACAGCGCGCGGCGCCAGGGCGTGTGATTGGTGCAGTAACATCCGGTGGTGCTGCCGGATACCTTCCTCGCCTACGTGATTCTGGGGCGGGATATCTCCg                                                                                                                                                                                           <  2:227022/251‑1 (MQ=255)
                                                      aTTAATCTCTCTGTGCCTGCGCGCAATTTTGTGGCGCTGTTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCga                                                                                                                                                                                          <  1:38589/251‑1 (MQ=255)
                                                        tAATCTCTCTGTGCCTTCGCGCACTTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCACAAGCACCCTCGCCCGATTATTGATGGGCTATTACCCGCTAACGGAAGGAGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCACCGAGACGCTGGGGCGGGatcgcaccga                                                                                                                                                                                          >  2:126498/1‑241 (MQ=255)
                                                         aaTCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTTATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATTTCCgaag                                                                                                                                                                                        <  2:155482/250‑1 (MQ=255)
                                                         aaTCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCgc                                                                                                                                                                                          >  1:403507/1‑247 (MQ=255)
                                                            ctctctGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTCGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAAcg                                                                                                                                                                                    >  1:158125/1‑251 (MQ=255)
                                                                  gtgCCTTCGCGCAATTTGGTGGCGCTGGTCTGGCATACCTGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGGCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCTTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGAGACCTTCCTCGCCATCGTGACGCTGGGGCGGGATATCTCCGAATAACGCGt                                                                                                                                                                                 <  1:374752/248‑1 (MQ=255)
                                                                                ttttGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCATTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAATAACGCGTCTGGCAGGCGCTGGaa                                                                                                                                                                 <  2:403507/250‑1 (MQ=255)
                                                                                 tgtgtgGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCGAACGGAAGGTGAGATGCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGGCACAGCGCGCTGCGCCAGGGCGTGGCAATGTTGCAGTAAGATCCGGTTGTGCTGGCGGAGTCGTGCCTCGCCAACGGGACGCGGGGTCGGGATAAGTCCGAAGATCGCGTCTGGCAGGCCCTGGaaa                                                                                                                                                                <  2:39221/248‑1 (MQ=255)
                                                                                            ggTCGGGCATCCCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGGGAGATTCGCCTTGAGGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGCAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCGGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACt                                                                                                                                                      <  1:425767/249‑1 (MQ=255)
                                                                                            ggTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACt                                                                                                                                                      <  1:115303/249‑1 (MQ=255)
                                                                                             gTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGc                                                                                                                                                   >  1:77780/1‑251 (MQ=255)
                                                                                                 ggCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCGGGAAACCGTGCAACTGGCGGAg                                                                                                                                               <  1:332773/251‑1 (MQ=255)
                                                                                                        cGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGgcgc                                                                                                                                       <  1:184152/251‑1 (MQ=255)
                                                                                                                   aaaGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGAGTCTGGCAGGCGCAGGAAACCGTGCCACTGGCGGAGCTGGCGCGTAGCATGAGc                                                                                                                            >  1:278578/1‑251 (MQ=255)
                                                                                                                         cccTCGCCATTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGt                                                                                                                      <  1:1546/251‑1 (MQ=255)
                                                                                                                              gCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGCAGAACGCGTCTGGCCGGCGCTGGAAACCGTGCAACTGGCGGGGCTGGCGCGTAGCATGGGCGGCGGTAttt                                                                                                                  >  2:24480/1‑250 (MQ=255)
                                                                                                                                 aGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGACTGCGCGGAGCCAGAGCGACGGTAttttaac                                                                                                              >  2:306863/1‑247 (MQ=255)
                                                                                                                                     tATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACAcgcc                                                                                                          <  1:495/250‑1 (MQ=255)
                                                                                                                                       ttGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGCAGCATGAGCGACGGTATTTacaacccgct                                                                                                        >  1:153455/1‑244 (MQ=255)
                                                                                                                                        tGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAGGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGGTATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGGCGGTATTTACACGCCGCTg                                                                                                      >  1:211128/1‑251 (MQ=255)
                                                                                                                                                 aTTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGTCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACTCCGCTGGGCGAGCAg                                                                                             <  1:326261/251‑1 (MQ=255)
                                                                                                                                                   tACCCGCTAACGGAAGGTGAGATTCGCCTCGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGCTCCGGTGGTGCTGGCGGATACCTTCCTAGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGGAACCGTGCAAATGGCGGCGCTGGCGCGGAGACTGAGCGACGGGATTTACAAGCCGCTGGGCGAGCAggg                                                                                           >  2:88933/1‑251 (MQ=255)
                                                                                                                                                           aaCGGATGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCCTGAGCGACGGGATTTACACGCCGCTGGGCGGGCAGGGGGATAAtc                                                                                   >  1:127022/1‑251 (MQ=255)
                                                                                                                                                             cGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGACACCTTCCTAGCCAACGTGACGCTGGGGAGGGGTATCTCCGAAGAACGCGTCTGGCAGGCGCTGGGAACCGTGCAACTGGCGGAGCTGGCGCGTAGGCTGAGCGACGGTCTTTAACAGCCGCTGGGCGAGCAGGGGAATAAactc                                                                                 >  2:172967/1‑251 (MQ=255)
                                                                                                                                                                 aGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAg                                                                             <  1:355211/251‑1 (MQ=255)
                                                                                                                                                                    tGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCATGGGAATAATCTCTCAGTTg                                                                          >  1:26628/1‑251 (MQ=255)
                                                                                                                                                                           cGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCaaaa                                                                   <  2:212092/251‑1 (MQ=255)
                                                                                                                                                                              cTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCTAGGGCGTGGCAACGGTGCAGTAAGCTCCGGTGGCGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGGAACCGGGCACCTGGCGCAGCTGGCGCGTCGCATGAGCGGCGGTATGTCAACACCGCGTGGCGGGCCGGGGAATAAACTCTCCATTGGGCAAAAGcc                                                                >  2:71517/1‑250 (MQ=255)
                                                                                                                                                                                 gATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAAct                                                             >  1:167143/1‑251 (MQ=255)
                                                                                                                                                                                 gATGGTCGTCCATTAAGTTCGCAAAGTCCCAGCGCGCGGCGCCAGGGCGTGGCAATGGGGCAGTAAGAGCCGGGGGTGCTGGCGTATACCTTCCTCGCCAACGTGACGCTGGGGCGGTATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCGGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAAct                                                             <  2:147327/251‑1 (MQ=255)
                                                                                                                                                                                  atggctgttcattaatgtCGCTAAGTCCCAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCGGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAActg                                                            <  1:6880/234‑1 (MQ=255)
                                                                                                                                                                                        gtcCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGcc                                                      >  1:2706/1‑250 (MQ=255)
                                                                                                                                                                                         tcCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCAc                                                     >  1:165253/1‑251 (MQ=255)
                                                                                                                                                                                          cctttAAGTGCGCTAAGGCACCGCGCGCGGCGCCAGGTCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAAGAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCAc                                                     <  2:117607/247‑1 (MQ=255)
                                                                                                                                                                                          ccATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACt                                                    >  1:345591/1‑251 (MQ=255)
                                                                                                                                                                                           cATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGGGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTg                                                   >  1:388372/1‑251 (MQ=255)
                                                                                                                                                                                                    cgctgtgtaacatgtcgcggggcgaGGGCGTGGCAATGGTGAAGTAAAAGCCGGGGGTGCGGTCGGATACCTTCCGCGCCAACTGGACGCTGGGGCTGGAGAGCTCCAAAGAACGCGTTTGGCAGGCGCTGGAAACCGTGCAACTGGCGTAGCGGGGGCGTGGCATGAGCGACGTTTTTTACACGCCGCTGGGCGAGCAGGGGAATAATCGCTCAGTTGGGCAAAAGCAACTGCTGGCTCTGGCGCGCGt                                           <  2:252509/226‑1 (MQ=255)
                                                                                                                                                                                                      cTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGCGCAGTAAGATCCGGTGGCGCTGGCGGATACCTTACTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCCTGAGCGACGGAATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCGTGCt                                        >  1:251321/1‑251 (MQ=255)
                                                                                                                                                                                                        aaGTCACAGCGCGCTGGGCAAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCGCGCCAGCGGGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCGCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCGTGCTg                                       <  1:245973/250‑1 (MQ=255)
                                                                                                                                                                                                          gTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATCCCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCAGGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCGTACTGGTc                                    <  1:116216/251‑1 (MQ=255)
                                                                                                                                                                                                           gcacggcgcgtTGCGCCGGGGCGTGGCAATGGGGCTGTAAGTTCCGGTGGTGCTGGCGGATACCTTCCGCGCCAACGGGACGCGGGGGCGGGATATCTCCGAAGACCGCGTTTGGAAGGCGCTGGAAACCGTGCATCTGGCGGAGCTGGCGCGTAGCATGAGCGACTGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCTTGCTGGTCg                                   <  2:380506/246‑1 (MQ=255)
                                                                                                                                                                                                             acaGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCGTGCTGGTCg                                   >  1:235784/1‑249 (MQ=255)
                                                                                                                                                                                                                       gtgactgggCGTGGCATTGGGGCAGTAGTATCCGGTGGTGCTGGCGGATGCCTTCCTCGCCATAGTGCAGCGGGGGCGGGTTAGCTCCGAAGAGCGCGTCTGGCAGGCGCTGGAAACCGGGCAACTGGCGGAGCTGGCGCGTAGCTGGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAAGCTCTCAGTTGGGCAAAAGGAACTGCTGGCACAGGCGCGCGTGCTTGACGATACACCGCaaa                      <  2:52282/245‑1 (MQ=255)
                                                                                                                                                                                                                        cgcCAGGGCGTGGGAATGTTGCAGTAAGATCCGGTGGTGCTGGCGGATCCCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAAAGCGTCTGGAAGGCGCTGGAAAACGTGCCACTGGCGGGGGTGGCGCGGAGCATGAGCGGAGGTATTTTCACGACGCGGGGCGAGGAGGGGAATAATCACTCCGTTGGGCAAAAACACCTGCTGGCACTTGCTCTCGTGCTGGTCGGAGCGCCGCaaa                      >  2:233555/1‑250 (MQ=255)
                                                                                                                                                                                                                             gggCGTGGCAATGGTGCAGTAAGATCCGGGGGTGCTGGCGGATCCCTTCCTCGCCAACGGGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCGGGCAGGCGCGGGAAACCGTGCAACGGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCGTGCTGGTCGAGACGCCGCAAATCCTGa                <  2:235784/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ggCAATGGTGCAGTAAGATCCGGTGGTGCTGGCGGCTACCTTCCTCGCCAACGTGACGCTGTGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGAATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTTCTGTCACTGGCGCGCGTGCTGGTCGAGACGCCGCAAATCCTGATCCTtg          <  2:24710/251‑1 (MQ=255)
                                                                                                                                                                                                                                     caagggTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCGGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCGTGCTGGTCGAGACGCCGCAAATCCTGATCCTtgat        <  1:153473/247‑1 (MQ=255)
                                                                                                                                                                                                                                       aTGGTGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGCGACGCTGGGGCGGGATATCTCCGAAGAACCCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCGTGCTGGTCGAGACGCCGCAAATCCTGATCCTtgatgg      >  1:182775/1‑250 (MQ=255)
                                                                                                                                                                                                                                           tGCAGTAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCGTGCTGGTCGAGACGCCGCAAATCCTGATCCTTGATGAGGc   >  1:30928/1‑250 (MQ=255)
                                                                                                                                                                                                                                            gcagtaAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCGGGGGCGGGATATCTCCGAAGAACGCGCCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCGTGCTGGTCGAGACGCCGCAAATCCTGATCCTTGATGAGGCa  <  1:314621/250‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        
GAAGTCGATAACGTGTCATTTGCTTATCGCGATGACAATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGGCTATTACCCGCTAACGGAAGGTGAGATTCGCCTTGATGGTCGTCCATTAAGTTCGCTAAGTCACAGCGCGCTGCGCCAGGGCGTGGCAATGGTGCAGCAAGATCCGGTGGTGCTGGCGGATACCTTCCTCGCCAACGTGACGCTGGGGCGGGATATCTCCGAAGAACGCGTCTGGCAGGCGCTGGAAACCGTGCAACTGGCGGAGCTGGCGCGTAGCATGAGCGACGGTATTTACACGCCGCTGGGCGAGCAGGGGAATAATCTCTCAGTTGGGCAAAAGCAACTGCTGGCACTGGCGCGCGTGCTGGTCGAGACGCCGCAAATCCTGATCCTTGATGAGGCA  >  NC_000913/471659‑472144

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