Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,325,285 G→C 6.4% intergenic (‑120/+13) phnC ← / ← yjdN phosphonate ABC transporter ATP binding subunit/conserved protein YjdN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,325,2850GC 6.4% 426.8 / 15.2 125intergenic (‑120/+13)phnC/yjdNphosphonate ABC transporter ATP binding subunit/conserved protein YjdN
Reads supporting (aligned to +/- strand):  ref base G (70/47);  new base C (4/4);  total (74/51)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 130 total aligned reads displayed.

GAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGATTGTCAA  >  NC_000913/4325136‑4325434
                                                                                                                                                     |                                                                                                                                                     
gAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAgggcggg                                                                                                                                                       <  1:1151384/150‑1 (MQ=255)
gAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAgggcggg                                                                                                                                                       <  1:1362601/150‑1 (MQ=255)
  gCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGcc                                                                                                                                                     <  2:1062865/150‑1 (MQ=255)
      ctcGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTccc                                                                                                                                                 <  2:65291/150‑1 (MQ=255)
      ctcGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTccc                                                                                                                                                 >  1:679914/1‑150 (MQ=255)
        cGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTTGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGCCCGCCCTCCcgc                                                                                                                                               <  2:476278/150‑2 (MQ=255)
        cGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCcgg                                                                                                                                               >  1:394280/1‑150 (MQ=255)
        cGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCcgg                                                                                                                                               >  2:1407824/1‑150 (MQ=255)
         gACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCcggc                                                                                                                                              >  1:716461/1‑150 (MQ=255)
           cacGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGCCCGCCCTCCcgcact                                                                                                                                            <  1:192594/150‑5 (MQ=255)
            acGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCggg                                                                                                                                           <  2:1342859/150‑1 (MQ=255)
            acGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCggg                                                                                                                                           <  1:1335293/150‑1 (MQ=255)
            acGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCggg                                                                                                                                           <  1:1098022/150‑1 (MQ=255)
            acGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTCGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCggg                                                                                                                                           <  1:233060/150‑1 (MQ=255)
             cGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGt                                                                                                                                          <  2:251426/150‑1 (MQ=255)
                 tAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTAtt                                                                                                                                      >  2:494265/1‑150 (MQ=255)
                   aTCGTCTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGc                                                                                                                                    <  1:597561/150‑1 (MQ=255)
                        ttGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttgg                                                                                                                               <  2:724561/150‑1 (MQ=255)
                        ttGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttgg                                                                                                                               >  2:1328827/1‑150 (MQ=255)
                        ttGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttgg                                                                                                                               <  2:20100/150‑1 (MQ=255)
                          gCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggtt                                                                                                                             >  1:1176858/1‑150 (MQ=255)
                          gCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGCGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggtt                                                                                                                             >  1:371538/1‑150 (MQ=255)
                           cATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggttg                                                                                                                            >  2:225745/1‑150 (MQ=255)
                           cATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggttg                                                                                                                            >  2:996825/1‑150 (MQ=255)
                           cATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggttg                                                                                                                            >  1:332633/1‑150 (MQ=255)
                            atgatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgt                                                                                                                           <  1:1172562/150‑1 (MQ=255)
                             tgatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGCCCGCCCTCCcgca                                                                                                                                              <  2:150357/130‑3 (MQ=255)
                              gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgttg                                                                                                                         >  1:860181/1‑150 (MQ=255)
                                tgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgttgtt                                                                                                                       <  1:441039/150‑1 (MQ=255)
                                tgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgttgtt                                                                                                                       >  2:1155800/1‑150 (MQ=255)
                                  cAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTTCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTg                                                                                                                     >  2:1073066/1‑150 (MQ=255)
                                  cAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTg                                                                                                                     <  2:1096215/150‑1 (MQ=255)
                                   aGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTga                                                                                                                    >  2:834926/1‑150 (MQ=255)
                                       tCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTgacgac                                                                                                                <  1:284678/150‑1 (MQ=255)
                                       tCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTgacgac                                                                                                                <  1:484528/150‑1 (MQ=255)
                                        cccTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACa                                                                                                               <  2:274793/150‑1 (MQ=255)
                                             aaaaaGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAat                                                                                                          >  2:599031/1‑150 (MQ=255)
                                             aaaaaGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAat                                                                                                          >  2:1275840/1‑150 (MQ=255)
                                             aaaaaGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAat                                                                                                          >  1:615671/1‑150 (MQ=255)
                                                 tgtgtCTCCATCGTGGAGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACATATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAatgatc                                                                                                      <  2:1052970/149‑1 (MQ=255)
                                                 aGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAatcatc                                                                                                      <  1:1296706/150‑1 (MQ=255)
                                                       tCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGt                                                                                                >  1:194292/1‑150 (MQ=255)
                                                           tCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGc                                                                                            <  2:716461/150‑1 (MQ=255)
                                                             gTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGAAGACATTAATCATCCACGGTACGCCa                                                                                          >  1:1266913/1‑150 (MQ=255)
                                                                      tattGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCg                                                                                 <  2:363686/150‑1 (MQ=255)
                                                                        ttGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGt                                                                               <  1:495593/150‑1 (MQ=255)
                                                                         tgtgACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTg                                                                              <  1:1155800/150‑1 (MQ=255)
                                                                           tgACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGaa                                                                            >  2:1334673/1‑149 (MQ=255)
                                                                           tgACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGAc                                                                            >  1:1064980/1‑150 (MQ=255)
                                                                             aCGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACtt                                                                          >  2:384060/1‑150 (MQ=255)
                                                                              cGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACttt                                                                         >  1:159979/1‑150 (MQ=255)
                                                                               gATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTg                                                                        <  2:802313/150‑1 (MQ=255)
                                                                               gATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTg                                                                        >  1:861506/1‑150 (MQ=255)
                                                                               gATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTg                                                                        >  2:904300/1‑150 (MQ=255)
                                                                                 tttGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGcc                                                                      <  2:679914/150‑1 (MQ=255)
                                                                                    ggTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCaaa                                                                   >  1:16242/1‑150 (MQ=255)
                                                                                    ggTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCaaa                                                                   >  2:757831/1‑150 (MQ=255)
                                                                                     gTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAg                                                                  >  2:270453/1‑150 (MQ=255)
                                                                                     gTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAg                                                                  >  1:812977/1‑150 (MQ=255)
                                                                                      ttaattaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGc                                                                 >  2:981255/1‑150 (MQ=255)
                                                                                       taattaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGcc                                                                <  2:860181/150‑1 (MQ=255)
                                                                                        aattaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCa                                                               >  1:1266275/1‑150 (MQ=255)
                                                                                        aattaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCa                                                               <  2:1427964/150‑1 (MQ=255)
                                                                                               aGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGccc                                                        >  2:1179274/1‑150 (MQ=255)
                                                                                               aGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGccc                                                        >  1:1052668/1‑150 (MQ=255)
                                                                                                  gACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGa                                                     >  2:185582/1‑150 (MQ=255)
                                                                                                    cTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaaa                                                   <  2:783673/150‑1 (MQ=255)
                                                                                                    cTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaaa                                                   >  1:935788/1‑150 (MQ=255)
                                                                                                      aaCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGt                                                 >  1:765695/1‑150 (MQ=255)
                                                                                                      aaCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACACCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGt                                                 >  1:1194913/1‑150 (MQ=255)
                                                                                                        cAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGttt                                               >  1:1028700/1‑150 (MQ=255)
                                                                                                           aTGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCt                                            <  1:165213/150‑1 (MQ=255)
                                                                                                            tGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTg                                           >  2:1238474/1‑150 (MQ=255)
                                                                                                              aaGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCTGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGcc                                         <  2:527302/150‑1 (MQ=255)
                                                                                                              aaGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGcc                                         >  2:1397434/1‑150 (MQ=255)
                                                                                                                     aaCGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCAt                                  >  1:495422/1‑150 (MQ=255)
                                                                                                                     aaCGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCAt                                  >  2:1254545/1‑150 (MQ=255)
                                                                                                                       cGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCAttt                                >  1:499500/1‑150 (MQ=255)
                                                                                                                       cGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCAttt                                >  1:1402095/1‑150 (MQ=255)
                                                                                                                           aaTTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGAt                            <  1:834926/150‑1 (MQ=255)
                                                                                                                             ttCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGAttt                          <  1:788580/150‑1 (MQ=255)
                                                                                                                               cGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGAttttt                        >  2:109929/1‑150 (MQ=255)
                                                                                                                                gATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTc                       <  1:1238474/150‑1 (MQ=255)
                                                                                                                                gATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTc                       >  2:1024858/1‑150 (MQ=255)
                                                                                                                                gATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTc                       <  1:694954/150‑1 (MQ=255)
                                                                                                                                  tGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCg                     <  2:1429748/150‑1 (MQ=255)
                                                                                                                                  tGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCg                     >  2:600384/1‑150 (MQ=255)
                                                                                                                                  tGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCg                     >  1:735165/1‑150 (MQ=255)
                                                                                                                                    accgTGCGGGAGGGCCGGCCTCCCGGCGGGTAATTGCGTTGGTGGTGGTTTGACGACATAAATCATCCACGGTACGCCAAATTAAGCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGtt                   <  2:752148/147‑1 (MQ=255)
                                                                                                                                     cAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGttg                  >  1:414325/1‑149 (MQ=255)
                                                                                                                                      ccgccGGGAGGCCCGCCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAg                 >  1:1157143/5‑150 (MQ=255)
                                                                                                                                       gTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGcc                                         <  2:12106/125‑1 (MQ=255)
                                                                                                                                       gTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGATGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGc                >  2:326612/1‑150 (MQ=255)
                                                                                                                                       cgccGGGAGGCCCGCCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGc                <  2:1157143/147‑1 (MQ=255)
                                                                                                                                       cgccGGGAGGCCCGCCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGc                >  1:325392/4‑150 (MQ=255)
                                                                                                                                          cGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGc             >  2:737453/1‑150 (MQ=255)
                                                                                                                                            ggAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCaa           >  2:1175819/1‑150 (MQ=255)
                                                                                                                                            ggAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCaa           >  1:583670/1‑150 (MQ=255)
                                                                                                                                               gggcgggcCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGAt        <  2:382652/150‑1 (MQ=255)
                                                                                                                                                     gcCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGATTGTCaa  <  2:499500/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGATTGTCAA  >  NC_000913/4325136‑4325434

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