Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 713,195 G→A 5.9% R70H (CGC→CAC)  seqA → 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.9% 268.6 / 4.3 84R70H (CGC→CAC) seqASeqA, negative modulator of initiation of replication
Reads supporting (aligned to +/- strand):  ref base G (42/37);  new base A (2/3);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

TATCGGCGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGTTT  >  NC_000913/713046‑713344
                                                                                                                                                     |                                                                                                                                                     
tatCGGCGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTcg                                                                                                                                                       <  2:1403543/150‑1 (MQ=255)
tatCGGCGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTcg                                                                                                                                                       <  2:1443662/150‑1 (MQ=255)
  tCGGCGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTcacg                                                                                                                                                     >  2:1391469/1‑150 (MQ=255)
   cGGCGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCgcgc                                                                                                                                                    <  2:1702812/150‑1 (MQ=255)
    ggCGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCa                                                                                                                                                   >  1:1586151/1‑150 (MQ=255)
      cGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAAt                                                                                                                                                 >  2:1149131/1‑150 (MQ=255)
       gagaGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATg                                                                                                                                                <  2:577352/150‑1 (MQ=255)
         gagCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCg                                                                                                                                              <  1:249196/150‑1 (MQ=255)
               aTCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCACGCATTGCGCGAAcc                                                                                                                                        <  2:1228109/150‑2 (MQ=255)
                tCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACtt                                                                                                                                       >  2:256371/1‑150 (MQ=255)
                 ccGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTc                                                                                                                                      >  1:799850/1‑150 (MQ=255)
                  cGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCt                                                                                                                                     >  2:1234312/1‑150 (MQ=255)
                  cGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCt                                                                                                                                     <  2:580028/150‑1 (MQ=255)
                   gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTg                                                                                                                                    >  1:1479665/1‑150 (MQ=255)
                   gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTg                                                                                                                                    >  2:581679/1‑150 (MQ=255)
                   gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTg                                                                                                                                    <  2:1031748/150‑1 (MQ=255)
                   gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCACGCATTGCGCGAAccttgt                                                                                                                                    <  1:1391469/150‑6 (MQ=255)
                   gACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCACGCATTGCGCGAAccttgt                                                                                                                                    <  1:1793270/150‑6 (MQ=255)
                         ttACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCg                                                                                                                              <  2:491536/150‑1 (MQ=255)
                          tACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCgg                                                                                                                             <  2:1269507/150‑1 (MQ=255)
                             ggCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATg                                                                                                                          <  2:1177626/150‑1 (MQ=255)
                                    gTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGc                                                                                                                   >  2:1656617/1‑150 (MQ=255)
                                       gAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCaga                                                                                                                <  1:845295/150‑1 (MQ=255)
                                        aaaTTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCagag                                                                                                               <  2:881167/150‑1 (MQ=255)
                                         aaTTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGc                                                                                                              >  1:954735/1‑150 (MQ=255)
                                         aaTTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGc                                                                                                              <  2:1304927/150‑1 (MQ=255)
                                         aaTTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGc                                                                                                              <  2:855824/150‑1 (MQ=255)
                                              tCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAag                                                                                                         <  1:629335/150‑1 (MQ=255)
                                                 gccgcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcga                                                                                                      >  1:1987950/1‑150 (MQ=255)
                                                 gccgcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcga                                                                                                      >  1:1565237/1‑150 (MQ=255)
                                                   cgcATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAagcgagc                                                                                                    >  1:1903910/1‑150 (MQ=255)
                                                     cATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCgg                                                                                                  >  1:72316/1‑150 (MQ=255)
                                                     cATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCgg                                                                                                  >  2:772555/1‑150 (MQ=255)
                                                       tCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTc                                                                                                >  2:1588975/1‑150 (MQ=255)
                                                       tCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACCTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTc                                                                                                <  1:1365738/150‑1 (MQ=255)
                                                         acaGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCaa                                                                                              >  1:1480353/1‑150 (MQ=255)
                                                              ctgctgCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCt                                                                                         <  2:258147/150‑1 (MQ=255)
                                                              ctgctgCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGATTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCt                                                                                         >  2:986989/1‑150 (MQ=255)
                                                                 ctgctCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTa                                                                                      <  2:473282/150‑1 (MQ=255)
                                                                 ctgctCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTa                                                                                      <  2:1224510/150‑1 (MQ=255)
                                                                 ctgctCCGGTGACGAAAGACGGTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGCTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTa                                                                                      <  2:1026846/150‑1 (MQ=255)
                                                                            aCGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttg                                                                           <  1:256371/150‑1 (MQ=255)
                                                                             cGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                          >  2:1685137/1‑150 (MQ=255)
                                                                             cGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                          >  1:1376045/1‑150 (MQ=255)
                                                                              gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                         <  2:1586151/150‑1 (MQ=255)
                                                                              gAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                         <  1:1707886/150‑1 (MQ=255)
                                                                                 agagGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGTTGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTac                                                                      <  2:1698900/150‑1 (MQ=255)
                                                                                 agagGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTac                                                                      <  1:1043435/150‑1 (MQ=255)
                                                                                 agagGTTCGCGTTACGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTac                                                                      >  2:1764775/1‑150 (MQ=255)
                                                                                   agGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTacac                                                                    <  1:1149131/150‑1 (MQ=255)
                                                                                       tCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTAtgc                                                                                   <  1:628507/131‑1 (MQ=255)
                                                                                       tCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTAtgc                                                                                   >  2:628507/1‑131 (MQ=255)
                                                                                       tCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACtata                                                                >  1:305259/1‑150 (MQ=255)
                                                                                       tCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACtata                                                                <  2:1565237/150‑1 (MQ=255)
                                                                                             tacGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTctct                                                          <  1:1764775/148‑1 (MQ=255)
                                                                                             tGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTctct                                                          >  2:799544/1‑150 (MQ=255)
                                                                                                 tCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGAc                                                      <  1:1685137/150‑1 (MQ=255)
                                                                                                  cACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACg                                                     >  2:1349784/1‑150 (MQ=255)
                                                                                                    ccTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGcc                                                   <  2:799850/150‑1 (MQ=255)
                                                                                                      tGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAGTACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCa                                                 >  1:1292649/1‑150 (MQ=255)
                                                                                                        cTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAgg                                               <  1:771225/150‑1 (MQ=255)
                                                                                                        cTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAgg                                               <  2:1376045/150‑1 (MQ=255)
                                                                                                           tcgtcgAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGt                                            >  2:1162013/1‑150 (MQ=255)
                                                                                                             gtcgAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGtt                                           >  1:962960/1‑149 (MQ=255)
                                                                                                             gtcgAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGtt                                           <  2:962960/149‑1 (MQ=255)
                                                                                                              tcgAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTg                                         <  2:954735/150‑1 (MQ=255)
                                                                                                                 aagcgaagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCg                                      >  2:383815/1‑150 (MQ=255)
                                                                                                                  agcgaagcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGa                                     >  1:395334/1‑150 (MQ=255)
                                                                                                                       agcCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCaa                                >  1:1128426/1‑150 (MQ=255)
                                                                                                                         ccGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACg                              <  2:1903910/150‑1 (MQ=255)
                                                                                                                            gTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGaa                           >  1:1848886/1‑150 (MQ=255)
                                                                                                                              cAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                         >  2:386805/1‑150 (MQ=255)
                                                                                                                              cAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                         >  2:430767/1‑150 (MQ=255)
                                                                                                                               aaaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCg                        >  1:1864848/1‑150 (MQ=255)
                                                                                                                                aaaCGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGt                       >  1:1339023/1‑150 (MQ=255)
                                                                                                                                     aTTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCAc                  <  2:952140/150‑1 (MQ=255)
                                                                                                                                     aTTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCAc                  >  1:1830973/1‑150 (MQ=255)
                                                                                                                                           gACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGt            >  1:882381/1‑150 (MQ=255)
                                                                                                                                            cagaagTTCACGCATTGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTa           >  2:22578/6‑150 (MQ=255)
                                                                                                                                            cagaagTTCACGCATTGCGCGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTa           >  1:620016/6‑150 (MQ=255)
                                                                                                                                            aCAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTa           >  1:808132/1‑150 (MQ=255)
                                                                                                                                             cAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCg             <  2:1533483/147‑1 (MQ=255)
                                                                                                                                             cAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCg             >  1:1533483/1‑147 (MQ=255)
                                                                                                                                             cAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTac          <  1:1162013/150‑1 (MQ=255)
                                                                                                                                             cAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTac          >  2:430715/1‑150 (MQ=255)
                                                                                                                                              aaGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTaca         <  2:305259/150‑1 (MQ=255)
                                                                                                                                                ggTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACAcg       >  1:967900/1‑150 (MQ=255)
                                                                                                                                                ggTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACAcg       >  2:1421944/1‑150 (MQ=255)
                                                                                                                                                     gcgcAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGttt  >  1:794608/1‑150 (MQ=255)
                                                                                                                                                     gcgcAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGttt  >  1:1888878/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATCGGCGAGAGCGCATCCGACATTTTACGGCGTATGTTGAAATTTTCCGCCGCATCACAGCCTGCTGCTCCGGTGACGAAAGAGGTTCGCGTTGCGTCACCTGCTATCGTCGAAGCGAAGCCGGTCAAAACGATTAAAGACAAGGTTCGCGCAATGCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACACGCGTTT  >  NC_000913/713046‑713344

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: