Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,325,281 G→C 6.4% intergenic (‑116/+17) 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,2810GC 6.4% 427.7 / 15.1 125intergenic (‑116/+17)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 = 8.07e-01

Only 100 of 132 total aligned reads displayed.

TGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGATT  >  NC_000913/4325132‑4325429
                                                                                                                                                     |                                                                                                                                                    
tGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGagggaggg                                                                                                                                                      >  1:787468/1‑150 (MQ=255)
tGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAggg                                                                                                                                                      <  1:286/150‑1 (MQ=255)
   cGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAgggcgg                                                                                                                                                   >  1:1052970/1‑150 (MQ=255)
    gAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAgggcggg                                                                                                                                                  <  1:1151384/150‑1 (MQ=255)
      gCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGcc                                                                                                                                                <  2:1062865/150‑1 (MQ=255)
          ctcGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTccc                                                                                                                                            >  1:679914/1‑150 (MQ=255)
          ctcGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTccc                                                                                                                                            <  2:65291/150‑1 (MQ=255)
            cGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCcgg                                                                                                                                          >  2:1407824/1‑150 (MQ=255)
             gACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCcggc                                                                                                                                         >  1:716461/1‑150 (MQ=255)
               cacGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGCCCGCCCTCCcgcact                                                                                                                                       <  1:192594/150‑5 (MQ=255)
                acGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCggg                                                                                                                                      <  1:1335293/150‑1 (MQ=255)
                acGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTCGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCggg                                                                                                                                      <  1:233060/150‑1 (MQ=255)
                 cGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGt                                                                                                                                     <  2:251426/150‑1 (MQ=255)
                  ggATAATCGTCTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGtt                                                                                                                                    >  2:597561/1‑150 (MQ=255)
                     tAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTAtt                                                                                                                                 >  2:494265/1‑150 (MQ=255)
                      aaTCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTg                                                                                                                                <  2:731105/150‑1 (MQ=255)
                       aTCGTCTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGc                                                                                                                               <  1:597561/150‑1 (MQ=255)
                          gTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgtt                                                                                                                            >  1:173974/1‑150 (MQ=255)
                          gTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgtt                                                                                                                            <  1:1168476/150‑1 (MQ=255)
                           tttGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCTTCCCGGCGGGTTATTGCgttg                                                                                                                           <  2:1338815/150‑1 (MQ=255)
                            ttGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttgg                                                                                                                          <  2:724561/150‑1 (MQ=255)
                            ttGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttgg                                                                                                                          <  2:20100/150‑1 (MQ=255)
                            ttGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttgg                                                                                                                          >  2:1328827/1‑150 (MQ=255)
                              gCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggtt                                                                                                                        >  1:1176858/1‑150 (MQ=255)
                              gCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGCGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggtt                                                                                                                        >  1:371538/1‑150 (MQ=255)
                               cATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggttg                                                                                                                       >  2:996825/1‑150 (MQ=255)
                               cATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggttg                                                                                                                       >  1:332633/1‑150 (MQ=255)
                               cATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggttg                                                                                                                       >  2:225745/1‑150 (MQ=255)
                                atgatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgt                                                                                                                      <  1:1172562/150‑1 (MQ=255)
                                 tgatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGCCCGCCCTCCcgca                                                                                                                                         <  2:150357/130‑3 (MQ=255)
                                 tgatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGCCCGCCCTCCcgca                                                                                                                                         >  1:150357/1‑128 (MQ=255)
                                    tgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgttgtt                                                                                                                  <  2:659155/150‑1 (MQ=255)
                                    tgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgttgtt                                                                                                                  >  2:1155800/1‑150 (MQ=255)
                                    tgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgttgtt                                                                                                                  <  1:441039/150‑1 (MQ=255)
                                      cAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTTCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTg                                                                                                                >  2:1073066/1‑150 (MQ=255)
                                       aGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTga                                                                                                               >  2:834926/1‑150 (MQ=255)
                                         ccTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTgacg                                                                                                             <  2:1421222/150‑1 (MQ=255)
                                           tCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTgacgac                                                                                                           <  1:284678/150‑1 (MQ=255)
                                           tCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTgacgac                                                                                                           <  2:1059566/150‑1 (MQ=255)
                                            cccTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACa                                                                                                          <  2:274793/150‑1 (MQ=255)
                                              cTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACAtt                                                                                                        >  2:1249343/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)
                                                          cTCCATCGTGGCGGATTAT‑‑TGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTa                                                                                          >  1:382652/1‑150 (MQ=255)
                                                             cATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTAc                                                                                         <  2:394280/150‑1 (MQ=255)
                                                               tCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGc                                                                                       <  2:716461/150‑1 (MQ=255)
                                                                   ggCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCaaa                                                                                   >  2:342988/1‑150 (MQ=255)
                                                                     cGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAAtt                                                                                 >  1:445698/1‑150 (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                                                                   >  2:904300/1‑150 (MQ=255)
                                                                                   gATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTg                                                                   >  1:861506/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)
                                                                                             attaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCAt                                                         >  1:1117218/1‑150 (MQ=255)
                                                                                                   aGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGccc                                                   >  1:1052668/1‑150 (MQ=255)
                                                                                                   aGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGccc                                                   >  2:1179274/1‑150 (MQ=255)
                                                                                                       aCTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaa                                               >  2:1155748/1‑150 (MQ=255)
                                                                                                        cTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaaa                                              >  1:935788/1‑150 (MQ=255)
                                                                                                        cTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaaa                                              <  2:783673/150‑1 (MQ=255)
                                                                                                          aaCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGt                                            >  1:765695/1‑150 (MQ=255)
                                                                                                          aaCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGt                                            >  1:862772/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)
                                                                                                                   agagTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCa                                   >  1:830485/1‑150 (MQ=255)
                                                                                                                         aaCGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCAt                             >  1:495422/1‑150 (MQ=255)
                                                                                                                         aaCGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCAt                             >  2:1254545/1‑150 (MQ=255)
                                                                                                                           cGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCAttt                           >  1:1402095/1‑150 (MQ=255)
                                                                                                                           cGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCAttt                           >  1:499500/1‑150 (MQ=255)
                                                                                                                               aaTTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGAt                       <  1:834926/150‑1 (MQ=255)
                                                                                                                                   cGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGAttttt                   >  2:109929/1‑150 (MQ=255)
                                                                                                                                    gATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTc                  <  1:1238474/150‑1 (MQ=255)
                                                                                                                                    gATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTc                  <  1:694954/150‑1 (MQ=255)
                                                                                                                                      tGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGGCAAAGCCATTGGCCCAGAAAGTTTCCTGCCATGCCATTTCGATTTTTcag                <  2:159979/150‑3 (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)
                                                                                                                                           gTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGcc                                    >  1:12106/1‑125 (MQ=255)
                                                                                                                                           gTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGATGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGc           >  2:326612/1‑150 (MQ=255)
                                                                                                                                           cgccGGGAGGCCCGCCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGc           >  1:325392/4‑150 (MQ=255)
                                                                                                                                               gggAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCa       <  2:1176858/150‑1 (MQ=255)
                                                                                                                                                ggAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCaa      >  1:583670/1‑150 (MQ=255)
                                                                                                                                                ggAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCaa      >  2:1175819/1‑150 (MQ=255)
                                                                                                                                                   gggcgggcCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGAt   <  2:382652/150‑1 (MQ=255)
                                                                                                                                                    ggcgggcCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGAtt  <  1:1328827/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGATT  >  NC_000913/4325132‑4325429

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: