Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 931,298 C→T 7.3% A251A (GCG→GCA trxB ← thioredoxin reductase, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913931,2980CT 7.3% 647.9 / 27.1 192A251A (GCG→GCAtrxBthioredoxin reductase, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (79/99);  new base T (7/7);  total (86/106)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

CGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCA  >  NC_000913/931150‑931446
                                                                                                                                                    |                                                                                                                                                    
cGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCg                                                                                                                                                     >  1:318657/1‑150 (MQ=255)
  gCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCa                                                                                                                                                   >  2:2121055/1‑150 (MQ=255)
  gCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCa                                                                                                                                                   >  2:1736832/1‑150 (MQ=255)
  gCCGAAGTAATCGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCa                                                                                                                                                   <  2:2000195/150‑1 (MQ=255)
         tAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCg                                                                                                                                            >  1:365804/1‑150 (MQ=255)
         tAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCg                                                                                                                                            >  2:442066/1‑150 (MQ=255)
           aTGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCggg                                                                                                                                          <  2:1721100/150‑1 (MQ=255)
           aTGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCggg                                                                                                                                          <  1:729291/150‑1 (MQ=255)
            tGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGc                                                                                                                                         >  1:1848703/1‑150 (MQ=255)
                cTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCtgtg                                                                                                                                     >  1:1450940/1‑150 (MQ=255)
                  ggCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGAc                                                                                                                                   <  2:1921506/150‑1 (MQ=255)
                  ggCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGAc                                                                                                                                   <  2:2698970/150‑1 (MQ=255)
                   gCGAGAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATGACCATGACTACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGAcc                                                                                                                                  <  2:2424044/150‑1 (MQ=255)
                    cGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCg                                                                                                                                 <  1:2506644/150‑1 (MQ=255)
                      aTAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGACAATCGCAGTATTCGGGCTGTGACCGAt                                                                                                                               >  1:2735235/1‑150 (MQ=255)
                      aTAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGAt                                                                                                                               <  2:2885414/150‑1 (MQ=255)
                      aTAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGAt                                                                                                                               <  2:1914757/150‑1 (MQ=255)
                      aTAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGAt                                                                                                                               >  2:692895/1‑150 (MQ=255)
                      aTAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACCGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGAt                                                                                                                               >  1:2871579/1‑150 (MQ=255)
                          aTGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGca                                                                                                                           >  1:1142962/1‑150 (MQ=255)
                           tgtgATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGcaa                                                                                                                          >  2:1703970/1‑150 (MQ=255)
                               aTCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCaacaaa                                                                                                                      <  2:906040/150‑1 (MQ=255)
                                 ccATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCaacaaaca                                                                                                                    >  1:431857/1‑150 (MQ=255)
                                  cATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAg                                                                                                                   <  1:2309359/150‑1 (MQ=255)
                                    tCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGAc                                                                                                                 <  1:1925442/150‑1 (MQ=255)
                                      aCGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCg                                                                                                               <  1:328731/150‑1 (MQ=255)
                                        gTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGAGAGCAACAAACAGACCGGc                                                                                                             >  1:573165/1‑150 (MQ=255)
                                         tCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCa                                                                                                            <  1:2825444/150‑1 (MQ=255)
                                          cGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCaa                                                                                                           <  1:668224/150‑1 (MQ=255)
                                           gCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTTCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAAc                                                                                                          <  1:509260/150‑1 (MQ=255)
                                            ccTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACg                                                                                                         >  1:763743/1‑150 (MQ=255)
                                               gCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCACCGTCg                                                                                                      >  1:1652350/1‑150 (MQ=255)
                                                cGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGa                                                                                                     >  2:185898/1‑150 (MQ=255)
                                                   cAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTg                                                                                                  >  1:77311/1‑150 (MQ=255)
                                                     aaGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATTGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGAc                                                                                                >  1:1879012/1‑150 (MQ=255)
                                                     aaGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATTGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGAc                                                                                                >  1:801398/1‑150 (MQ=255)
                                                     aaGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATTGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGAc                                                                                                >  1:801466/1‑150 (MQ=255)
                                                     aaGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATTGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGAc                                                                                                >  1:1010806/1‑150 (MQ=255)
                                                        aCGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCg                                                                                             >  2:2259886/1‑150 (MQ=255)
                                                            cAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACATCGAGTCACTCGACGt                                                                                         <  2:1200594/150‑1 (MQ=255)
                                                             aGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGtt                                                                                        >  2:2768863/1‑150 (MQ=255)
                                                                aaTGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATc                                                                                     <  2:95397/150‑1 (MQ=255)
                                                                  tGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGc                                                                                   >  2:2597201/1‑150 (MQ=255)
                                                                   gCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCt                                                                                  >  1:1157189/1‑150 (MQ=255)
                                                                     tGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGt                                                                                <  1:2880342/150‑1 (MQ=255)
                                                                     tGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGt                                                                                <  2:494950/150‑1 (MQ=255)
                                                                       gTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGttt                                                                              >  2:887427/1‑150 (MQ=255)
                                                                          tGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGc                                                                           >  2:1515194/1‑150 (MQ=255)
                                                                             gTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTa                                                                        >  1:2491187/1‑150 (MQ=255)
                                                                                 cATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATcgc                                                                    <  1:753081/150‑1 (MQ=255)
                                                                                 cATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATcgc                                                                    >  2:1235023/1‑150 (MQ=255)
                                                                                   ttACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCgcgc                                                                  >  1:1086722/1‑150 (MQ=255)
                                                                                    tACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTTGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCa                                                                 <  2:1018984/150‑1 (MQ=255)
                                                                                    tACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCa                                                                 <  2:1019001/150‑1 (MQ=255)
                                                                                    tACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCa                                                                 <  2:303054/150‑1 (MQ=255)
                                                                                    tACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCa                                                                 <  1:2489157/150‑1 (MQ=255)
                                                                                    tACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCa                                                                 <  1:2585675/150‑1 (MQ=255)
                                                                                    tACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCa                                                                 <  2:1297870/150‑1 (MQ=255)
                                                                                     aCCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAg                                                                <  1:1235734/150‑1 (MQ=255)
                                                                                      ccATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGa                                                               <  1:1902024/150‑1 (MQ=255)
                                                                                        aTGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACg                                                             >  1:360421/1‑150 (MQ=255)
                                                                                         tGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGa                                                            >  2:949716/1‑150 (MQ=255)
                                                                                          gAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGaa                                                           <  1:185898/150‑1 (MQ=255)
                                                                                           aaTACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAAc                                                          >  1:2905799/1‑150 (MQ=255)
                                                                                               cccGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCa                                                      <  2:719867/150‑1 (MQ=255)
                                                                                               cccGACTGTACTTGGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCa                                                      <  2:573165/150‑1 (MQ=255)
                                                                                                ccGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAg                                                     >  1:707671/1‑150 (MQ=255)
                                                                                                   aCTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGa                                                  <  2:2905119/150‑1 (MQ=255)
                                                                                                   aCTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGa                                                  <  1:1736832/150‑1 (MQ=255)
                                                                                                   aCTGTAATTTGTTTTTGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGAAACGTCGTGTGACTCGATGTTATCGCTTTTTTGCGTATCGCGCAGACGAACGCCAGTGa                                                  <  1:2121055/150‑1 (MQ=255)
                                                                                                    cTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGac                                                 <  2:77303/150‑1 (MQ=255)
                                                                                                    cTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGac                                                 <  2:77311/150‑1 (MQ=255)
                                                                                                          tttGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCAt                                           >  1:182963/1‑150 (MQ=255)
                                                                                                          tttGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCAt                                           >  2:838918/1‑150 (MQ=255)
                                                                                                               tGTAGCCTTTTTCCTGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGAt                                      <  1:2792139/150‑1 (MQ=255)
                                                                                                               tGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGAt                                      <  1:2636868/150‑1 (MQ=255)
                                                                                                                gTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATc                                     <  2:1396168/150‑1 (MQ=255)
                                                                                                                  aGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGc                                   <  2:936749/150‑1 (MQ=255)
                                                                                                                  aGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGc                                   <  2:936804/150‑1 (MQ=255)
                                                                                                                  aGCCGTTTTCCAGTTCCAGCTGCCCGTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGc                                   <  2:847049/150‑1 (MQ=255)
                                                                                                                   gCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGcc                                  >  2:1413589/1‑150 (MQ=255)
                                                                                                                    ccGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCg                                 >  2:2356374/1‑150 (MQ=255)
                                                                                                                            cagttccagCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCActtc                         >  1:1584942/1‑150 (MQ=255)
                                                                                                                            cagttccagCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCActtc                         >  2:1768325/1‑150 (MQ=255)
                                                                                                                            cagttccagCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCActtc                         >  1:90788/1‑150 (MQ=255)
                                                                                                                             agttccagCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCActtct                        >  2:1211109/1‑150 (MQ=255)
                                                                                                                                  cagCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAg                   <  2:90790/150‑1 (MQ=255)
                                                                                                                                      tGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTa               <  1:1045520/150‑1 (MQ=255)
                                                                                                                                      tGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTa               <  1:1373664/150‑1 (MQ=255)
                                                                                                                                      tGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTa               <  1:254011/150‑1 (MQ=255)
                                                                                                                                      tGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTa               <  1:348972/150‑1 (MQ=255)
                                                                                                                                      tGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTa               <  2:1652350/150‑1 (MQ=255)
                                                                                                                                          cTTCGAAAATTGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGt           <  2:801358/150‑1 (MQ=255)
                                                                                                                                          cTTCGAAAATTGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGt           <  2:801447/150‑1 (MQ=255)
                                                                                                                                          cTTCGAAAATTGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGt           <  2:1879012/150‑1 (MQ=255)
                                                                                                                                            tCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTg         >  2:2558938/1‑150 (MQ=255)
                                                                                                                                            tCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTg         >  2:201293/1‑150 (MQ=255)
                                                                                                                                            tCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTg         >  1:827505/1‑150 (MQ=255)
                                                                                                                                                 aaTCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGtgtg    >  2:1344323/1‑150 (MQ=255)
                                                                                                                                                   tCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCa  <  1:2447763/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCA  >  NC_000913/931150‑931446

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