Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,438,536 A→C 5.5% L750R (CTG→CGG)  ydbK ← putative pyruvate‑flavodoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,438,5360AC 5.5% 707.1 / 26.0 218L750R (CTG→CGG) ydbKputative pyruvate‑flavodoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (104/102);  new base C (4/8);  total (108/110)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

Only 100 of 224 total aligned reads displayed.

GTTCCAGTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTG  >  NC_000913/1438387‑1438684
                                                                                                                                                     |                                                                                                                                                    
  tCCAGTTTGCTACGGTCGATTTCTGGCGGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAgg                                                                                                                                                    <  1:672024/150‑1 (MQ=255)
   ccAGTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGt                                                                                                                                                   >  2:290802/1‑150 (MQ=255)
    cAGTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGtt                                                                                                                                                  >  2:1717775/1‑150 (MQ=255)
      gTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTAc                                                                                                                                                >  2:1498258/1‑150 (MQ=255)
      gTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTAc                                                                                                                                                <  1:1370731/150‑1 (MQ=255)
      gTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTAc                                                                                                                                                <  2:103823/150‑1 (MQ=255)
           cTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCg                                                                                                                                           >  1:77393/1‑150 (MQ=255)
                gTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcg                                                                                                                                      <  1:425062/150‑2 (MQ=255)
                gTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCa                                                                                                                                      <  2:1940262/150‑1 (MQ=255)
                gTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCa                                                                                                                                      <  1:271734/150‑1 (MQ=255)
                   gATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATc                                                                                                                                   >  1:1583815/1‑150 (MQ=255)
                   gATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATc                                                                                                                                   <  2:508398/150‑1 (MQ=255)
                    aTTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgtcga                                                                                                                                  <  1:525428/150‑6 (MQ=255)
                    aTTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTGTAACCTGTGcgtcga                                                                                                                                  >  1:641783/1‑145 (MQ=255)
                    aTTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCt                                                                                                                                  <  1:604372/150‑1 (MQ=255)
                       tCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTcc                                                                                                                               <  1:143498/150‑1 (MQ=255)
                       tCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTcc                                                                                                                               <  1:962731/150‑1 (MQ=255)
                       tCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTcc                                                                                                                               <  1:1167425/150‑1 (MQ=255)
                        cTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCg                                                                                                                              >  1:1651398/1‑150 (MQ=255)
                        cTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCg                                                                                                                              <  2:162587/150‑1 (MQ=255)
                         tGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTg                                                                                                                                        <  1:1331729/139‑1 (MQ=255)
                         tGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCgg                                                                                                                             >  1:492210/1‑150 (MQ=255)
                          ggCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGt                                                                                                                            >  1:339266/1‑150 (MQ=255)
                             aGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGcc                                                                                                                         <  2:169104/150‑1 (MQ=255)
                             aGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGcc                                                                                                                         <  1:365759/150‑1 (MQ=255)
                              ggTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGGCGGTCAAACTTCTACTCAGAGGTTACAACCGGTGCAATCTTCCGGTgtta                                                                                                                        <  2:1389196/150‑4 (MQ=255)
                              ggTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCa                                                                                                                        >  1:720275/1‑150 (MQ=255)
                                ttGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCAcc                                                                                                                      >  1:1834222/1‑150 (MQ=255)
                                ttGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCAcc                                                                                                                      <  1:909566/150‑1 (MQ=255)
                                  gAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTg                                                                                                                    >  1:887126/1‑150 (MQ=255)
                                   aGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGc                                                                                                                                       <  2:1201419/130‑1 (MQ=255)
                                    ggAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCa                                                                                                                  >  2:783526/1‑150 (MQ=255)
                                     gaagaaATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCaa                                                                                                                 >  2:886266/1‑150 (MQ=255)
                                        gaaATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGAc                                                                                                              <  2:448283/150‑1 (MQ=255)
                                         aaaTCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACa                                                                                                             <  2:368450/150‑1 (MQ=255)
                                          aaTCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACat                                                                                                            <  2:1487924/150‑1 (MQ=255)
                                               tAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTc                                                                                                       <  2:708358/150‑1 (MQ=255)
                                                aaTTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGCCGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCt                                                                                                      >  2:1002041/1‑150 (MQ=255)
                                                 aTTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTg                                                                                                     <  1:290802/150‑1 (MQ=255)
                                                 aTTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTg                                                                                                     <  1:873159/150‑1 (MQ=255)
                                                     aTTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCcg                                                                                                 >  1:40089/1‑150 (MQ=255)
                                                     aTTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCcg                                                                                                 >  2:190368/1‑150 (MQ=255)
                                                         tctcTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCa                                                                                             >  2:1960654/1‑150 (MQ=255)
                                                          ctctTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTAc                                                                                                                                                >  1:1487972/1‑98 (MQ=255)
                                                           tcttctTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCata                                                                                           >  2:536824/1‑150 (MQ=255)
                                                           tcttctTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCata                                                                                           <  2:465307/150‑1 (MQ=255)
                                                           tcttctTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCAGGGCAAACTTAGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCata                                                                                           >  2:1356030/1‑150 (MQ=255)
                                                            cttcttCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCatat                                                                                          >  2:1875498/1‑150 (MQ=255)
                                                              tcttcGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTg                                                                                                                                        >  2:1750154/1‑102 (MQ=255)
                                                              tcttcGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCa                                                                                        >  1:1635682/1‑150 (MQ=255)
                                                                 tcGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACat                                                                                                            >  1:438710/1‑127 (MQ=255)
                                                                     cATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGAtt                                                                                 <  2:339266/150‑1 (MQ=255)
                                                                        gTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTca                                                                              >  1:1346720/1‑150 (MQ=255)
                                                                        gTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTca                                                                              <  1:483763/150‑1 (MQ=255)
                                                                          tCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTcaca                                                                            >  1:1252496/1‑150 (MQ=255)
                                                                           ccAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACAt                                                                           >  1:1279529/1‑150 (MQ=255)
                                                                            cAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATc                                                                          <  1:798537/150‑1 (MQ=255)
                                                                             aGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATcc                                                                         <  1:1450541/150‑1 (MQ=255)
                                                                              ggCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAg                                                                                                                <  2:1091092/110‑1 (MQ=255)
                                                                               gCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAg                                                                       >  1:83187/1‑150 (MQ=255)
                                                                                  agaCATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGa                                                                    <  2:1120461/150‑1 (MQ=255)
                                                                                    aCATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAAt                                                                  >  2:1820951/1‑150 (MQ=255)
                                                                                    aCATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAAt                                                                  >  1:1085232/1‑150 (MQ=255)
                                                                                      atcatATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGc                                                                <  2:434192/150‑1 (MQ=255)
                                                                                          tatTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGc                                                            <  2:1279529/150‑1 (MQ=255)
                                                                                           atTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCt                                                           >  2:16611/1‑150 (MQ=255)
                                                                                              gATGGCTTTAATCTCTGGATTCTGCCGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATACAGCGAATGCAGGCTGGc                                                        >  2:590012/1‑150 (MQ=255)
                                                                                                   gtttAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAgggg                                                   <  2:725674/149‑1 (MQ=255)
                                                                                                   cTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAgggg                                                   >  1:1202732/1‑150 (MQ=255)
                                                                                                   cTTTAATAACTGGATTCTAACAGTCTTTCTCAGTGAAAAATTCGACGCACAGCTTAAAAAATATGAAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCTGGCTGGCAgggg                                                   >  2:814493/1‑150 (MQ=255)
                                                                                                    tttAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGc                                                  >  2:275956/1‑150 (MQ=255)
                                                                                                    tttAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGc                                                  <  2:1559235/150‑1 (MQ=255)
                                                                                                     ttAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCg                                                 <  1:1554023/150‑1 (MQ=255)
                                                                                                     ttAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCg                                                 <  1:229239/150‑1 (MQ=255)
                                                                                                      tAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGt                                                >  1:466961/1‑150 (MQ=255)
                                                                                                         tctcTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGtttt                                             <  1:703366/150‑1 (MQ=255)
                                                                                                          ctctGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTc                                            >  2:1766425/1‑150 (MQ=255)
                                                                                                            ctGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCa                                          <  2:40089/150‑1 (MQ=255)
                                                                                                             tGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCAt                                         <  1:1147071/150‑1 (MQ=255)
                                                                                                                  tCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCtt                                    >  2:562039/1‑150 (MQ=255)
                                                                                                                     gACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAg                                 <  1:1680337/150‑1 (MQ=255)
                                                                                                                       cGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGACACCTGCAAGACATATTTCTGCCCGCGCATATAACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAggc                               >  2:1238852/1‑150 (MQ=255)
                                                                                                                        ggTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATg                                                                 >  2:225153/1‑115 (MQ=255)
                                                                                                                        ggTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAggcg                              <  2:1272888/150‑1 (MQ=255)
                                                                                                                        ggTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAggcg                              >  1:1966998/1‑150 (MQ=255)
                                                                                                                         gTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAggcgg                             <  2:93880/150‑1 (MQ=255)
                                                                                                                          tCTTTCGCCGGGCAAACTTCGACGCAGAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAggctgc                            <  1:1102314/150‑1 (MQ=255)
                                                                                                                          tCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAggcggc                            <  2:1150084/150‑1 (MQ=255)
                                                                                                                                 ccGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACtt                     <  2:1187409/150‑1 (MQ=255)
                                                                                                                                      cAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCCCCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTGTTgcg                >  2:1783409/1‑150 (MQ=255)
                                                                                                                                      cAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTgcg                <  1:566833/150‑1 (MQ=255)
                                                                                                                                        aaCTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGcgcg              <  1:1757560/150‑1 (MQ=255)
                                                                                                                                            agattgCACCGGTTGTAACCTGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGTAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATa          >  1:1016562/6‑150 (MQ=255)
                                                                                                                                            agattgCACCGGTTGTAACCTGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATa          >  1:178111/6‑150 (MQ=255)
                                                                                                                                             cGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAg         >  2:65832/1‑150 (MQ=255)
                                                                                                                                              gACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGc        <  1:1212221/150‑1 (MQ=255)
                                                                                                                                              gACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGc        >  1:119630/1‑150 (MQ=255)
                                                                                                                                                cGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTg      >  2:256493/1‑150 (MQ=255)
                                                                                                                                                  cacaGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAg    <  1:174937/150‑1 (MQ=255)
                                                                                                                                                  cacaGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAg    <  2:1443494/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTCCAGTTTGCTACGGTCGATTTCTGGCAGGTTGAGGAAGAAATCGTAATTGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTG  >  NC_000913/1438387‑1438684

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