Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,315,100 G→C 6.9% intergenic (‑118/+15) phnC ← / ← yjdN phosphonate ABC transporter ATP‑binding subunit/metalloprotein superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,315,1000GC 6.9% 179.6 / 3.6 61intergenic (‑118/+15)phnC/yjdNphosphonate ABC transporter ATP‑binding subunit/metalloprotein superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base C (2/3);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGATT  >  CP009273/4314952‑4315246
                                                                                                                                                    |                                                                                                                                                  
ccaGCTTCTCGACACGGATAATCGGTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAgggcg                                                                                                                                                    <  2:305315/147‑1 (MQ=255)
 gAGCTTCTCGACCCGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAgggcgg                                                                                                                                                   <  1:172374/149‑1 (MQ=255)
    cTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTTGATGACAGTGCGGGAGGGCGGGcc                                                                                                                                                <  1:628237/149‑1 (MQ=255)
     ttCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCt                                                                                                                                               >  1:616532/1‑149 (MQ=255)
     ttCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCAATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGACGAGTTAACGGGAATTCGGTGACGGTGCGGGCGGGCCGGCCt                                                                                                                                               >  1:103597/1‑149 (MQ=255)
          gACACGGATAATGGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGCCCGCCCTCCcgc                                                                                                                                          <  1:184820/149‑2 (MQ=255)
          gACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTACAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCcgg                                                                                                                                          <  1:1221442/149‑1 (MQ=255)
          gACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCcgg                                                                                                                                          >  2:974131/1‑149 (MQ=255)
          gACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAAAAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCcgg                                                                                                                                          >  2:220568/1‑149 (MQ=255)
          gACACGGATAATCGTTTGCATGATCCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGCCCGCCCTCCcgc                                                                                                                                          <  2:637047/149‑2 (MQ=255)
          gACAAGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGCCCGCCCTCCcgc                                                                                                                                          <  2:796269/149‑2 (MQ=255)
                gATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCCTCGTGGCCGATTATTGGGACGATTTGGCTACTTAAAAGCGACTAACAGATGAAGACTTGACAGGAATTCGACGTCCGTGCGGCAGAGCGCGCCTCGCGGCggggt                                                                                                                                    >  2:110742/1‑147 (MQ=255)
                 aTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTa                                                                                                                                   >  2:538789/1‑149 (MQ=255)
                 aTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTa                                                                                                                                   >  2:126476/1‑149 (MQ=255)
                    aTCGTTTGCATCATGCAGCCGCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTACCAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTg                                                                                                                                <  1:684832/149‑1 (MQ=255)
                          tGCATGATGCCGCCTCCCTAATACAGTGGCTCCATCGTAGCAGCATATTGTTACGACCGGGTTAATTAAAAGTGACTAACAGATGAAGAGTTACCTGGAACTCCATGACAAGGTGGGAGGGCTGGCCTCCCGGCGGGTAATTGCgttgg                                                                                                                          <  2:529494/149‑1 (MQ=255)
                          tGCAGGATGCCGCCCCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttgg                                                                                                                          <  2:463749/149‑1 (MQ=255)
                            cATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGTAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggtt                                                                                                                        <  1:132123/149‑1 (MQ=255)
                            aatgatgCAGCCTCCCTAAAAAAGTGGCTCCGTCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCgttggtt                                                                                                                        >  2:906489/2‑149 (MQ=255)
                              tgatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGTTTTATAAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTAGAGGACAGTGCGGGAGGGCCGGCCTCCAGGCGGGTTAGTGCGTAGgttgt                                                                                                                      >  2:1025389/1‑149 (MQ=255)
                              tgatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCCCCCGGCGGGGTATTGCGCTGgttgt                                                                                                                      >  1:308388/1‑149 (MQ=255)
                               gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTTACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgtt                                                                                                                     <  1:1150165/149‑1 (MQ=255)
                                 tgCAGCCTCCCTAAACAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgttgt                                                                                                                   <  2:976739/149‑1 (MQ=255)
                                 tgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGgttgttgt                                                                                                                   <  2:662408/149‑1 (MQ=255)
                                  gCAGCCTCCCTAAAAAAGTGGCTCCATCATGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCAATGACAGTGCGGGAGGGAGGGCCTCCCGGCGGGTTATTGCGTTGgttgttgtt                                                                                                                  >  1:698190/1‑149 (MQ=255)
                                   cAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGttt                                                                                                                 <  1:1134669/149‑1 (MQ=255)
                                   cAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGGGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGttt                                                                                                                 <  2:681030/149‑1 (MQ=255)
                                     gCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTga                                                                                                               <  2:1209278/149‑1 (MQ=255)
                                     gCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTga                                                                                                               <  2:1209045/149‑1 (MQ=255)
                                     gCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAAATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGGTGTTGTTTga                                                                                                               >  2:1202784/1‑149 (MQ=255)
                                          ccTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACa                                                                                                          <  2:207534/149‑1 (MQ=255)
                                             aaaaaaGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGTCTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGACTCCCAGCGGGTTATTGCGTTGGTTGTAGTTTGACGACATTa                                                                                                       >  1:874435/1‑149 (MQ=255)
                                              aaaaaGTGGCTCCATCGTGGCTGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTaa                                                                                                      <  1:556662/149‑1 (MQ=255)
                                                   gTGGCCCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATGCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAatcatc                                                                                                 <  1:116725/149‑1 (MQ=255)
                                                       cTCCATCGTGGCGGATTATTGTGCCGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACg                                                                                             >  1:688261/1‑149 (MQ=255)
                                                       cTCCATCGTGGCGAATTATTGTGACGATTTGGTTAATTACAAGTGACTAACAGATGAAGAGTAAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGGTGTTGTTTGACGACATTAATCATCCCCg                                                                                             >  2:760554/1‑149 (MQ=255)
                                                        tCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACgg                                                                                            <  2:176146/149‑1 (MQ=255)
                                                        tCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCAGCCACgg                                                                                            >  1:332224/1‑149 (MQ=255)
                                                          cATCGTGGCGGATTATTGTGAAGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGCCGCCATTAATGCTCCACGGTa                                                                                          >  2:168995/1‑149 (MQ=255)
                                                                 gCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTGGTTGTTTGACGCCATTAATCATCCACGGTACGCCaag                                                                                   >  1:209890/1‑148 (MQ=255)
                                                                            tgACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGa                                                                        <  1:771645/149‑1 (MQ=255)
                                                                            tgACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGa                                                                        >  1:17548/1‑149 (MQ=255)
                                                                               ctaTTTGGTTATTTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACtt                                                                     <  2:637437/147‑1 (MQ=255)
                                                                                    tGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCa                                                                >  2:550114/1‑149 (MQ=255)
                                                                                        taattaaAAGTGACTAACAGATGAAGGGTGAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGc                                                            <  1:232786/149‑1 (MQ=255)
                                                                                            taaAAGTGACTGACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATg                                                        <  2:690399/149‑1 (MQ=255)
                                                                                               aaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATgcgc                                                     <  1:538789/149‑1 (MQ=255)
                                                                                               aaGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATgcgc                                                     <  1:126476/149‑1 (MQ=255)
                                                                                                  tGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGTCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCa                                                  <  1:958039/149‑1 (MQ=255)
                                                                                                      tAACAGATGAAGAGTTAACGGGAATTCGGTGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCAAGCGCCCAGaaa                                              >  2:697246/1‑149 (MQ=255)
                                                                                                      tAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGGTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGaaa                                              >  1:459527/1‑149 (MQ=255)
                                                                                                              gAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTccag                                      >  1:429617/1‑147 (MQ=255)
                                                                                                               aaGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGc                                     >  1:109479/1‑149 (MQ=255)
                                                                                                                    ttAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCGTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCGAAATTTATCGGTGGCTTTGCCAAGGCCATGCGCGCAGAAAGTTTCCTGCCAGGc                                >  2:231651/1‑149 (MQ=255)
                                                                                                                      aaCGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAGAGCCATGCGCCCAGAAAGTTTCCTCCCAGGCCa                              <  1:906489/149‑1 (MQ=255)
                                                                                                                          ggAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTc                          >  2:1204795/1‑149 (MQ=255)
                                                                                                                            aaTTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTCGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGa                        <  2:889247/149‑1 (MQ=255)
                                                                                                                                  aTGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTc                  >  2:362696/1‑149 (MQ=255)
                                                                                                                                       agcgCGGGAGGGTGCGCCTCCGCGGGGGTTATTGCGTTGGTTGCTGCTTGACGACATTAATCATCCTCGGTACGCCAAACTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTa             <  1:576033/146‑1 (MQ=255)
                                                                                                                                         tGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCCGGCCATTTCGATTTTTCCGTTAGc           >  2:203486/1‑149 (MQ=255)
                                                                                                                                         tGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGc           >  1:849670/1‑149 (MQ=255)
                                                                                                                                          gCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCg          <  1:47274/149‑1 (MQ=255)
                                                                                                                                           cGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCgg         <  1:1003095/149‑1 (MQ=255)
                                                                                                                                                 ggcgggcCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGAt   >  1:254033/1‑149 (MQ=255)
                                                                                                                                                  gcgggcCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCCTTTCGATTTTTCCGTTAGCGGCCAGAtt  >  2:156499/1‑149 (MQ=255)
                                                                                                                                                  gcgggcCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATGTCGATTTTTCCGTTAGCGGCAAGAtt  >  1:1150397/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCTCCCGGCGGGTTATTGCGTTGGTTGTTGTTTGACGACATTAATCATCCACGGTACGCCAAATTTATCGGTGACTTTGCCAAAGCCATGCGCCCAGAAAGTTTCCTGCCAGGCCATTTCGATTTTTCCGTTAGCGGCAAGATT  >  CP009273/4314952‑4315246

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