Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 713,195 G→A 9.8% 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 9.8% 266.3 / 16.1 92R70H (CGC→CAC) seqAnegative modulator of initiation of replication
Reads supporting (aligned to +/- strand):  ref base G (42/41);  new base A (4/5);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGTT  >  NC_000913/713053‑713343
                                                                                                                                              |                                                                                                                                                    
gagaGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAAt                                                                                                                                                <  1:178849/149‑1 (MQ=255)
gagaGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAAt                                                                                                                                                <  4:452259/149‑1 (MQ=255)
 agagCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATg                                                                                                                                               >  3:426153/1‑149 (MQ=255)
       cATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGaa                                                                                                                                         <  3:480409/149‑1 (MQ=255)
        aTCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAAc                                                                                                                                        <  4:115942/149‑1 (MQ=255)
        aTCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAAc                                                                                                                                        <  3:731056/149‑1 (MQ=255)
        aTCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCACGCATTGCGCGAAc                                                                                                                                        <  2:360899/149‑1 (MQ=255)
         tCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCACGCATTGCGCGAAcc                                                                                                                                       <  2:563367/149‑2 (MQ=255)
          ccGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCACGCATTGCGCGAAc                                                                                                                                        >  3:133539/1‑147 (MQ=255)
          ccGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCACGCATTGCGCGAAc                                                                                                                                        <  4:133539/147‑1 (MQ=255)
            gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCt                                                                                                                                    >  2:48471/1‑149 (MQ=255)
            gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCt                                                                                                                                    >  4:5761/1‑149 (MQ=255)
            gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCt                                                                                                                                    <  4:515430/149‑1 (MQ=255)
              cATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGc                                                                                                                                  <  4:503382/149‑1 (MQ=255)
              cATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCACGCATTGCGCGAAccttgtc                                                                                                                                  <  3:102366/149‑7 (MQ=255)
                  ttACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTc                                                                                                                              <  4:89445/149‑1 (MQ=255)
                  ttACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTc                                                                                                                              <  4:247242/149‑1 (MQ=255)
                  ttACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTc                                                                                                                              <  1:253882/149‑1 (MQ=255)
                          tATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAAt                                                                                                                      >  3:485106/1‑149 (MQ=255)
                             gTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACg                                                                                                                   >  2:324751/1‑149 (MQ=255)
                              ttGAAATTTTCCGCCGCATCACTGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTTAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGc                                                                                                                  >  1:567536/1‑149 (MQ=255)
                              ttGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGATAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGATCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTATCGGTTGAATACGc                                                                                                                  >  4:384887/1‑149 (MQ=255)
                               tGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCa                                                                                                                 <  2:52753/149‑1 (MQ=255)
                                    ttttCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCa                                                                                                            >  1:393993/1‑149 (MQ=255)
                                    ttttCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCa                                                                                                            >  4:594689/1‑149 (MQ=255)
                                       tCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCaaaa                                                                                                         <  2:670243/149‑1 (MQ=255)
                                       tCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCaaaa                                                                                                         <  2:286416/149‑1 (MQ=255)
                                          tccgcATCACAGCGTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGACAAAACGATTAAAGTCAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCACAagcg                                                                                                      >  4:491060/2‑149 (MQ=255)
                                           ccgcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCTAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcga                                                                                                     <  2:689436/149‑1 (MQ=255)
                                           ccgcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcga                                                                                                     >  3:432888/1‑149 (MQ=255)
                                           ccgcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcga                                                                                                     <  1:324751/149‑1 (MQ=255)
                                                  acaGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCa                                                                                              >  3:261129/1‑149 (MQ=255)
                                                           tgctCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTa                                                                                     <  4:311916/149‑1 (MQ=255)
                                                            gctCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTAt                                                                                    <  1:485914/149‑1 (MQ=255)
                                                                  gTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATgctgct                                                                              >  1:55962/1‑149 (MQ=255)
                                                                       gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                         >  2:332874/1‑149 (MQ=255)
                                                                       gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGGCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGTATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                         >  3:376288/1‑149 (MQ=255)
                                                                        aaaGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                        >  3:278109/1‑149 (MQ=255)
                                                                        aaaGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                        <  1:74949/149‑1 (MQ=255)
                                                                        aaaGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                        <  1:119173/149‑1 (MQ=255)
                                                                         aaGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCt                                                                       >  2:391311/1‑149 (MQ=255)
                                                                         aaGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCt                                                                       >  1:746398/1‑149 (MQ=255)
                                                                          agagGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTa                                                                      <  1:567108/149‑1 (MQ=255)
                                                                            agGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTaca                                                                    <  4:426153/149‑1 (MQ=255)
                                                                             ggTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTacac                                                                   <  3:610344/149‑1 (MQ=255)
                                                                                   cgTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTTAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATAtt                                                             <  2:567536/149‑1 (MQ=255)
                                                                                    gTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATtc                                                            >  4:311050/1‑149 (MQ=255)
                                                                                      tGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATtctc                                                          <  1:762684/149‑1 (MQ=255)
                                                                                       gCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTctct                                                         <  3:573948/149‑1 (MQ=255)
                                                                                        cGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                        <  2:693323/149‑1 (MQ=255)
                                                                                          tCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGa                                                      >  1:186252/1‑149 (MQ=255)
                                                                                           cACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGAc                                                     <  1:452163/149‑1 (MQ=255)
                                                                                            aCCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGTTGAATACGCAGAGCAAAagc                                                                                                       <  4:732320/98‑1 (MQ=255)
                                                                                            aCCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagc                                                                                                       >  3:732320/1‑98 (MQ=255)
                                                                                            aCCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACg                                                    >  1:105465/1‑149 (MQ=255)
                                                                                            aCCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACg                                                    >  2:83751/1‑149 (MQ=255)
                                                                                              cTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCt                                                                                        <  1:14868/111‑1 (MQ=255)
                                                                                              cTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCt                                                                                        >  2:14868/1‑111 (MQ=255)
                                                                                              cTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGcc                                                  >  4:662455/1‑149 (MQ=255)
                                                                                              cTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGcc                                                  >  4:61363/1‑149 (MQ=255)
                                                                                               tGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGccc                                                 >  1:76002/1‑149 (MQ=255)
                                                                                                  tATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGTACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACTCCCAgg                                              >  4:679240/1‑149 (MQ=255)
                                                                                                   aTCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGc                                             >  2:342424/1‑149 (MQ=255)
                                                                                                       tcgAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATtctc                                                          <  2:177896/132‑1 (MQ=255)
                                                                                                       tcgAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATtctc                                                          >  1:177896/1‑132 (MQ=255)
                                                                                                       tcgAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGttt                                         <  4:212278/149‑1 (MQ=255)
                                                                                                         gaagcgaagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGc                                       >  1:10321/1‑149 (MQ=255)
                                                                                                         gaagcgaagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGc                                       >  1:225934/1‑149 (MQ=255)
                                                                                                          aagcgaagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGcc                                      <  2:461384/149‑1 (MQ=255)
                                                                                                          aagcgaagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTGCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGcc                                      <  2:81245/149‑1 (MQ=255)
                                                                                                                        aaaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                        >  3:589688/1‑117 (MQ=255)
                                                                                                                        aaaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                        <  4:589688/117‑1 (MQ=255)
                                                                                                                        aaaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                        >  3:29380/1‑149 (MQ=255)
                                                                                                                        aaaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                        <  3:115463/149‑1 (MQ=255)
                                                                                                                         aaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCg                       >  2:282780/1‑149 (MQ=255)
                                                                                                                          aaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGttt                                         >  2:489681/1‑130 (MQ=255)
                                                                                                                          aaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGttt                                         <  1:489681/130‑1 (MQ=255)
                                                                                                                          aaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGt                      >  1:479227/1‑149 (MQ=255)
                                                                                                                          aaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGt                      >  2:18980/1‑149 (MQ=255)
                                                                                                                           aCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGtt                     <  2:235582/149‑1 (MQ=255)
                                                                                                                             gATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTACGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGc                   >  1:66320/1‑149 (MQ=255)
                                                                                                                                   aGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTc             <  2:396991/149‑1 (MQ=255)
                                                                                                                                   aGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTc             <  2:396876/149‑1 (MQ=255)
                                                                                                                                      cAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTa          >  3:86382/1‑149 (MQ=255)
                                                                                                                                      agaagTTCACGCATTGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGtt                     <  4:13494/134‑1 (MQ=255)
                                                                                                                                      agaagTTCACGCATTGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGtt                     >  3:13494/5‑138 (MQ=255)
                                                                                                                                       aaGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTac         >  4:104305/1‑149 (MQ=255)
                                                                                                                                       aaGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTac         >  1:325350/1‑149 (MQ=255)
                                                                                                                                        aGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTaca        <  4:376288/149‑1 (MQ=255)
                                                                                                                                          gTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACAcg      <  3:433321/149‑1 (MQ=255)
                                                                                                                                          gTTCACGCATTGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACAcg      >  2:339392/1‑149 (MQ=255)
                                                                                                                                          gTTCACGCATTGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACAcg      >  4:200815/1‑149 (MQ=255)
                                                                                                                                              gcgcAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGtt  <  1:18980/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGTT  >  NC_000913/713053‑713343

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