Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,748,110 C→A 8.2% L343M (CTG→ATG)  lyxK → L‑xylulose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,748,1100CA 8.2% 252.5 / 12.2 93L343M (CTG→ATG) lyxKL‑xylulose kinase
Reads supporting (aligned to +/- strand):  ref base C (40/45);  new base A (2/6);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

GTGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTACT  >  NC_000913/3747961‑3748259
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                                            >  1:766435/1‑150 (MQ=255)
           gTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgcagc                                                                                                                                            >  2:1397263/1‑150 (MQ=255)
           gTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGgccgca                                                                                                                                           >  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:897896/1‑150 (MQ=255)
                          ccTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCgag                                                                                                                             <  2:558165/150‑1 (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                                                                                                          >  1:478414/1‑150 (MQ=255)
                                             tCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACg                                                                                                          <  2:200803/150‑1 (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)
                                                                                                                                                    ccTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTAc   >  2:1049687/1‑150 (MQ=255)
                                                                                                                                                     cTGTCCGGCAGCAACGCCGGACTCGAGATGTCCAGTGGTTTCTACGGGATGCAGGCCATTCACTCCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGCGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTACt  <  2:1596882/150‑1 (MQ=255)
                                                                                                                                                     cTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTACt  >  1:1853773/1‑150 (MQ=255)
                                                                                                                                                     cTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTACt  <  1:1738667/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGAATTTATCGTTCACGAAGCCAGCCCTACCTCTTCCGGCAACCTCGAATGGTTTACCGCACAGTGGGGAGAAATCTCGTTTGATGAGATCAATCAGGCCGTTGCCAGCTTGCCGAAGGCTGGGGGCGATCTCTTTTTCCTGCCGTTCCTGTACGGCAGCAACGCCGGACTCGAGATGACCAGTGGTTTCTACGGGATGCAGGCCATTCACACCCGCGCGCACCTGTTGCAGGCCATCTATGAAGGGGTGGTGTTCAGCCATATGACCCACCTCAACCGAATGCGCGAACGTTTTACT  >  NC_000913/3747961‑3748259

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: