Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 713,195 G→A 5.2% R70H (CGC→CAC)  seqA → negative modulator of initiation of replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913713,1950GA 5.2% 370.8 / 2.8 110R70H (CGC→CAC) seqAnegative modulator of initiation of replication
Reads supporting (aligned to +/- strand):  ref base G (48/56);  new base A (3/3);  total (51/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

TCGGCGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGT  >  NC_000913/713048‑713342
                                                                                                                                                   |                                                                                                                                                   
tCGGCGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTcgc                                                                                                                                                    <  1:1023360/149‑1 (MQ=255)
tCGGCGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTcgc                                                                                                                                                    <  3:725201/149‑1 (MQ=255)
    cGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACTAGGTTCGCGCaa                                                                                                                                                <  3:228865/149‑1 (MQ=255)
     gagaGCGCATCTGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACTATTAAAGACAAGGTTCGCGCAAt                                                                                                                                               >  4:583849/1‑149 (MQ=255)
     gagaGCGCATCCGACATTTTACGGCGTATGTTGGAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAAt                                                                                                                                               <  2:389786/149‑1 (MQ=255)
     gagaGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAAt                                                                                                                                               <  1:934155/149‑1 (MQ=255)
             aTCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAAc                                                                                                                                       <  3:642066/149‑1 (MQ=255)
              tCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACt                                                                                                                                      >  4:870021/1‑149 (MQ=255)
              tCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACt                                                                                                                                      >  4:729049/1‑149 (MQ=255)
               ccGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCACGCATTGCGCGAAcct                                                                                                                                     <  4:26925/149‑3 (MQ=255)
                cGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTc                                                                                                                                    <  1:1321/149‑1 (MQ=255)
                cGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTc                                                                                                                                    >  4:152336/1‑149 (MQ=255)
                 gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCt                                                                                                                                   >  3:325099/1‑149 (MQ=255)
                   cATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGc                                                                                                                                 <  3:1025738/149‑1 (MQ=255)
                    aTTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCt                                                                                                                                >  1:293177/1‑149 (MQ=255)
                          cGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGa                                                                                                                          <  1:358705/149‑1 (MQ=255)
                          cGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGa                                                                                                                          <  4:75994/149‑1 (MQ=255)
                           ggCGTATGTTGAAATTTTCCTCCGCATCACAGCCTGCTGCTCCGGTGATGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGGCAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGAt                                                                                                                         <  2:851860/149‑1 (MQ=255)
                           ggCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGAt                                                                                                                         <  1:493120/149‑1 (MQ=255)
                           ggCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGAt                                                                                                                         <  1:738183/149‑1 (MQ=255)
                             cGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGa                                                                                                                       >  2:608320/1‑149 (MQ=255)
                              gTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGaa                                                                                                                      >  1:136885/1‑149 (MQ=255)
                               tATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAAt                                                                                                                     <  2:892391/149‑1 (MQ=255)
                               tATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAAt                                                                                                                     <  3:974735/149‑1 (MQ=255)
                               tATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAAt                                                                                                                     >  4:253380/1‑149 (MQ=255)
                                   ttGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGc                                                                                                                 >  2:285524/1‑149 (MQ=255)
                                    tGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCa                                                                                                                <  4:470761/149‑1 (MQ=255)
                                    tGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCa                                                                                                                <  2:721822/149‑1 (MQ=255)
                                      aaaTTTTCGGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCATCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATGCGCaga                                                                                                              >  3:967145/1‑149 (MQ=255)
                                      aaaTTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCaga                                                                                                              <  2:322508/149‑1 (MQ=255)
                                         ttttCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCa                                                                                                           >  3:338850/1‑149 (MQ=255)
                                                ccgcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcga                                                                                                    <  4:429898/149‑1 (MQ=255)
                                                 cgcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcgag                                                                                                   <  1:14804/149‑1 (MQ=255)
                                                 cgcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcgag                                                                                                   >  2:788084/1‑149 (MQ=255)
                                                  gcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcgagc                                                                                                  <  2:110217/149‑1 (MQ=255)
                                                     tCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGt                                                                                               >  2:362400/1‑149 (MQ=255)
                                                     tCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGt                                                                                               <  3:549270/149‑1 (MQ=255)
                                                        caGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCaa                                                                                            <  1:622112/149‑1 (MQ=255)
                                                          gCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATc                                                                                          <  1:32940/149‑1 (MQ=255)
                                                            ctgctgCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGc                                                                                        <  4:399254/149‑1 (MQ=255)
                                                             tgctgcTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGTAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCt                                                                                       <  1:788084/149‑1 (MQ=255)
                                                             tgctgcTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCt                                                                                       <  1:1050366/149‑1 (MQ=255)
                                                             tgctgcTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCt                                                                                       <  2:601360/149‑1 (MQ=255)
                                                               ctgctCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGa                                                                                                                          <  3:729546/112‑1 (MQ=255)
                                                               ctgctCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGa                                                                                                                          >  4:729546/1‑112 (MQ=255)
                                                                tgctCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTa                                                                                    >  1:2125/1‑149 (MQ=255)
                                                                       gTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAAGGCTTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATgctgct                                                                             <  2:41932/149‑1 (MQ=255)
                                                                        tGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGctgctg                                                                            >  3:479966/1‑149 (MQ=255)
                                                                            gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                        >  2:485459/1‑149 (MQ=255)
                                                                            gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                        >  4:809916/1‑149 (MQ=255)
                                                                            gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                        >  2:1004389/1‑149 (MQ=255)
                                                                               agagGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATAAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTa                                                                     <  3:224489/149‑1 (MQ=255)
                                                                                   gTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACt                                                                 >  2:723105/1‑149 (MQ=255)
                                                                                    ttCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACta                                                                >  4:336098/1‑149 (MQ=255)
                                                                                       gcgTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTatat                                                             >  1:347599/1‑149 (MQ=255)
                                                                                       gcgTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTatat                                                             >  4:660343/1‑149 (MQ=255)
                                                                                         gTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATtc                                                           <  2:493084/149‑1 (MQ=255)
                                                                                            gCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTctct                                                        <  1:331746/149‑1 (MQ=255)
                                                                                             cGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                       <  4:359925/149‑1 (MQ=255)
                                                                                             cGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                       <  2:136885/149‑1 (MQ=255)
                                                                                                 aCCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGAt                                                                                                                         >  4:329626/1‑79 (MQ=255)
                                                                                                 aCCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGAt                                                                                                                         <  3:329626/79‑1 (MQ=255)
                                                                                                   cTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTacac                                                                  >  2:763477/1‑132 (MQ=255)
                                                                                                   cTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTacac                                                                  <  1:763477/132‑1 (MQ=255)
                                                                                                        aTCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCCCTTGACGCCCAGGc                                            <  2:957059/149‑1 (MQ=255)
                                                                                                           gtcgAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGtt                                         <  1:817867/149‑1 (MQ=255)
                                                                                                            tcgAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGttt                                        >  2:133569/1‑149 (MQ=255)
                                                                                                            tcgAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGttt                                        <  4:24451/149‑1 (MQ=255)
                                                                                                               aagcgaagaCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGcc                                     >  4:90887/1‑149 (MQ=255)
                                                                                                                 gcgaagcCGGTCAAAACGATTAAAGACGAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGa                                   >  3:703655/1‑149 (MQ=255)
                                                                                                                  cgaagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGaa                                  <  3:996118/149‑1 (MQ=255)
                                                                                                                  cgaagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGaa                                  <  3:983456/149‑1 (MQ=255)
                                                                                                                  cgaagcCGGTCAAAACGATTAAAGACAAGGTTCACGCAATGCGTGAACTTCTGCTTGCGGATGAATACGCAGAGCAAAAGCGAGCGGACAATCGCTTTATGCTGCTGTTGTCTACTCTATACTCTCTTGACGCACAGCCGTTTGCCGaa                                  <  2:780095/149‑1 (MQ=255)
                                                                                                                    aagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                       <  3:645101/126‑1 (MQ=255)
                                                                                                                    aagcCGGTCAAAACGATTAAAGACAAGGGTCACGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTAATGCTGCTGTTGTCTACACTATATTCTCtt                                                       >  4:645101/1‑126 (MQ=255)
                                                                                                                      gcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATAAGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCaa                              <  3:646629/149‑1 (MQ=255)
                                                                                                                        cGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACg                            <  3:152336/149‑1 (MQ=255)
                                                                                                                          gTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGa                          >  2:512346/1‑149 (MQ=255)
                                                                                                                             aaaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                       <  3:336249/149‑1 (MQ=255)
                                                                                                                             aaaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                       >  4:269555/1‑149 (MQ=255)
                                                                                                                             aaaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                       >  4:314342/1‑149 (MQ=255)
                                                                                                                             aaaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                       <  1:285524/149‑1 (MQ=255)
                                                                                                                               aaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                       <  1:404432/115‑1 (MQ=255)
                                                                                                                               aaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                       >  2:404432/1‑115 (MQ=255)
                                                                                                                               aaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGt                     >  2:570097/1‑149 (MQ=255)
                                                                                                                               aaCGATTAAAGACAAGGTTCACGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACACCCAGGCGTTTGCCGAAGCAACGGAATCGt                     >  1:849724/1‑149 (MQ=255)
                                                                                                                                  gATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGt                     >  4:621866/1‑146 (MQ=255)
                                                                                                                                   aTTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCa                 >  2:346560/1‑149 (MQ=255)
                                                                                                                                       aaGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                       >  1:353029/1‑139 (MQ=255)
                                                                                                                                       aaGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                       <  2:353029/139‑1 (MQ=255)
                                                                                                                                       aaGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGt             >  1:748189/1‑149 (MQ=255)
                                                                                                                                        aGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTc            <  2:698622/149‑1 (MQ=255)
                                                                                                                                        aGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTc            <  1:1004389/149‑1 (MQ=255)
                                                                                                                                           cAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTa         <  3:336098/149‑1 (MQ=255)
                                                                                                                                            gaagTTCACGCATTGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTac        >  2:1020302/4‑149 (MQ=255)
                                                                                                                                              ggTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTacac      >  1:378385/1‑149 (MQ=255)
                                                                                                                                               gTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACAcg     >  2:396519/1‑149 (MQ=255)
                                                                                                                                                ttCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACAcgc    >  3:816708/1‑149 (MQ=255)
                                                                                                                                                  cgcgCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGt  <  1:643867/149‑1 (MQ=255)
                                                                                                                                                  cgcgCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGt  <  1:542152/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGGCGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGT  >  NC_000913/713048‑713342

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