Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,708,128 C→T 5.2% L209L (TTG→TTA rseB ← anti‑sigma E factor, binds RseA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,708,1280CT 5.2% 245.9 / 2.3 77L209L (TTG→TTArseBanti‑sigma E factor, binds RseA
Reads supporting (aligned to +/- strand):  ref base C (42/31);  new base T (2/2);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

TAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCGCATCGGTAATTTCGATTCGGTGT  >  NC_000913/2707981‑2708273
                                                                                                                                                   |                                                                                                                                                 
taatCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCa                                                                                                                                                  <  6:15323/149‑1 (MQ=255)
             tAGAGACGTGATTCGATAGGCATGTGGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCTCTAAAACCCTGTGGCAACCACGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAAtt                                                                                                                                     <  4:37872/149‑1 (MQ=255)
               gagaCGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAAttt                                                                                                                                    <  8:122391/148‑1 (MQ=255)
                   cGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCtt                                                                                                                               <  6:241160/149‑1 (MQ=255)
                       aTTCGATATGCATGTTTTACATCGTCGTTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTTGCAACCAGGTTGTCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATATGcc                                                                                                                                 <  6:357253/143‑1 (MQ=255)
                       aTTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGATAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCCACGGCGGCAAATTTGcc                                                                                                                                 >  5:357253/1‑143 (MQ=255)
                        ttCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCa                                                                                                                          >  2:323800/1‑149 (MQ=255)
                             tAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCACGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTc                                                                                                                     >  4:262486/1‑149 (MQ=255)
                                gCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGc                                                                                                                  <  4:211688/149‑1 (MQ=255)
                                   tgttgtCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATa                                                                                                               <  7:335185/149‑1 (MQ=255)
                                    gttgtCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATAc                                                                                                              >  1:40376/1‑149 (MQ=255)
                                     ttgtCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATAc                                                                                                              <  8:156962/148‑1 (MQ=255)
                                       gtCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATActgc                                                                                                           >  2:286314/1‑149 (MQ=255)
                                         ccATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATActgctg                                                                                                         >  7:379426/1‑149 (MQ=255)
                                         ccATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATActgctg                                                                                                         >  1:187003/1‑149 (MQ=255)
                                           aTCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATActgctgct                                                                                                       >  3:128/1‑149 (MQ=255)
                                           aTCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATActgctgct                                                                                                       >  8:375294/1‑149 (MQ=255)
                                            tcgtcgGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGTAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATActgctgct                                                                                                       >  1:389840/1‑148 (MQ=255)
                                               tcgGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGata                                                                                                   >  3:133814/1‑149 (MQ=255)
                                                   tAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCt                                                                                               >  4:201214/1‑149 (MQ=255)
                                                     gCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGa                                                                                             <  3:300407/149‑1 (MQ=255)
                                                       ggACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCCACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGAtt                                                                                           >  7:202323/1‑149 (MQ=255)
                                                        gACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGTAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGAtt                                                                                           >  7:358574/1‑148 (MQ=255)
                                                        gACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATATAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGACAGCGTCTGCATACTGCTGCTGATATCCTGATTg                                                                                          >  1:159874/1‑149 (MQ=255)
                                                          cGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGAc                                                                                        <  6:69556/149‑1 (MQ=255)
                                                                   ctGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGACAGCACCGGCGGCAAATTTGCCTTCGCCAGCGTCTTCATACTGCTGCTGATATCCTGATTGCCGTTAAAAg                                                                                >  1:207793/1‑148 (MQ=255)
                                                                   ctGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTaaaa                                                                                 <  1:59046/147‑1 (MQ=255)
                                                                   ctGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAg                                                                                <  2:50273/148‑1 (MQ=255)
                                                                   ctGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAg                                                                                >  7:86751/1‑148 (MQ=255)
                                                                    tGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGc                                                                               >  2:150633/1‑148 (MQ=255)
                                                                    tGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGc                                                                               >  1:272689/1‑148 (MQ=255)
                                                                    tGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGc                                                                               >  5:250619/1‑148 (MQ=255)
                                                                      gAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGAc                                                                                        >  6:230540/1‑137 (MQ=255)
                                                                      gAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGAc                                                                                        <  5:230540/137‑1 (MQ=255)
                                                                      gAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCaa                                                                             <  6:250619/148‑1 (MQ=255)
                                                                      gAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAAt                                                                            >  3:197147/1‑149 (MQ=255)
                                                                      gAAACTTCGCTAAAACCCTGTGGCAACCAGGGTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTTCTGATATCCTGATTGACGTTAAAAGCAAt                                                                            >  6:311231/1‑149 (MQ=255)
                                                                               cTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGaaa                                                                   >  8:384944/1‑149 (MQ=255)
                                                                                tAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAAt                                                                  >  7:290686/1‑149 (MQ=255)
                                                                                tAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAAt                                                                  >  3:365443/1‑149 (MQ=255)
                                                                                  aaaCCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCCACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGATGCTGATATCCTGAGTGACGTTAAAAGCAATCACGCGAAATTg                                                                >  2:16351/1‑149 (MQ=255)
                                                                                  aaaCCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTATTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGCCGTTAAAAGCAATCACGCGAAAt                                                                  <  2:22861/147‑1 (MQ=255)
                                                                                   aaCCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGt                                                               <  1:286314/149‑1 (MQ=255)
                                                                                      ccTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTcc                                                            <  2:40376/149‑1 (MQ=255)
                                                                                      ccTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTc                                                             >  2:132349/1‑148 (MQ=255)
                                                                                               aCCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCa                                                   >  3:165366/1‑149 (MQ=255)
                                                                                                 cAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCAc                                                  >  2:170401/1‑148 (MQ=255)
                                                                                                  aGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGTAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCa                                                <  2:389840/149‑1 (MQ=255)
                                                                                                       ggCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATcgc                                            <  4:133814/148‑1 (MQ=255)
                                                                                                          gTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGAt                                         >  5:212330/1‑148 (MQ=255)
                                                                                                          gTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATc                                        <  1:95909/149‑1 (MQ=255)
                                                                                                            ccAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCAcgc                                                                       <  1:259795/116‑1 (MQ=255)
                                                                                                            ccAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCAcgc                                                                       >  2:259795/1‑116 (MQ=255)
                                                                                                             cAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAg                                     <  1:93326/149‑1 (MQ=255)
                                                                                                              aggtGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCTTCTGCATACTGCTGCTGATATCCTGATTGACGTTTAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAgg                                    <  3:182882/146‑1 (MQ=255)
                                                                                                              aGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACAATCGCGATCAAg                                     >  2:316032/1‑148 (MQ=255)
                                                                                                                   aaTTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGgaga                                 <  2:272689/147‑1 (MQ=255)
                                                                                                                    aTTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGATGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGAACAAGGAGATc                               >  8:331414/1‑148 (MQ=255)
                                                                                                                     tttAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCaa                             >  6:324705/1‑149 (MQ=255)
                                                                                                                      ttAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAAc                            >  4:108324/1‑149 (MQ=255)
                                                                                                                         gCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCg                         >  6:114310/1‑149 (MQ=255)
                                                                                                                         gCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACATTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAAccc                          <  3:201214/148‑1 (MQ=255)
                                                                                                                             tttCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCGCa                       <  1:150633/147‑1 (MQ=255)
                                                                                                                                  ccTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGCCGTTAAACGCAATCACGCGACATTGTTCCAGCGTGTCACCCTCGCGATCAAGGAGATCAACCCGCATCGGTaa                >  4:225059/1‑149 (MQ=255)
                                                                                                                                   cTACAGGAACAGAAAGTAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGTAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCGCATCGGTAAt               <  8:358574/149‑1 (MQ=255)
                                                                                                                                    tACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCGCATCGGTAAtt              >  1:322775/1‑149 (MQ=255)
                                                                                                                                     tcAGGAACAGAAATCAACGGAGGCTAATTTGCCTTCGCCAGCGTCTGCATACTGCTGGTGATATCCTGATTGACGTTTAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGATATCAAGCCGCATCGGTAAttt             <  3:33277/148‑1 (MQ=255)
                                                                                                                                     aCAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCGCATCGGTAAttt             >  7:225972/1‑149 (MQ=255)
                                                                                                                                       aGGAACAGAAAGCAACGTCCGCAAATTTTCCTTCTCCAGCTTCTGCATACTGCTTCTCATATAATGATTGAAGTAAAAAGCAATCACGCGAAATTGTTCCAGCGTATCACCATCTCCATCAAGGAGATCAACCACCATCGGGAAATTc            >  4:33277/1‑148 (MQ=255)
                                                                                                                                       aGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCGCATCGGTAATTTc            >  7:131037/1‑148 (MQ=255)
                                                                                                                                          aaCAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCTTACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCCCGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCGCATCGGTAATTTCGAtt        >  3:230292/1‑149 (MQ=255)
                                                                                                                                            cAGAAAGCAACGGCGGCAAATTTGCCGTCGCCAGCGGCTGCAGACTGCTGCGGATATCCTGCTTGACGGAAAACGCAATCACGCGAAATTGTTCCAGCGTTTCACCCTCGCGATCAAGGAGATCACCCCGCCTCGGTAATTTCGATTCg      <  2:207793/149‑1 (MQ=255)
                                                                                                                                             aGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCGCATCGGTAATTTCGATTCgg     <  4:197147/149‑1 (MQ=255)
                                                                                                                                              gAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCAGCATCGGTAATTTCGATTCGgt    >  1:87628/1‑149 (MQ=255)
                                                                                                                                               aaaGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCa                                                           <  8:321224/93‑1 (MQ=255)
                                                                                                                                               aaaGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCGCATCGGTAATTTCGATTCGgtg   <  8:47464/149‑1 (MQ=255)
                                                                                                                                                aaGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCGCATCGGTAATTTCGATTCGgtgt  <  1:323800/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGTGGCAACCAGGTTGGCGTCCAGCTGAATTTAGCTTTTTCACCTACAGGAACAGAAAGCAACGGCGGCAAATTTGCCTTCGCCAGCGTCTGCATACTGCTGCTGATATCCTGATTGACGTTAAAAGCAATCACGCGAAATTGTTCCAGCGTTTCACCATCGCGATCAAGGAGATCAACCCGCATCGGTAATTTCGATTCGGTGT  >  NC_000913/2707981‑2708273

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