Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,998,152 T→G 6.3% I57S (ATT→AGT)  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,1520TG 6.3% 296.0 / 10.0 95I57S (ATT→AGT) uvrDssDNA translocase and dsDNA helicase ‑ DNA helicase II
Reads supporting (aligned to +/- strand):  ref base T (38/51);  new base G (2/4);  total (40/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCC  >  NC_000913/3998004‑3998297
                                                                                                                                                    |                                                                                                                                                 
gACAGCCTTAATGACACACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTGTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGCTGAGCGTGGAAAACTGCTCGCCATACTCGAtt                                                                                                                                                  <  2:470940/150‑1 (MQ=255)
gACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGAtt                                                                                                                                                  <  2:1204507/150‑1 (MQ=255)
  cAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAATCGAGTAt                                                                                                                                                <  2:465047/150‑1 (MQ=255)
    gCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTGTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATgg                                                                                                                                              >  1:759694/1‑150 (MQ=255)
    gCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATgg                                                                                                                                              >  2:1222399/1‑150 (MQ=255)
    gCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATgg                                                                                                                                              >  2:578706/1‑150 (MQ=255)
     ccTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACAGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGc                                                                                                                                             <  1:172900/150‑1 (MQ=255)
       ttAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCgg                                                                                                                                           <  1:1192043/150‑1 (MQ=255)
        tAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAATCGAGTATGGCgag                                                                                                                                          <  1:1068614/150‑3 (MQ=255)
         aaTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTg                                                                                                                                         <  1:1229737/150‑1 (MQ=255)
         aaTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTg                                                                                                                                         >  2:1132167/1‑150 (MQ=255)
         aaTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAATCGAGTATGGCgagc                                                                                                                                         <  2:1346654/150‑4 (MQ=255)
         aaTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAATCGAGTATGGCgagc                                                                                                                                         <  2:1197016/150‑4 (MQ=255)
         aaTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCCGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCAGGTTGAATAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTg                                                                                                                                         >  2:1612533/1‑150 (MQ=255)
          aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGa                                                                                                                                        <  2:840263/150‑1 (MQ=255)
          aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGa                                                                                                                                        <  2:834511/150‑1 (MQ=255)
             aCAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATTCTCGATTATGGCGGTGACGt                                                                                                                                     <  2:1548050/150‑1 (MQ=255)
                aaCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTa                                                                                                                                  <  1:809043/150‑1 (MQ=255)
                aaCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTa                                                                                                                                  >  2:841824/1‑150 (MQ=255)
                 aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTAc                                                                                                                                 <  1:211959/150‑1 (MQ=255)
                  cAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTAcc                                                                                                                                <  1:711837/150‑1 (MQ=255)
                    gcgcGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACcaa                                                                                                                              <  1:1606656/150‑1 (MQ=255)
                     cgcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACcaac                                                                                                                             <  2:1303951/150‑1 (MQ=255)
                       cgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACcaacaa                                                                                                                           >  1:787047/1‑150 (MQ=255)
                       cgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGATCGCCATACTCGATTATGGAGGTGACGTTTacaaacaa                                                                                                                           >  2:699716/1‑150 (MQ=255)
                       cgAAGCGGTGGCCGCGCCAAGCAGCAACATTCTGGTGCTGGCGGGCGAGGGCAGTGGTAAGACGCCCGTACAGGTTCATAGTATCGCCAGGTTGAAGAGAGTGGAAAACAGCTAGCCAAACACGATTATGTCGCTGACCTTTaaaaacca                                                                                                                           >  2:699722/1‑143 (MQ=255)
                         aaGCGGTGGCCGCGCCACGCAGCAAGCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACGAACAAAg                                                                                                                         <  2:622617/150‑1 (MQ=255)
                         aaGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAg                                                                                                                         <  2:362541/150‑1 (MQ=255)
                             ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcggc                                                                                                                     >  2:1584833/1‑150 (MQ=255)
                                ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcggcggc                                                                                                                  >  2:46765/1‑150 (MQ=255)
                                 gCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcggcggcg                                                                                                                 >  1:494782/1‑150 (MQ=255)
                                   cgcgCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTCGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGga                                                                                                               <  2:766980/150‑1 (MQ=255)
                                    gcgcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGgag                                                                                                              <  2:1567163/150‑1 (MQ=255)
                                         aCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGGTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCg                                                                                                         <  1:234640/150‑1 (MQ=255)
                                          cGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGt                                                                                                        <  1:195971/150‑1 (MQ=255)
                                            cagcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCa                                                                                                      >  1:1093906/1‑150 (MQ=255)
                                              gcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATc                                                                                                    >  2:189499/1‑150 (MQ=255)
                                              gcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATc                                                                                                    <  1:389763/150‑1 (MQ=255)
                                                     tCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCgg                                                                                             <  2:368911/150‑1 (MQ=255)
                                                     tCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCgg                                                                                             <  1:522650/150‑1 (MQ=255)
                                                      cTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCggg                                                                                            <  1:958480/150‑1 (MQ=255)
                                                      cTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCggg                                                                                            <  2:1376810/150‑1 (MQ=255)
                                                        ggTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCa                                                                                          <  2:1093906/150‑1 (MQ=255)
                                                              ggcgggcgCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGAt                                                                                    <  1:113040/150‑1 (MQ=255)
                                                               gcgggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATg                                                                                   <  1:1129530/150‑1 (MQ=255)
                                                               cgtggcgcggccAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATg                                                                                   >  1:1061218/4‑150 (MQ=255)
                                                                    cgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCAc                                                                              >  1:373362/1‑150 (MQ=255)
                                                                    cgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCAc                                                                              >  1:227034/1‑150 (MQ=255)
                                                                     gcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACg                                                                             <  1:48133/150‑1 (MQ=255)
                                                                       gggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAg                                                                           <  2:869395/150‑1 (MQ=255)
                                                                         gcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGcc                                                                         <  1:841824/150‑1 (MQ=255)
                                                                           aGTGGTAAGACGCGCGTACTGGTGCCTCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAg                                                                       <  2:1486301/150‑1 (MQ=255)
                                                                           aGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAg                                                                       <  2:994420/150‑1 (MQ=255)
                                                                             gggTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAggg                                                                     <  2:1556331/149‑1 (MQ=255)
                                                                              ggTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGc                                                                    <  2:770880/150‑1 (MQ=255)
                                                                               gTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCg                                                                   <  1:1611428/150‑1 (MQ=255)
                                                                                tAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCgg                                                                  <  1:1584833/150‑1 (MQ=255)
                                                                                  aGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTa                                                                <  1:1222399/150‑1 (MQ=255)
                                                                                   gACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTAt                                                               <  2:787047/150‑1 (MQ=255)
                                                                                     cgcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTAtgt                                                             >  1:1010496/1‑150 (MQ=255)
                                                                                      gcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTAtgtg                                                            <  2:494782/150‑1 (MQ=255)
                                                                                      gcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTAtgtg                                                            <  2:1407553/150‑1 (MQ=255)
                                                                                       cgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTgg                                                           >  1:1001558/1‑150 (MQ=255)
                                                                                          gTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTc                                                        >  2:1411365/1‑150 (MQ=255)
                                                                                           tACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCg                                                       <  2:1061218/150‑1 (MQ=255)
                                                                                           tACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCg                                                       >  2:164925/1‑150 (MQ=255)
                                                                                            aCTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCgg                                                      <  2:1059705/150‑1 (MQ=255)
                                                                                               ggTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCAc                                                   >  1:125632/1‑150 (MQ=255)
                                                                                               ggTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCAc                                                   >  1:755786/1‑150 (MQ=255)
                                                                                               ggTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCAc                                                   >  2:80618/1‑150 (MQ=255)
                                                                                                  gCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCtt                                                >  2:1212517/1‑150 (MQ=255)
                                                                                                     tCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCa                                             >  1:360316/1‑150 (MQ=255)
                                                                                                        tATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACgg                                          >  1:646206/1‑150 (MQ=255)
                                                                                                         aTCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACggg                                         <  1:908908/150‑1 (MQ=255)
                                                                                                         aTCGCATCGTTGATGAGCGTTGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGCTATGTGGGTCGGCACCTTCCACggg                                         >  2:1566410/1‑150 (MQ=255)
                                                                                                              cTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgc                                    <  1:1612533/150‑1 (MQ=255)
                                                                                                                   tgatgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCAcc                               <  1:639725/150‑1 (MQ=255)
                                                                                                                    gatgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCg                              >  1:1504342/1‑150 (MQ=255)
                                                                                                                      tgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGtt                            >  2:961217/1‑150 (MQ=255)
                                                                                                                      tgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGtt                            >  2:1277419/1‑150 (MQ=255)
                                                                                                                        aGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtg                          >  1:816293/1‑150 (MQ=255)
                                                                                                                            gggAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGGCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgctgc                      <  1:1428905/149‑1 (MQ=255)
                                                                                                                             ggAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcg                     <  1:1132167/150‑1 (MQ=255)
                                                                                                                              gAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcgt                    >  2:1195313/1‑150 (MQ=255)
                                                                                                                              gAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcgt                    <  1:46765/150‑1 (MQ=255)
                                                                                                                                   cTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCa               <  1:1212517/150‑1 (MQ=255)
                                                                                                                                   cTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCa               <  2:1342904/150‑1 (MQ=255)
                                                                                                                                    tcaccGCCATAATCGAGTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGACACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCAc              >  1:164034/5‑150 (MQ=255)
                                                                                                                                     gCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCAcc             <  1:164925/150‑1 (MQ=255)
                                                                                                                                      cTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCa            <  2:373362/150‑1 (MQ=255)
                                                                                                                                       tCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCat           >  1:419132/1‑150 (MQ=255)
                                                                                                                                       tCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCat           >  2:1362640/1‑150 (MQ=255)
                                                                                                                                       ccGCCATAATCGAGTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCat           >  2:398620/2‑150 (MQ=255)
                                                                                                                                           cATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATgg       >  2:398492/1‑150 (MQ=255)
                                                                                                                                           cATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATgg       >  2:300744/1‑150 (MQ=255)
                                                                                                                                               cTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGc   >  2:1024489/1‑150 (MQ=255)
                                                                                                                                                tCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGcc  <  2:755786/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCC  >  NC_000913/3998004‑3998297

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: