Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 713,205 T→C 6.0% R73R (CGT→CGC 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,2050TC 6.0% 636.9 / 18.4 186R73R (CGT→CGCseqAnegative modulator of initiation of replication
Reads supporting (aligned to +/- strand):  ref base T (85/89);  new base C (6/5);  total (91/95)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

GCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGTTTACTTTGCGG  >  NC_000913/713057‑713353
                                                                                                                                                    |                                                                                                                                                    
gcgcATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGt                                                                                                                                                      <  3:831824/149‑1 (MQ=255)
    aTCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAAc                                                                                                                                                  >  3:251937/1‑149 (MQ=255)
    aTCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAAc                                                                                                                                                  <  4:78300/149‑1 (MQ=255)
     tCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACt                                                                                                                                                 >  4:223795/1‑149 (MQ=255)
       cGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTc                                                                                                                                               <  2:1569813/149‑1 (MQ=255)
       cGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTc                                                                                                                                               >  3:1159402/1‑149 (MQ=255)
        gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCt                                                                                                                                              >  1:1474494/1‑149 (MQ=255)
        gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCt                                                                                                                                              <  2:189978/149‑1 (MQ=255)
        gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCt                                                                                                                                              >  1:40751/1‑149 (MQ=255)
        gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCt                                                                                                                                              >  1:496252/1‑149 (MQ=255)
        gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCACGCATTGCGCGAAc                                                                                                                                                  >  1:1614809/1‑145 (MQ=255)
        gACATTTTACGGCGTATGTTGAAATCTTCCGCAGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAATGTTCGCGCAATGCGTGAACTTCt                                                                                                                                              >  3:648804/1‑149 (MQ=255)
                 cGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGa                                                                                                                                     <  4:251937/149‑1 (MQ=255)
                     gTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGaa                                                                                                                                 <  2:496252/149‑1 (MQ=255)
                                ttttCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCa                                                                                                                      <  3:223795/149‑1 (MQ=255)
                                 tttCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCaa                                                                                                                     >  3:72837/1‑149 (MQ=255)
                                       ccgcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcga                                                                                                               <  2:40751/149‑1 (MQ=255)
                                         gcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcgagc                                                                                                             <  3:185158/149‑1 (MQ=255)
                                           aTCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGATTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCgg                                                                                                           >  2:1205354/1‑149 (MQ=255)
                                            tCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTGACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGGTTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGt                                                                                                          >  3:1584277/1‑149 (MQ=255)
                                              acaGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCTTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGc                                                                                                                                            <  3:83757/113‑1 (MQ=255)
                                              acaGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGc                                                                                                                                            >  4:83757/1‑113 (MQ=255)
                                               caGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCaa                                                                                                       >  3:1140103/1‑149 (MQ=255)
                                                   ctgctgCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGc                                                                                                   >  4:918359/1‑149 (MQ=255)
                                                            cGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTAtgctg                                                                                          >  3:1072842/1‑149 (MQ=255)
                                                            cGGTAACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTAtgctg                                                                                          <  4:1072842/149‑1 (MQ=255)
                                                             ggTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCt                                                                                                  >  2:1550462/1‑140 (MQ=255)
                                                              gTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATgctgct                                                                                        >  1:1266023/1‑149 (MQ=255)
                                                              gTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATgctgct                                                                                        <  2:389027/149‑1 (MQ=255)
                                                                   gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAAt                                                                                                                                >  3:1505895/1‑104 (MQ=255)
                                                                   gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAAt                                                                                                                                <  4:1505895/104‑1 (MQ=255)
                                                                   gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGt                                                                                                          >  2:1540828/1‑126 (MQ=255)
                                                                   gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                                   >  1:542456/1‑149 (MQ=255)
                                                                   gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                                   <  3:918359/149‑1 (MQ=255)
                                                                    aaaGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCa                                                                                                        <  1:1692823/127‑1 (MQ=255)
                                                                    aaaGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCa                                                                                                        >  2:1692823/1‑127 (MQ=255)
                                                                      agagGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcgag                                                                                                              >  3:213866/1‑119 (MQ=255)
                                                                      agagGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTAtgctg                                                                                          <  1:1194011/139‑1 (MQ=255)
                                                                      agagGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTa                                                                                <  4:1293015/149‑1 (MQ=255)
                                                                          gTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCt                                                                                                                                              >  1:602244/1‑83 (MQ=255)
                                                                           ttCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTc                                                                                                                                        >  1:182288/1‑88 (MQ=255)
                                                                           ttCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTc                                                                                                                                        <  2:182288/88‑1 (MQ=255)
                                                                               cgTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATAtt                                                                       <  1:1540931/149‑1 (MQ=255)
                                                                                gTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATtc                                                                      >  2:1673509/1‑149 (MQ=255)
                                                                                  tGCGTCACCTGCTATCGTCGAAGCGAAGCCTGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATtctc                                                                    <  4:1584277/149‑1 (MQ=255)
                                                                                  tGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATtctc                                                                    >  4:1254643/1‑149 (MQ=255)
                                                                                    cGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                                  >  1:705497/1‑149 (MQ=255)
                                                                                    cGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                                  <  4:1140103/149‑1 (MQ=255)
                                                                                      tCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCt                                                                                 <  2:1049384/132‑1 (MQ=255)
                                                                                      tCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCt                                                                                 >  1:1049384/1‑132 (MQ=255)
                                                                                      tCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATAtt                                                                       >  2:989240/1‑142 (MQ=255)
                                                                                      tCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATAtt                                                                       <  1:989240/142‑1 (MQ=255)
                                                                                         ccTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGc                                                             >  3:44112/1‑149 (MQ=255)
                                                                                          cTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGa                                                                >  1:1205736/1‑145 (MQ=255)
                                                                                          cTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGcc                                                            >  3:35113/1‑149 (MQ=255)
                                                                                          cTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGcc                                                            <  2:542456/149‑1 (MQ=255)
                                                                                           tGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                                  <  1:196363/142‑1 (MQ=255)
                                                                                           tGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGcc                                                            >  1:5144/1‑148 (MQ=255)
                                                                                           tGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGcc                                                            <  2:5144/148‑1 (MQ=255)
                                                                                                   tagaagcgaagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGttt                                                   >  4:1296990/3‑149 (MQ=255)
                                                                                                    cgaagcgaagCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTg                                                  <  3:1296990/149‑1 (MQ=255)
                                                                                                    cgaagcgaagCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTg                                                  <  4:44112/149‑1 (MQ=255)
                                                                                                     gaagcgaagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGc                                                 >  3:1895495/1‑149 (MQ=255)
                                                                                                       agcgaagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCg                                               <  3:1254643/149‑1 (MQ=255)
                                                                                                                 gTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTAt                                                                                              <  3:89352/92‑1 (MQ=255)
                                                                                                                 gTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTAt                                                                                              >  4:89352/1‑92 (MQ=255)
                                                                                                                     aaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                                   >  4:1407961/1‑99 (MQ=255)
                                                                                                                     aaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                                   <  3:1407961/99‑1 (MQ=255)
                                                                                                                      aaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGaa                                             <  4:477793/136‑1 (MQ=255)
                                                                                                                      aaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGaa                                             >  3:477793/1‑136 (MQ=255)
                                                                                                                       aCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGAAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGtt                               <  4:196961/149‑1 (MQ=255)
                                                                                                                         gATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAAt                                   <  1:1734561/143‑1 (MQ=255)
                                                                                                                           ttAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACg                                       >  4:1622443/1‑137 (MQ=255)
                                                                                                                           ttAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACg                                       <  3:1622443/137‑1 (MQ=255)
                                                                                                                           ttAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCAc                           <  4:35113/149‑1 (MQ=255)
                                                                                                                              aaGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACAGTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGt                        >  2:1291267/1‑149 (MQ=255)
                                                                                                                               aGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTacac                                                                             <  4:184257/95‑1 (MQ=255)
                                                                                                                                gACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTa                                                                                >  2:668141/1‑91 (MQ=255)
                                                                                                                                 aCAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGtt                                                    <  3:1848192/118‑1 (MQ=255)
                                                                                                                                   aaGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTacac                                                                             >  2:1242352/1‑91 (MQ=255)
                                                                                                                                   aaGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTg                                                  <  1:517081/118‑1 (MQ=255)
                                                                                                                                   aaGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTg                                                  >  2:517081/1‑118 (MQ=255)
                                                                                                                                   aaGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                                  <  2:137420/134‑1 (MQ=255)
                                                                                                                                    aagTTCACGCATTGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGATGTCTCCACTATATTCTCTTGACGCCCAGGCGTTTGCCGaa                                             >  4:72775/3‑122 (MQ=255)
                                                                                                                                      gTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACAcg                >  1:812809/1‑149 (MQ=255)
                                                                                                                                      gTTCACGCATTGCTCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTAt                                                                                              <  2:469795/71‑1 (MQ=255)
                                                                                                                                      gTTCACGCATTGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTAt                                                                                              >  1:469795/1‑71 (MQ=255)
                                                                                                                                      gTTCACGCATTGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTacac                 >  4:1875900/1‑148 (MQ=255)
                                                                                                                                       ttCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAGGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACAcgc               >  1:1358877/1‑149 (MQ=255)
                                                                                                                                       ttCACGCATTGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgtt                                                                                     >  2:1193113/1‑79 (MQ=255)
                                                                                                                                          gcgcAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGtt            >  4:115353/1‑149 (MQ=255)
                                                                                                                                          gcgcAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGt             >  1:1140708/1‑148 (MQ=255)
                                                                                                                                              aaTGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGTTTACt        >  2:287789/1‑149 (MQ=255)
                                                                                                                                              aaTGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGTTTACt        >  2:1781975/1‑149 (MQ=255)
                                                                                                                                                 gCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACAGTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGTTTACTTTg     <  1:1291267/149‑1 (MQ=255)
                                                                                                                                                   gTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                                   >  1:427215/1‑69 (MQ=255)
                                                                                                                                                   gTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                                   <  2:427215/69‑1 (MQ=255)
                                                                                                                                                    tGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTc                       >  4:379096/1‑128 (MQ=255)
                                                                                                                                                    tGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTc                       <  3:379096/128‑1 (MQ=255)
                                                                                                                                                    tGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGTTTACTTTGCgg  >  3:1816608/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGTTTACTTTGCGG  >  NC_000913/713057‑713353

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: