Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,325,283 G→C 7.1% intergenic (‑118/+15) 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,2830GC 7.1% 418.2 / 16.1 125intergenic (‑118/+15)phnC/yjdNphosphonate ABC transporter ATP binding subunit/conserved protein YjdN
Reads supporting (aligned to +/- strand):  ref base G (70/46);  new base C (4/5);  total (74/51)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

Only 100 of 130 total aligned reads displayed.

CGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGATT  >  NC_000913/4325135‑4325429
                                                                                                                                                    |                                                                                                                                                  
cGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAgggcgg                                                                                                                                                   >  1:1052970/1‑150 (MQ=255)
 gAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAgggcggg                                                                                                                                                  <  1:1151384/150‑1 (MQ=255)
 gAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAgggcggg                                                                                                                                                  <  1:1362601/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)
         cGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTTGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGCCCGCCCTCCcgc                                                                                                                                          <  2:476278/150‑2 (MQ=255)
         cGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCcgg                                                                                                                                          >  1:394280/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)
             acGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTCGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCggg                                                                                                                                      <  1:233060/150‑1 (MQ=255)
                  tAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTAtt                                                                                                                                 >  2:494265/1‑150 (MQ=255)
                    aTCGTCTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGc                                                                                                                               <  1:597561/150‑1 (MQ=255)
                       gTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgtt                                                                                                                            <  1:1168476/150‑1 (MQ=255)
                       gTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgtt                                                                                                                            >  1:173974/1‑150 (MQ=255)
                        tttGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCTTCCCGGCGGGTTATTGCgttg                                                                                                                           <  2:1338815/150‑1 (MQ=255)
                         ttGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttgg                                                                                                                          <  2:20100/150‑1 (MQ=255)
                         ttGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttgg                                                                                                                          <  2:724561/150‑1 (MQ=255)
                         ttGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttgg                                                                                                                          >  2:1328827/1‑150 (MQ=255)
                           gCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggtt                                                                                                                        >  1:1176858/1‑150 (MQ=255)
                            cATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggttg                                                                                                                       >  2:225745/1‑150 (MQ=255)
                            cATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggttg                                                                                                                       >  1:332633/1‑150 (MQ=255)
                            cATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggttg                                                                                                                       >  2:996825/1‑150 (MQ=255)
                             atgatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgt                                                                                                                      <  1:1172562/150‑1 (MQ=255)
                              tgatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGCCCGCCCTCCcgca                                                                                                                                         >  1:150357/1‑128 (MQ=255)
                               gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgttg                                                                                                                    >  1:860181/1‑150 (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)
                                      ccTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTgacg                                                                                                             <  2:1421222/150‑1 (MQ=255)
                                        tCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTgacgac                                                                                                           <  2:1059566/150‑1 (MQ=255)
                                        tCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTgacgac                                                                                                           <  1:484528/150‑1 (MQ=255)
                                         cccTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACa                                                                                                          <  2:274793/150‑1 (MQ=255)
                                           cTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACAtt                                                                                                        >  2:1249343/1‑150 (MQ=255)
                                              aaaaaGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAat                                                                                                     >  2:1275840/1‑150 (MQ=255)
                                              aaaaaGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAat                                                                                                     >  2:599031/1‑150 (MQ=255)
                                                      gCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGATGTTGTTTTACGACATTAATCATCCACg                                                                                             >  1:315740/1‑150 (MQ=255)
                                                       cTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACgg                                                                                            <  1:716760/150‑1 (MQ=255)
                                                       cTCCATCGTGGCGGATTAT‑‑TGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTa                                                                                          >  1:382652/1‑150 (MQ=255)
                                                        tCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGt                                                                                           >  1:194292/1‑150 (MQ=255)
                                                          cATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTAc                                                                                         <  2:394280/150‑1 (MQ=255)
                                                           aTCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACg                                                                                        >  1:662250/1‑150 (MQ=255)
                                                            tCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGc                                                                                       <  2:716461/150‑1 (MQ=255)
                                                              gTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGAAGACATTAATCATCCACGGTACGCCa                                                                                     >  1:1266913/1‑150 (MQ=255)
                                                               tGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTAGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCaa                                                                                    >  1:524398/1‑150 (MQ=255)
                                                                ggCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCaaa                                                                                   >  2:342988/1‑150 (MQ=255)
                                                                  cGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAAtt                                                                                 >  1:445698/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)
                                                                           gtgACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGa                                                                        >  2:26594/1‑150 (MQ=255)
                                                                            tgACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGAc                                                                       >  1:1064980/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)
                                                                                  tttGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGcc                                                                 <  2:679914/150‑1 (MQ=255)
                                                                                     ggTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCaaa                                                              >  2:757831/1‑150 (MQ=255)
                                                                                     ggTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCaaa                                                              >  1:16242/1‑150 (MQ=255)
                                                                                      gTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAg                                                             >  2:270453/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)
                                                                                          attaaAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCAt                                                         >  1:1117218/1‑150 (MQ=255)
                                                                                                aGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGccc                                                   >  1:1052668/1‑150 (MQ=255)
                                                                                                aGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGccc                                                   >  2:1179274/1‑150 (MQ=255)
                                                                                                   gACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGa                                                >  2:185582/1‑150 (MQ=255)
                                                                                                    aCTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaa                                               >  2:1155748/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)
                                                                                                       aaCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGt                                            >  1:862772/1‑150 (MQ=255)
                                                                                                         cAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGttt                                          >  1:1028700/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)
                                                                                                                            aaTTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGAt                       <  1:834926/150‑1 (MQ=255)
                                                                                                                              ttCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGAttt                     <  1:788580/150‑1 (MQ=255)
                                                                                                                                 gATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTc                  <  1:694954/150‑1 (MQ=255)
                                                                                                                                 gATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTc                  <  1:1238474/150‑1 (MQ=255)
                                                                                                                                 gATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTc                  >  2:1024858/1‑150 (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                                    <  2:12106/125‑1 (MQ=255)
                                                                                                                                        gTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGcc                                    >  1:12106/1‑125 (MQ=255)
                                                                                                                                        cgccGGGAGGCCCGCCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGc           >  1:325392/4‑150 (MQ=255)
                                                                                                                                        cgccGGGAGGCCCGCCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGc           <  2:1157143/147‑1 (MQ=255)
                                                                                                                                           cGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGc        >  2:737453/1‑150 (MQ=255)
                                                                                                                                            gggAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCa       <  2:1176858/150‑1 (MQ=255)
                                                                                                                                             ggAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCaa      >  1:583670/1‑150 (MQ=255)
                                                                                                                                                gggcgggcCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGAt   <  2:382652/150‑1 (MQ=255)
                                                                                                                                                 ggcgggcCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGAtt  <  1:1328827/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGATT  >  NC_000913/4325135‑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: