Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,748,108 T→A 8.1% F342Y (TTC→TAC)  lyxK → L‑xylulose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,748,1080TA 8.1% 244.9 / 7.4 94F342Y (TTC→TAC) lyxKL‑xylulose kinase
Reads supporting (aligned to +/- strand):  ref base T (41/45);  new base A (2/6);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TGGTGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAA  >  NC_000913/3747959‑3748250
                                                                                                                                                     |                                                                                                                                              
tGGTGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGGTTGATGAGATCAATCAGGCCGTTGCCAGCTCGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGtt                                                                                                                                                <  2:1380044/150‑1 (MQ=255)
  gTGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGTGGCGATCTCTTTTTCCTGCCGTTcc                                                                                                                                              >  2:1672551/1‑150 (MQ=255)
    gAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTg                                                                                                                                            <  2:1652867/150‑1 (MQ=255)
       tttATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTAc                                                                                                                                         <  1:678002/150‑1 (MQ=255)
         tATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACgg                                                                                                                                       >  2:847829/1‑150 (MQ=255)
          ttCGTTCTCGAAGCCAGCCCTACCTCTTCCGTCAAACTCGTATGGTTTACCGCACAGTGGTGAGAAATCTAGTTTGATTAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgc                                                                                                                                      <  2:1854306/149‑1 (MQ=255)
           tCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgca                                                                                                                                     <  1:833478/150‑1 (MQ=255)
            cGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgcag                                                                                                                                    >  1:407877/1‑150 (MQ=255)
            cGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTACAGGAACGgcag                                                                                                                                    <  1:224379/150‑1 (MQ=255)
             gTTCCCGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTGTGATGAGATCAATCAGGCAGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTACAGGAACGgcagg                                                                                                                                   <  2:1193820/150‑2 (MQ=255)
             gTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTTATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgcagc                                                                                                                                   >  2:345409/1‑150 (MQ=255)
             gTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgcagc                                                                                                                                   >  2:1397263/1‑150 (MQ=255)
             gTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgcagc                                                                                                                                   >  1:766435/1‑150 (MQ=255)
             gTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCT‑GGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgccgca                                                                                                                                  >  1:774892/1‑150 (MQ=255)
             gTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTAACGCACAGAGGGGAGAAATCTCGTTTGATTAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGCCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgcagc                                                                                                                                   >  1:1586423/1‑150 (MQ=255)
              ttCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTACAGGAACGgcagga                                                                                                                                  <  2:672862/150‑3 (MQ=255)
              aacACGAAGCCAGCCCTACCTCTTCCCCCAACCTCGAATGGTTTACCGCACATTGGGGAGAAAACTCGTTTGATGATATCAATCAGGGCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTATTTCCTGCCGTACAGGAACGgcagga                                                                                                                                  <  2:952132/148‑3 (MQ=255)
                cACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAACACTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgcagcagc                                                                                                                                >  2:903174/1‑148 (MQ=255)
                cACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCCAc                                                                                                                                >  1:1363360/1‑150 (MQ=255)
                cACGAAGCCAGCCCTACCTCATCCGGCAACCTCGAATGGTTTACCGCACAGGGGGGAGCAATCTCGTTTGATGAGATCAATCAGCCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTACAGGAACGGCAGGaaa                                                                                                                                <  1:1387993/150‑2 (MQ=255)
                  cGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGc                                                                                                                              <  1:481687/150‑1 (MQ=255)
                   gAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTACAGGAACGGCAGGaaaaag                                                                                                                             <  2:1088645/150‑5 (MQ=255)
                   gAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAAGCTCGTGTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGcc                                                                                                                             <  1:11236/150‑1 (MQ=255)
                    aaGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGCGAAATCTCGTTTGATGAGATCAATCAGGCCGGTGCCAGCGTGCCGAAGGCGGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCg                                                                                                                            <  1:1119701/150‑1 (MQ=255)
                    aaGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCg                                                                                                                            <  1:1334855/150‑1 (MQ=255)
                       ccagccCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCCGCAACGCCGGAc                                                                                                                         >  1:913798/1‑150 (MQ=255)
                           ccgttccTCTTCCGGCATACTCGAATTGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATGAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCATGCCGTTCCTGTACGGCAGCAACGCCGGACTCga                                                                                                                     <  2:1842867/145‑1 (MQ=255)
                            ccTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCgag                                                                                                                    <  2:558165/150‑1 (MQ=255)
                            ccTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCgag                                                                                                                    >  2:897896/1‑150 (MQ=255)
                               aCCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATg                                                                                                                 <  2:1050780/150‑1 (MQ=255)
                               aCCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGGGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCATCTTGCCGAAGGCTGGGGTCGATCTCTTTTTCCTGCCGTTCCTGTACGTCAGCAACGCCGGACTCGAGATg                                                                                                                 <  1:1527171/150‑1 (MQ=255)
                                 ctctTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGCTCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGAc                                                                                                               >  1:1063344/1‑150 (MQ=255)
                                     tCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGt                                                                                                           >  2:1232966/1‑150 (MQ=255)
                                         gCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGtt                                                                                                       >  1:15567/1‑150 (MQ=255)
                                          cAACCTCGAATGGTTTACCGCACAGTGGGGAGAAAGCTCGGGTGATGAGATCAATCAGGCCGTTGCCAGCGTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGttt                                                                                                      <  1:286989/150‑1 (MQ=255)
                                           aaCCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGTGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTc                                                                                                     <  1:471745/150‑1 (MQ=255)
                                            aCCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCt                                                                                                    <  1:1766738/150‑1 (MQ=255)
                                             ccTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTa                                                                                                   >  2:1317973/1‑150 (MQ=255)
                                             ccTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTCGTTTCTa                                                                                                   >  2:224606/1‑150 (MQ=255)
                                               tCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACg                                                                                                 <  2:200803/150‑1 (MQ=255)
                                               tCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACg                                                                                                 >  1:478414/1‑150 (MQ=255)
                                                 gAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTAGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACggg                                                                                               <  2:590152/150‑1 (MQ=255)
                                                  aaTGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGa                                                                                              <  1:197909/150‑1 (MQ=255)
                                                     ggTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGc                                                                                           <  2:1121455/150‑1 (MQ=255)
                                                      gTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCa                                                                                          <  2:1494287/150‑1 (MQ=255)
                                                        ttACCGCACAGTGGGGAGATCTCTCGTTTGATGAGATCAATCAGGCCGGTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAgg                                                                                        <  1:141878/150‑1 (MQ=255)
                                                        ttACCGCACAGTGGGGAGAAACCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAgg                                                                                        >  1:773745/1‑150 (MQ=255)
                                                         tttaccacaGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCCGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGc                                                                                       <  2:1707374/145‑1 (MQ=255)
                                                         tACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGc                                                                                       <  1:146308/150‑1 (MQ=255)
                                                          aCCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGCCCGTTGACACCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGTACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGcc                                                                                      >  1:987278/1‑150 (MQ=255)
                                                            cGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCAt                                                                                    >  1:780466/1‑150 (MQ=255)
                                                            cGCACAGTGGGGAGAAATCTCGTTTGATGAGAACAATCACGCCGATGCCAGCTTGCCGAAGGCTGGCGGCGATATCATTTTCCTGCCGAACCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTATCTACGGGATGCAGGCCAt                                                                                    >  1:1259803/1‑150 (MQ=255)
                                                                 aggggggAGAAATCTCGGTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTcaca                                                                               <  2:930476/147‑1 (MQ=255)
                                                                 aGTGGGGAGAAATCTCGGTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACGCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTcaca                                                                               <  2:223669/150‑1 (MQ=255)
                                                                  gTGGGGAGAAATCTCGATTGATGAGATCAATCAGGCCGGAGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCacac                                                                              <  2:305405/150‑1 (MQ=255)
                                                                   gggggAGAAATCTCGGTAGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACAcc                                                                             <  1:679245/149‑1 (MQ=255)
                                                                       gagaAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCcgcg                                                                         <  1:1601099/150‑1 (MQ=255)
                                                                         gaAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACCCCcgcgcg                                                                       >  2:1790872/1‑150 (MQ=255)
                                                                         gaAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCcgcgcg                                                                       <  1:517395/150‑1 (MQ=255)
                                                                          aaaTCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCgcgcgc                                                                      <  2:890708/150‑1 (MQ=255)
                                                                             actcGTTGGAAGAGATGAATCAGGCCGTGGCCAGCTTGCCGAAGGCGTGGGGCGATATCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCAcc                                                                   <  2:252488/149‑1 (MQ=255)
                                                                              ctcGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCt                                                                  >  2:1403038/1‑150 (MQ=255)
                                                                                      atgaGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAgg                                                          <  1:1545100/150‑1 (MQ=255)
                                                                                        gagaTCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGCCCAGTGGTTTCTACGGGATGCAGGCCCTTCACACCCGCGCGCCCCTGTTGCAGGcc                                                        >  2:768396/1‑150 (MQ=255)
                                                                                        gagaTCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCCCCTGTTGCAGGcc                                                        >  2:1634088/1‑150 (MQ=255)
                                                                                         agaTCAATCAGGCCGTTGCCAGCTTGCCGTCGGCTGGGGGCGATCGCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCa                                                       <  1:220746/150‑1 (MQ=255)
                                                                                          gaTCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGGTGCAGGCCAt                                                      >  1:27774/1‑150 (MQ=255)
                                                                                            tcaatcaGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGCTGCCCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCAGCt                                                    >  1:1634079/1‑150 (MQ=255)
                                                                                              aatcaGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTAt                                                  <  2:766435/150‑1 (MQ=255)
                                                                                                  aGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAg                                              <  1:1397263/150‑1 (MQ=255)
                                                                                                   ggCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCTAGATGACCAGTGGTGTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCAGCTATGAAgg                                             >  2:1226585/1‑150 (MQ=255)
                                                                                                      cGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACCACCGCGCGCACCTGTTGCAGGCCATCTATGAAGGggt                                          >  2:1312681/1‑150 (MQ=255)
                                                                                                       gTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGggtg                                         <  1:345409/150‑1 (MQ=255)
                                                                                                         tGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTAATGTACGGCAGCAACGCCGGACTCGAGATGTCAAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCCCCTGTTGCAGGCCATCTATGAAGGcgtggg                                       >  2:283337/1‑149 (MQ=255)
                                                                                                            cAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGACGGACTCGAGATGACCAGTTGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCGATCTATGAAGGGGTGGTGtt                                    >  1:538664/1‑150 (MQ=255)
                                                                                                             aGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTc                                   <  1:577043/150‑1 (MQ=255)
                                                                                                              gCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCCCCGGTTGCAGGCCATCTATGAAGGGGTGGTGTTCa                                  >  2:552846/1‑150 (MQ=255)
                                                                                                              gCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCa                                  >  2:601769/1‑150 (MQ=255)
                                                                                                                 tGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGcc                               >  2:1026392/1‑150 (MQ=255)
                                                                                                                  gCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCa                              >  1:205875/1‑150 (MQ=255)
                                                                                                                    cGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCata                            <  1:170928/150‑1 (MQ=255)
                                                                                                                    cGAAGGCTGGGGGCGAGCTCTTTTTCCTGCCGTTCCTGTACTGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCata                            <  2:500338/150‑1 (MQ=255)
                                                                                                                      caGGCTGGGGGCGCTCTCGATTTCCTGCCGTTCCTGTACGGCAGCAACCGCGGACTCGAGCTGACCAGTGGTTGCTACGGGAGGCAGGCCATGCTCACCCGGGCGATCCTGCTGCAGGCCATCTCTGAAGGGGTGGTGTTCAGCCATATg                          <  1:242338/149‑1 (MQ=255)
                                                                                                                       aGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGa                         <  1:1679320/150‑1 (MQ=255)
                                                                                                                          cTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGAccc                      >  2:391897/1‑150 (MQ=255)
                                                                                                                          cTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGAAGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGAccc                      <  1:1113077/150‑1 (MQ=255)
                                                                                                                          cTGGGGGCGATCTCTTTATCCTGCCGTTCCTGGACGGCAGCAACGCCTGACTCGAGATGACCAGTGGTTTCTACTTGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATTTGAccg                      <  2:548607/150‑2 (MQ=255)
                                                                                                                             ggggCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCAcc                   <  1:1317973/150‑1 (MQ=255)
                                                                                                                              gggCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACACCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCt                  >  2:808239/1‑150 (MQ=255)
                                                                                                                               ggCGATCTCTTTTTCCTGACGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTATCATCGGGATGCAGGCCATTCACACCCGCTCTCACCTGTTGCAGGACATCAATGAATGGGTGTTGATCCGCCATATGCCCCACCTc                 >  1:72113/1‑150 (MQ=255)
                                                                                                                                  gATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCCCCTCAAc              >  2:1643471/1‑150 (MQ=255)
                                                                                                                                    tctcTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCg            >  1:237708/1‑150 (MQ=255)
                                                                                                                                        gcgttgcTGCCGTACAGGAACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATg        >  1:617740/7‑150 (MQ=255)
                                                                                                                                            tCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGcgcg    <  2:407877/150‑1 (MQ=255)
                                                                                                                                              cTGCCGGTCCTGTACGGCAGTAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACCCCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTGCAGCCAAATGACCCACCGCAACCGAATGCGCGaa  >  1:1455010/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TGGTGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAA  >  NC_000913/3747959‑3748250

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