Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,998,147 C→A 5.6% Y55* (TAC→TAA uvrD → ssDNA translocase and dsDNA helicase ‑ DNA helicase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,998,1470CA 5.6% 862.0 / 25.0 250Y55* (TAC→TAAuvrDssDNA translocase and dsDNA helicase ‑ DNA helicase II
Reads supporting (aligned to +/- strand):  ref base C (123/113);  new base A (6/8);  total (129/121)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

Only 100 of 253 total aligned reads displayed.

CTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGC  >  NC_000913/3997998‑3998296
                                                                                                                                                     |                                                                                                                                                     
cTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAc                                                                                                                                                       >  1:138333/1‑150 (MQ=255)
  gCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTc                                                                                                                                                     <  2:774874/150‑1 (MQ=255)
   cTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCCACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTTATGAGCCTGGGAAACTGCTCGCCATAcacc                                                                                                                                                    >  2:1560742/1‑147 (MQ=255)
           ccTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGc                                                                                                                                            <  1:943443/150‑1 (MQ=255)
           ccTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGc                                                                                                                                            >  1:390852/1‑150 (MQ=255)
            cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCg                                                                                                                                           <  1:1228373/150‑1 (MQ=255)
             ttAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCgg                                                                                                                                          >  2:1779510/1‑150 (MQ=255)
             ttAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCgg                                                                                                                                          <  2:1043644/150‑1 (MQ=255)
              tAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGGGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAATCGAGTATGGCgag                                                                                                                                         <  2:623217/150‑3 (MQ=255)
               aaTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAATCGAGTATGGCgagc                                                                                                                                        <  2:863862/150‑4 (MQ=255)
                 tGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAc                                                                                                                                      >  1:1312454/1‑150 (MQ=255)
                 tGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAATCGAGTATGGCgagcag                                                                                                                                      <  1:619997/150‑6 (MQ=255)
                  gACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACg                                                                                                                                     <  1:1831388/150‑1 (MQ=255)
                      aaCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTa                                                                                                                                 >  1:1111255/1‑150 (MQ=255)
                      aaCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTa                                                                                                                                 >  1:1215810/1‑150 (MQ=255)
                      aaCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCATCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTa                                                                                                                                 >  2:691315/1‑150 (MQ=255)
                       aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGGGCGGGCAGTGTTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTGGATTTTGGCGGTGACGTTTAc                                                                                                                                <  2:1464228/150‑1 (MQ=255)
                       aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTAc                                                                                                                                <  1:1511103/150‑1 (MQ=255)
                        cAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTAcc                                                                                                                               <  1:889672/150‑1 (MQ=255)
                             cgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTTCcaacaa                                                                                                                          >  1:819533/1‑150 (MQ=255)
                             cgAAGCGGTGGCCCCGCCACGCAGCAACCTTCTGGTGCAGGCGGGCCCGGGCAGTGGTAAGACGCGCGTACTGGTGCAACGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGAATATGGCGGTGACGTTTACcaacaa                                                                                                                          >  2:1204171/1‑150 (MQ=255)
                                  cGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcgg                                                                                                                     >  1:1022392/1‑150 (MQ=255)
                                      ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcggcggc                                                                                                                 <  1:61539/150‑1 (MQ=255)
                                        CCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcg                                                                                                                      <  2:417974‑M1/143‑1 (MQ=255)
                                        CCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcg                                                                                                                      >  1:417974‑M1/1‑143 (MQ=255)
                                        CCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGcggcggcgg                                                                                                               <  2:49372‑M1/150‑1 (MQ=255)
                                              cccgccagcaccagaaggttgctgcgtggcgcgGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGc                                                                                                         >  1:731321‑M2/34‑150 (MQ=255)
                                              CACGCAGCAACCTTCTGGTGCAGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATATTATAGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGc                                                                                                         <  2:1117191‑M1/150‑1 (MQ=255)
                                               aCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCg                                                                                                        <  2:257609/150‑1 (MQ=255)
                                               aCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCg                                                                                                        <  1:842771/150‑1 (MQ=255)
                                                 gcagcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGAATCGTATCGCCTGCTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTc                                                                                                      >  2:1492638/1‑150 (MQ=255)
                                                     caACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGTCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCg                                                                                                  >  1:510066/1‑150 (MQ=255)
                                                             tGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGc                                                                                          <  2:178152/150‑1 (MQ=255)
                                                              ggTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAATACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAGGCTGCGGCGGAGATGCGTCTTCGTATCGGGCa                                                                                         <  2:449527/150‑1 (MQ=255)
                                                              ggTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCa                                                                                         >  1:1373020/1‑150 (MQ=255)
                                                              ggTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCa                                                                                         <  1:691315/150‑1 (MQ=255)
                                                                 gCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGt                                                                                                       <  1:1452398/133‑1 (MQ=255)
                                                                   tGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGa                                                                                    >  1:578786/1‑150 (MQ=255)
                                                                    ggcgggcgCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAg                                                                                                                        >  1:1814745/1‑113 (MQ=255)
                                                                    ggcgggcgCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAg                                                                                                                        <  2:1814745/113‑1 (MQ=255)
                                                                      cgggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATgg                                                                                 >  1:310224/1‑150 (MQ=255)
                                                                           gcggGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACg                                                                            <  1:1615528‑M2/146‑1 (MQ=255)
                                                                           gcggGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACg                                                                            <  2:1970927‑M2/146‑1 (MQ=255)
                                                                             ggGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAg                                                                          <  2:1219667‑M2/148‑1 (MQ=255)
                                                                                 aGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAg                                                                      <  1:1560742/150‑1 (MQ=255)
                                                                                  gTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAGAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAgg                                                                     >  1:1283868/1‑150 (MQ=255)
                                                                                  gTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAgg                                                                     >  1:1160637/1‑150 (MQ=255)
                                                                                  gTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAgg                                                                     <  2:54472/150‑1 (MQ=255)
                                                                                    ggTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGc                                                                   >  1:995404/1‑150 (MQ=255)
                                                                                    ggTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGc                                                                   <  2:654402/150‑1 (MQ=255)
                                                                                     gTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCg                                                                  >  1:1963143/1‑150 (MQ=255)
                                                                                       aaGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGt                                                                <  2:1284784/150‑1 (MQ=255)
                                                                                        aGACGCGCGCACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTa                                                               >  1:1479074/1‑150 (MQ=255)
                                                                                             cgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCg                                                                  >  2:858040/1‑142 (MQ=255)
                                                                                              gcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTggg                                                         <  1:791126/150‑1 (MQ=255)
                                                                                              gcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTggg                                                         <  2:1090200/150‑1 (MQ=255)
                                                                                                  aCTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCgg                                                     <  2:808962/150‑1 (MQ=255)
                                                                                                    tGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCa                                                   <  2:1832537/150‑1 (MQ=255)
                                                                                                          aTCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTcc                                             >  1:1705016/1‑150 (MQ=255)
                                                                                                          aTCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTcc                                             <  2:337703/150‑1 (MQ=255)
                                                                                                          aTCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTcc                                             <  2:138333/150‑1 (MQ=255)
                                                                                                              tATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACgg                                         <  2:1408376/150‑1 (MQ=255)
                                                                                                               aTCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACggg                                        >  2:217783/1‑150 (MQ=255)
                                                                                                               aTCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACggg                                        >  2:374396/1‑150 (MQ=255)
                                                                                                                tCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGc                                       >  1:1752047/1‑150 (MQ=255)
                                                                                                                     tGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgcg                                  <  1:1873356/150‑1 (MQ=255)
                                                                                                                     tGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgcg                                  <  2:1963143/150‑1 (MQ=255)
                                                                                                                      ggTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGAGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGAGGTATGTGGGTCGGCACCTTCCACGGGCTGgcgc                                 >  1:220500/1‑150 (MQ=255)
                                                                                                                        ttGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTgg                                    >  1:861648/1‑145 (MQ=255)
                                                                                                                         tgatgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCAcc                              >  1:328164/1‑150 (MQ=255)
                                                                                                                          gatgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCg                             <  1:1120360/150‑1 (MQ=255)
                                                                                                                            tgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGtt                           >  1:1662046/1‑150 (MQ=255)
                                                                                                                            tgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGtt                           <  1:765705/150‑1 (MQ=255)
                                                                                                                              aGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtg                         >  1:1240262/1‑150 (MQ=255)
                                                                                                                               gCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgc                        <  2:310224/150‑1 (MQ=255)
                                                                                                                               gCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGACGCACCGTTtgc                        <  2:607044/150‑1 (MQ=255)
                                                                                                                                cGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgct                       >  2:763485/1‑150 (MQ=255)
                                                                                                                                 gTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgctg                      <  2:1229360/150‑1 (MQ=255)
                                                                                                                                  tGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgctgc                     >  2:1136209/1‑150 (MQ=255)
                                                                                                                                  tGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgctgc                     >  1:1870054/1‑150 (MQ=255)
                                                                                                                                  tGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGAGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgctgc                     >  2:1402872/1‑150 (MQ=255)
                                                                                                                                   ggAAAACTGCTCGCTATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcg                    <  2:1881553/150‑1 (MQ=255)
                                                                                                                                   ggAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcg                    >  2:1443508/1‑150 (MQ=255)
                                                                                                                                    gAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcgt                   >  1:713816/1‑150 (MQ=255)
                                                                                                                                     aaaaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCAcc                                                 <  1:106963/119‑1 (MQ=255)
                                                                                                                                       aattgttCGCAATACTCGTTTTTGGGGGTTACTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTgcgc               <  2:1312454/144‑1 (MQ=255)
                                                                                                                                       aaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGGCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCACGGGCTGGCGCACCGTTTGCTGCGTgcgc               <  1:1302961/150‑1 (MQ=255)
                                                                                                                                        aCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTgcgc               <  2:1705016/150‑1 (MQ=255)
                                                                                                                                         gtcaccGCCATAATCGAGTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCa              >  2:1604821/6‑150 (MQ=255)
                                                                                                                                         gtcaccGCCATAATCGAGTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCa              >  1:974728/6‑150 (MQ=255)
                                                                                                                                         cTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCa              >  2:1888817/1‑150 (MQ=255)
                                                                                                                                          tcaccGCCATAATCGAGTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGAGGTATGTGGGTCGGCCCCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCAc             >  2:1078885/5‑150 (MQ=255)
                                                                                                                                           gCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCAcc            >  1:480129/1‑150 (MQ=255)
                                                                                                                                            cTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCa           >  2:407742/1‑150 (MQ=255)
                                                                                                                                            cTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCa           <  2:713816/150‑1 (MQ=255)
                                                                                                                                             tCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCat          >  1:1538952/1‑150 (MQ=255)
                                                                                                                                             ccGCCATAATCGAGTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCat          >  2:112504/2‑150 (MQ=255)
                                                                                                                                              cGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCata         <  2:1283868/150‑1 (MQ=255)
                                                                                                                                                ccATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCAGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATg       >  1:1086758/1‑150 (MQ=255)
                                                                                                                                                    aCTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACg   <  2:1752047/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGC  >  NC_000913/3997998‑3998296

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