Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,613,622 C→A 9.5% R48S (CGT→AGT)  yjjV → putative DNase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,613,6220CA 9.5% 212.0 / 4.4 73R48S (CGT→AGT) yjjVputative DNase
Reads supporting (aligned to +/- strand):  ref base C (16/49);  new base A (6/2);  total (22/51)
Fisher's exact test for biased strand distribution p-value = 7.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAGGCGCTGGAAAGGCGTCCGGCGAAGGTGGT  >  NC_000913/4613483‑4613758
                                                                                                                                           |                                                                                                                                        
gATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCTACAGAGGCGGAAAATTTTGccc                                                                                                                                          >  2:357530/1‑140 (MQ=255)
gATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGccc                                                                                                                                          <  1:95249/140‑1 (MQ=255)
 aTTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGTAAATTTTGCCCg                                                                                                                                         >  2:153915/1‑140 (MQ=255)
      ccGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGACGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCgtgtgt                                                                                                                                    >  2:656174/1‑140 (MQ=255)
                   acCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATGTTGCCCGTGTGTTGGCATTAGCggt                                                                                                                       >  1:1123038/1‑139 (MQ=255)
                           ccATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGTGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATc                                                                                                               <  2:1306434/140‑1 (MQ=255)
                           ccATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATc                                                                                                               <  2:38744/140‑1 (MQ=255)
                             aTTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTTTTCCGGCAACAGAGGCGGAAAATTTTGCCAGTTTTTTTTAATTAGCGGAAAATTATCaa                                                                                                             <  1:901422/140‑1 (MQ=255)
                             aTTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCaa                                                                                                             <  2:906526/140‑1 (MQ=255)
                             aTTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCaa                                                                                                             <  1:1055447/140‑1 (MQ=255)
                             aTTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGACAGTGTGTTGTCATTAGCGGAAAATTATCaa                                                                                                             <  1:900642/140‑1 (MQ=255)
                             aTTTTGAGTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCaa                                                                                                             <  1:672613/140‑1 (MQ=255)
                                ttGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCg                                                                                                          <  1:40863/140‑1 (MQ=255)
                                  gATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCt                                                                                                        <  1:746254/140‑1 (MQ=255)
                                     ttCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTa                                                                                                     >  1:426426/1‑140 (MQ=255)
                                       cccGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATg                                                                                                   <  2:1145225/140‑1 (MQ=255)
                                       cccGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATg                                                                                                   <  2:1145201/140‑1 (MQ=255)
                                       accgccgTTTAGTGGCGATGTAGAGGCCAGCCTGCAACGCGCGTCTCTAGCGGGCTTAGGAAAGATCATTGTTCCGGCAACAGAGTCGGAAAATTTTTCCCGTTTGTTGGCATTAGCGGAAATTTATCATCATCTGTATg                                                                                                   <  2:1145741/139‑1 (MQ=255)
                                            cgTTTAGTGGCGATGAAGAGGCCAGACTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCa                                                                                              <  1:1508495/140‑1 (MQ=255)
                                             gTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCAt                                                                                             <  1:1245033/140‑1 (MQ=255)
                                             gTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCAt                                                                                             <  1:1496305/140‑1 (MQ=255)
                                             gTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCAt                                                                                             <  2:265421/140‑1 (MQ=255)
                                             gTTTAGTGGAGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGGTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCAt                                                                                             <  1:1496339/140‑1 (MQ=255)
                                                  gTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGc                                                                                        <  2:1035562/140‑1 (MQ=255)
                                                  gTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGc                                                                                        <  2:1035587/140‑1 (MQ=255)
                                                     gCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTg                                                                                     <  1:814264/140‑1 (MQ=255)
                                                      cGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCAGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGc                                                                                    >  2:1535389/1‑140 (MQ=255)
                                                      cGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCAGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGc                                                                                    >  2:1538043/1‑140 (MQ=255)
                                                      cGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCAGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGc                                                                                    >  2:82511/1‑140 (MQ=255)
                                                      cGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCAGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGc                                                                                    >  2:1577827/1‑140 (MQ=255)
                                                      cGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCAGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGc                                                                                    >  2:1576828/1‑140 (MQ=255)
                                                       gATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCa                                                                                   <  1:1467051/140‑1 (MQ=255)
                                                        aTGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCAt                                                                                  >  1:976668/1‑140 (MQ=255)
                                                                gCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTAt                                                                          >  1:188878/1‑140 (MQ=255)
                                                                      cTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGa                                                                    <  1:606378/140‑1 (MQ=255)
                                                                         cAACGCGCGGCTCATGCGGGAGTAGGCAAGCTCCTTGTTCCGGCAACAGAGGCGGAATATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAAACGCTGTATGCCGCATTGGGCTTGCCTCCTGGTATGTTGGaaaa                                                                 <  1:717527/140‑1 (MQ=255)
                                                                         cAACGCGCGGCACAAGCGGGCGTAGGCTAGATCATTGTTCCGGCAACAGAGGCGGAAAATATTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGaaaa                                                                 <  2:1570305/140‑1 (MQ=255)
                                                                         cAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGaaaa                                                                 <  1:789439/140‑1 (MQ=255)
                                                                         cAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGaaaa                                                                 <  1:77111/140‑1 (MQ=255)
                                                                              cgcgGCACACGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTATGCAAGTGTTTTGGCATTAGCGGAAAAATATCAACCGATGTATGCCGCAGTGGGCTTGCATCCGGGTATGTTGGAAAAACATa                                                            >  1:555490/1‑140 (MQ=255)
                                                                              cgcgGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATa                                                            >  1:555501/1‑140 (MQ=255)
                                                                              cgcgGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGGTGGAAAATCATa                                                            >  1:52272/1‑140 (MQ=255)
                                                                               gcgGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAg                                                           <  1:794256/140‑1 (MQ=255)
                                                                                cgGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGc                                                          <  1:395799/140‑1 (MQ=255)
                                                                                cgGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCAGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGc                                                          <  2:1049227/140‑1 (MQ=255)
                                                                                cgGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCAGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGc                                                          <  2:1049235/140‑1 (MQ=255)
                                                                                 ggCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCg                                                         <  2:1048298/140‑1 (MQ=255)
                                                                                 ggCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCg                                                         <  2:1047859/140‑1 (MQ=255)
                                                                                 cccacaAGCTGGCGTACCCAACAACATACTACCCCCAACAGACGCCGAAAAAATTGCCCGTGTGATGGCATTAGCGGAAAATTAACAACCGATGTATGCCGCATTGCGCTTGCATCCTCGTAAGATGGAAAAACAAAGCg                                                         <  2:1300228/138‑1 (MQ=255)
                                                                                  gCACCAGCGGGCGTAGGCTAGATCATTGGTCCGGCAACAGAGGCGGAAACTTGTGCCCGAGTGTTGGCATGAGCGGAAAAGTATCAACCGCTGTCTGCCGCATTGGGCTTGCATCCTGGAATGTTGGAAAAACATAGCGa                                                        <  2:7366/140‑1 (MQ=255)
                                                                                  gCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGa                                                        <  2:801726/140‑1 (MQ=255)
                                                                                   cacaAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGAt                                                       >  2:736610/1‑140 (MQ=255)
                                                                                     caAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGAtgt                                                     <  2:480644/140‑1 (MQ=255)
                                                                                     caAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGAtgt                                                     <  2:480083/140‑1 (MQ=255)
                                                                                     caAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGAtgt                                                     <  2:1596292/140‑1 (MQ=255)
                                                                                        gCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGATGGAAAAACATAGCGATGTGtc                                                  >  2:659524/1‑140 (MQ=255)
                                                                                              gTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGa                                            <  2:696676/140‑1 (MQ=255)
                                                                                                      gATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAATTTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTAc                                    >  2:1233248/1‑140 (MQ=255)
                                                                                                      gATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTAc                                    >  2:129863/1‑140 (MQ=255)
                                                                                                      gATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCGGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTAc                                    >  2:1548841/1‑140 (MQ=255)
                                                                                                      gATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCAGTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTAc                                    >  2:596501/1‑140 (MQ=255)
                                                                                                          aTTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTAcagca                                <  2:1613833/140‑1 (MQ=255)
                                                                                                          aTTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTAcagca                                <  2:1422008/140‑1 (MQ=255)
                                                                                                            tggtCCGGCAACAGAGGCGGTAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATGATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACCGCAgg                              <  2:30046/137‑1 (MQ=255)
                                                                                                            tGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAgg                              <  1:72753/140‑1 (MQ=255)
                                                                                                            tGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAgg                              <  2:1254601/140‑1 (MQ=255)
                                                                                                                   gCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAGGCGCTGGa                       >  2:1101868/1‑140 (MQ=255)
                                                                                                                          agGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAGGCGCTGGAAAGGCGt                <  2:196362/140‑1 (MQ=255)
                                                                                                                          agGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAGGCGCTGGAAAGGCGt                <  2:1630455/140‑1 (MQ=255)
                                                                                                                          agGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAGGCGCTGGAAAGGCGt                <  2:1420095/140‑1 (MQ=255)
                                                                                                                          agGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAGGCGCTGGAAAGGCGt                <  1:1690732/140‑1 (MQ=255)
                                                                                                                           ggCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAGGCGCTGGAAAGGCGTc               >  1:1591089/1‑140 (MQ=255)
                                                                                                                            gCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAGGCGCTGGAAAGGCGTcc              <  2:508172/140‑1 (MQ=255)
                                                                                                                                aaaaTTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAGGCGCTGGAAAGGCGTCCGGCg          <  2:834482/140‑1 (MQ=255)
                                                                                                                                        gCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAATCATAGCGATGTGTCTCTTGAGCAGCTACAGCAGGCGCTGGAAAGGCGTCCGGCGAAggtggt  <  2:1275315/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
GATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTATGCCGCATTGGGCTTGCATCCTGGTATGTTGGAAAAACATAGCGATGTGTCTCTTGAGCAGCTACAGCAGGCGCTGGAAAGGCGTCCGGCGAAGGTGGT  >  NC_000913/4613483‑4613758

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