Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,520,486 G→A 5.0% T73I (ACT→ATT)  hofQ ← DNA utilization protein HofQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,520,4860GA 5.0% 322.4 / 4.1 99T73I (ACT→ATT) hofQDNA utilization protein HofQ
Reads supporting (aligned to +/- strand):  ref base G (49/45);  new base A (3/2);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

CTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACT  >  NC_000913/3520342‑3520635
                                                                                                                                                |                                                                                                                                                     
cTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                  <  1:1232408/150‑1 (MQ=255)
cTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                  >  1:167390/1‑150 (MQ=255)
 tGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                  >  2:1152437/1‑149 (MQ=255)
 tGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                                                                                                                                  <  1:1152437/149‑1 (MQ=255)
  gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAg                                                                                                                                                >  1:855365/1‑150 (MQ=255)
  gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGAGGGCTATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAg                                                                                                                                                >  2:1026156/1‑150 (MQ=255)
   cGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGt                                                                                                                                               >  1:1789677/1‑150 (MQ=255)
    gATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTg                                                                                                                                              >  2:838860/1‑150 (MQ=255)
     aTTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAAACTGTAGTga                                                                                                                                             <  1:828885/150‑2 (MQ=255)
      ttttCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGcc                                                                                                                                            <  2:167390/150‑1 (MQ=255)
       tttCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAAACTGTAGTgaaa                                                                                                                                           <  1:467620/150‑4 (MQ=255)
             gCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTcc                                                                                                                                     >  2:1354496/1‑150 (MQ=255)
             gCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTcc                                                                                                                                     >  2:1105302/1‑150 (MQ=255)
               ggCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCGGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAg                                                                                                                                   >  1:1116413/1‑150 (MQ=255)
               ggCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAg                                                                                                                                   >  1:254643/1‑150 (MQ=255)
                gCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAgg                                                                                                                                  <  2:1531150/150‑1 (MQ=255)
                 cAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAggg                                                                                                                                 >  2:1967158/1‑150 (MQ=255)
                 cAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAggg                                                                                                                                 >  1:1539587/1‑150 (MQ=255)
                 cAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAggg                                                                                                                                 <  2:1157054/150‑1 (MQ=255)
                 cAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAggg                                                                                                                                 <  2:1539587/150‑1 (MQ=255)
                    aTTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAAc                                                                                                                              <  1:471127/150‑1 (MQ=255)
                       tGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGTACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATc                                                                                                                           <  1:636738/150‑1 (MQ=255)
                       tGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATc                                                                                                                           <  2:1477704/150‑1 (MQ=255)
                         ccTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTg                                                                                                                         >  2:615765/1‑150 (MQ=255)
                           tGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                                                                                                                       <  2:326640/150‑1 (MQ=255)
                               cccGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGAt                                                                                                                   >  2:1314634/1‑150 (MQ=255)
                                ccGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATg                                                                                                                  <  2:734258/150‑1 (MQ=255)
                                 cgcgCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGt                                                                                                                 <  2:1789677/150‑1 (MQ=255)
                                 cgcgCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGt                                                                                                                 >  2:979813/1‑150 (MQ=255)
                                    gcCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAAc                                                                                                              <  2:1862084/150‑1 (MQ=255)
                                      cTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGa                                                                                                            >  1:232853/1‑150 (MQ=255)
                                        gCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATc                                                                                                                           >  2:531911/1‑133 (MQ=255)
                                        gCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATc                                                                                                                           <  1:531911/133‑1 (MQ=255)
                                        gCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTGCCAGGGAACATCTGTTAGATGTAACGacc                                                                                                          >  2:67474/1‑149 (MQ=255)
                                             gCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTa                                                                                                     >  1:1457313/1‑150 (MQ=255)
                                              cctcctGGCTGGCGATATTGTTATTCTGCCAGGGAATGGAATGCACTGATAGAATGTTGCCTTCCTGCCGCGTTATCAGTGCTGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCATGGAACATCTGTTATATGTAACGACACCGTAc                                                                                                    <  2:196566/150‑1 (MQ=255)
                                                   tGGCGGGCGATATTGTTTTTGTGCCAGGCAATGTAATGCACTGAGTGATTGTTTCCTTCCTGCCGCGTTATCAGTCGGGCGCTTTTAACTAGAGTTTGTAGTGGCTGCTTCCATGGAACATCTTTTTGATGTAACGACACAGTACCGCTg                                                                                               <  2:1840267/150‑1 (MQ=255)
                                                   tGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTg                                                                                               <  1:110169/150‑1 (MQ=255)
                                                     gcgggcgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGAc                                                                                             <  1:1967158/150‑1 (MQ=255)
                                                         gcgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCt                                                                                         >  2:1110503/1‑150 (MQ=255)
                                                          cgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTg                                                                                        >  1:83409/1‑150 (MQ=255)
                                                            atatTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGc                                                                                      <  2:254643/150‑1 (MQ=255)
                                                            atatTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGc                                                                                      <  2:5839/150‑1 (MQ=255)
                                                               ttGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGac                                                                                   <  1:1845489/150‑1 (MQ=255)
                                                                  ttATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACg                                                                                >  1:24539/1‑150 (MQ=255)
                                                                    aTTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGAc                                                                              <  2:1359369/150‑1 (MQ=255)
                                                                       cTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCAGGCGACACGACCAg                                                                           >  1:554699/1‑150 (MQ=255)
                                                                        tGCCAGGCAATGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGt                                                                         <  1:400950/150‑1 (MQ=255)
                                                                         gCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGt                                                                         >  2:438388/1‑150 (MQ=255)
                                                                          ccAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGtt                                                                        <  2:621880/150‑1 (MQ=255)
                                                                           cAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTc                                                                       >  1:1887060/1‑150 (MQ=255)
                                                                            aGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCa                                                                      <  2:442131/150‑1 (MQ=255)
                                                                              gCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGGGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                                                    <  2:1060853/150‑1 (MQ=255)
                                                                              gCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                                                    <  1:1773254/150‑1 (MQ=255)
                                                                               cAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACt                                                                   >  1:88114/1‑150 (MQ=255)
                                                                               cAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACt                                                                   <  2:855365/150‑1 (MQ=255)
                                                                                 aTGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtc                                                                 >  1:1588002/1‑150 (MQ=255)
                                                                                    gAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTcc                                                              <  2:1494804/150‑1 (MQ=255)
                                                                                    gAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTcc                                                              <  2:1730015/150‑1 (MQ=255)
                                                                                    gAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTcc                                                              <  1:286180/150‑1 (MQ=255)
                                                                                       tGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGt                                                           <  1:438388/150‑1 (MQ=255)
                                                                                       tGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGt                                                           >  1:624722/1‑150 (MQ=255)
                                                                                         cACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTc                                                         >  2:1475272/1‑150 (MQ=255)
                                                                                         cACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTc                                                         <  1:615765/150‑1 (MQ=255)
                                                                                          aCTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCAGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTca                                                        >  1:286330/1‑150 (MQ=255)
                                                                                           cTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcag                                                       >  1:1179925/1‑150 (MQ=255)
                                                                                           cTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcag                                                       >  2:911088/1‑150 (MQ=255)
                                                                                             gagagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAg                                                                                                                     >  1:1746010/1‑86 (MQ=255)
                                                                                             gagagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAg                                                                                                                     <  2:1746010/86‑1 (MQ=255)
                                                                                             gagagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagcc                                                     <  1:1445371/150‑1 (MQ=255)
                                                                                                   aTGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGGCAGCGcc                                               >  1:465757/1‑150 (MQ=255)
                                                                                                   aTGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGcc                                               <  2:1588002/150‑1 (MQ=255)
                                                                                                    tGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCt                                              >  2:660902/1‑150 (MQ=255)
                                                                                                       ggCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGGGCCTGCTTCCAGGGAACAGCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCa                                           <  2:1444809/149‑1 (MQ=255)
                                                                                                         ccTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcag                                                       <  2:1543078/136‑1 (MQ=255)
                                                                                                         ccTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcag                                                       >  1:1543078/1‑136 (MQ=255)
                                                                                                             ccTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCt                                     <  2:1887060/150‑1 (MQ=255)
                                                                                                             ccTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCt                                     >  2:615574/1‑150 (MQ=255)
                                                                                                              cTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTg                                    <  1:660902/150‑1 (MQ=255)
                                                                                                               tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGa                                   <  1:979813/150‑1 (MQ=255)
                                                                                                                gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg                                  >  2:215886/1‑150 (MQ=255)
                                                                                                                gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg                                  >  1:1467900/1‑150 (MQ=255)
                                                                                                                 ccGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGc                                 >  1:260530/1‑150 (MQ=255)
                                                                                                                  cgcgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCt                                <  2:1467900/150‑1 (MQ=255)
                                                                                                                              ccGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCa                    >  2:720074/1‑150 (MQ=255)
                                                                                                                                ggCGCTTTTCACTACAGTTTGTAGTGACTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATc                  >  2:1649979/1‑150 (MQ=255)
                                                                                                                                  ggcTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATcca                <  1:1026156/149‑1 (MQ=255)
                                                                                                                                    gcaggcACTACAAACTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCAGAACGTCATccacc              >  1:1577985/6‑150 (MQ=255)
                                                                                                                                     caggcACTACAAACTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATccacca             >  2:153182/5‑150 (MQ=255)
                                                                                                                                      aggcACTACAAACTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCAt            >  1:843689/4‑150 (MQ=255)
                                                                                                                                        tCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCa                                           >  2:122360/1‑117 (MQ=255)
                                                                                                                                        tCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCa                                           <  1:122360/117‑1 (MQ=255)
                                                                                                                                         cACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCCCCATCAg         >  1:899947/1‑150 (MQ=255)
                                                                                                                                          aCTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCCTCCACCATCAGc        >  2:1143523/1‑150 (MQ=255)
                                                                                                                                           cTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCg       >  2:1337920/1‑150 (MQ=255)
                                                                                                                                           cTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCg       >  1:1751847/1‑150 (MQ=255)
                                                                                                                                            tACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGt      >  1:363604/1‑150 (MQ=255)
                                                                                                                                              cAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCa    <  1:215886/150‑1 (MQ=255)
                                                                                                                                               aGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCAc   >  2:1790929/1‑150 (MQ=255)
                                                                                                                                                gTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACt  >  1:860763/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACT  >  NC_000913/3520342‑3520635

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: