Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,315,907 +A coding (436/759 nt) phnK ← ATP‑binding component of phosphonate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,315,9031.A98.5% 994.9 / 4.7 270coding (440/759 nt)phnKATP‑binding component of phosphonate transport
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base A (147/119);  total (149/121)
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.87e-01

Only 100 of 301 total aligned reads displayed.

CATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCG  >  NC_000913/4315811‑4315989
                                                                                              |                                                                                      
caTCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCgg                                                                                         <  2:537269/93‑1 (MQ=255)
 aTCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGa                                                                                       <  2:6480237/93‑1 (MQ=255)
   ccAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaa                                                                                     <  1:3205881/93‑1 (MQ=255)
     aGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaa                                                                                   <  1:584196/93‑2 (MQ=255)
      gCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaag                                                                                  <  2:6517432/93‑3 (MQ=255)
       cccGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaagg                                                                                 <  2:3065076/93‑4 (MQ=255)
        ccgccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaaggt                                                                                <  2:2480392/93‑5 (MQ=255)
        ccgccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaaggt                                                                                >  2:2935180/1‑89 (MQ=255)
        ccgccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaaggt                                                                                <  1:4634364/93‑5 (MQ=255)
        ccgccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaaggt                                                                                >  1:4338180/1‑89 (MQ=255)
        ccgccgGTCGGTGCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTCCATACCGCCGGaaaaaggt                                                                                <  1:1748408/93‑5 (MQ=255)
           ccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAggtggt                                                                              <  2:5643834/93‑1 (MQ=255)
           ccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCATCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAggtggt                                                                              <  1:5696648/93‑1 (MQ=255)
            cggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTc                                                                             <  2:4582723/93‑1 (MQ=255)
            cggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTc                                                                             <  1:295861/93‑1 (MQ=255)
               tcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGc                                                                          <  2:1727594/93‑1 (MQ=255)
               tcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGc                                                                          <  2:1527556/93‑1 (MQ=255)
               tcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGc                                                                          <  1:6344135/93‑1 (MQ=255)
                cggtACATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATTCCGCCGGAAAAAGGTGGTCGGCa                                                                         <  1:6103590/93‑1 (MQ=255)
                 ggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAg                                                                        >  1:835377/1‑93 (MQ=255)
                  gtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAgg                                                                       >  1:707559/1‑93 (MQ=255)
                  gtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAgg                                                                       >  1:1504889/1‑93 (MQ=255)
                   ttCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgt                                                                      >  2:3579434/1‑93 (MQ=255)
                   ttCATCCATAAACACCAGCTTCGGATGCGTCACCAGGGTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGGGGTCGGCAggc                                                                      >  2:1263832/1‑92 (MQ=255)
                     cATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcg                                                                    <  1:4090630/93‑1 (MQ=255)
                      aTCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgt                                                                   <  1:6308189/93‑1 (MQ=255)
                      aTCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgt                                                                   >  2:4971207/1‑93 (MQ=255)
                       tCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgtc                                                                  >  1:3601853/1‑93 (MQ=255)
                       tCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgtc                                                                  >  1:989310/1‑93 (MQ=255)
                       tCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgtc                                                                  >  1:3706205/1‑93 (MQ=255)
                         cATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGa                                                                >  1:2836714/1‑93 (MQ=255)
                         cATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGa                                                                >  1:656744/1‑93 (MQ=255)
                          aTAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGAt                                                               >  2:1672175/1‑93 (MQ=255)
                           tAAACACCAGCTTTGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATc                                                              <  2:3798348/93‑1 (MQ=255)
                            aaaCACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATcc                                                             >  1:2283177/1‑93 (MQ=255)
                            aaaCACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATcc                                                             <  1:1772851/93‑1 (MQ=255)
                             aaCACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCg                                                            <  1:4811074/93‑1 (MQ=255)
                              acacCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCgg                                                           <  1:1932687/93‑1 (MQ=255)
                               cacCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGt                                                          >  1:5231605/1‑93 (MQ=255)
                               cacCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGt                                                          <  2:225594/93‑1 (MQ=255)
                                acCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGAAAAAAGGTGGTCGGCAGGTCGTCGATCCGGtt                                                         <  1:4545207/93‑1 (MQ=255)
                                acCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGtt                                                         <  1:6453057/93‑1 (MQ=255)
                                  cAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTgg                                                       <  2:3444856/93‑1 (MQ=255)
                                  cAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTgg                                                       >  2:3259392/1‑93 (MQ=255)
                                     cTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCgg                                                    >  2:806960/1‑93 (MQ=255)
                                      ttCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCggg                                                   <  1:4933891/93‑1 (MQ=255)
                                       tCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGa                                                  >  2:3983061/1‑93 (MQ=255)
                                        cGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGaa                                                 <  2:5162752/93‑1 (MQ=255)
                                        cGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGaa                                                 <  1:3936065/93‑1 (MQ=255)
                                         ggATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAt                                                <  2:5260957/93‑1 (MQ=255)
                                          gATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtc                                               >  2:6182487/1‑93 (MQ=255)
                                          gATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtc                                               <  1:3857473/93‑1 (MQ=255)
                                          gATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtc                                               <  1:2665415/93‑1 (MQ=255)
                                           aTGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtct                                              >  2:2579941/1‑93 (MQ=255)
                                             gCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTcc                                            >  2:1827426/1‑93 (MQ=255)
                                               gTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCAc                                          <  1:4163760/93‑1 (MQ=255)
                                               gTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCAc                                          <  1:2768902/93‑1 (MQ=255)
                                                  aCCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctc                                       <  1:5979926/93‑1 (MQ=255)
                                                   ccAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctct                                      <  2:5839828/93‑1 (MQ=255)
                                                   ccAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctct                                      <  2:1788475/93‑1 (MQ=255)
                                                     aGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTc                                    >  1:1335924/1‑93 (MQ=255)
                                                     aGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTc                                    <  1:3157159/93‑1 (MQ=255)
                                                     aGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTc                                    >  1:3302884/1‑93 (MQ=255)
                                                     aGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACATCTTc                                    >  2:613158/1‑93 (MQ=255)
                                                      ggTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTcc                                   <  1:3018015/93‑1 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGTCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                  >  2:1311372/1‑93 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                  <  2:2029790/93‑1 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                  <  1:3404070/93‑1 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                  <  2:4929089/93‑1 (MQ=255)
                                                        ttGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAg                                 <  1:5348582/93‑1 (MQ=255)
                                                         tGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGc                                <  2:661806/93‑1 (MQ=255)
                                                         tGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGc                                <  2:5418791/93‑1 (MQ=255)
                                                          gCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGcc                               <  1:3539438/93‑1 (MQ=255)
                                                            gggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCAc                             >  2:6288740/1‑93 (MQ=255)
                                                            gggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGCCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCATCCAc                             >  2:1468112/1‑93 (MQ=255)
                                                             ggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACt                            >  2:596410/1‑93 (MQ=255)
                                                             ggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACt                            >  1:221545/1‑93 (MQ=255)
                                                              gCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACtt                           >  1:4754377/1‑93 (MQ=255)
                                                                aaTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCt                         >  1:122542/1‑93 (MQ=255)
                                                                aaTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCt                         <  2:4686521/93‑1 (MQ=255)
                                                                 aTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTg                        <  2:5197232/93‑1 (MQ=255)
                                                                  tCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgc                       >  1:4474809/1‑93 (MQ=255)
                                                                  tCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgc                       >  2:4061779/1‑93 (MQ=255)
                                                                   cTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgcg                      >  1:1593496/1‑93 (MQ=255)
                                                                   cTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgcg                      >  1:1065426/1‑93 (MQ=255)
                                                                       aaaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGt                  >  1:2022223/1‑93 (MQ=255)
                                                                       aaaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGt                  >  1:6639880/1‑93 (MQ=255)
                                                                        aaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTg                 >  2:3865090/1‑93 (MQ=255)
                                                                        aaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTg                 >  2:703523/1‑93 (MQ=255)
                                                                        aaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGAAATCTCCACCTCTTCCAGCGACTTCTGGGCGGTg                 >  2:4734649/1‑93 (MQ=255)
                                                                         aCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTgg                <  2:4202993/93‑1 (MQ=255)
                                                                         aCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTgg                <  1:1312336/93‑1 (MQ=255)
                                                                              gctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGa           >  1:2156863/1‑93 (MQ=255)
                                                                              gctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGa           <  2:3617345/93‑1 (MQ=255)
                                                                              gctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGa           <  2:1901380/93‑1 (MQ=255)
                                                                                tgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAat         >  2:2825357/1‑93 (MQ=255)
                                                                                tgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAat         >  1:3576745/1‑93 (MQ=255)
                                                                                 gcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAata        >  2:1401501/1‑93 (MQ=255)
                                                                                   aTACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATc      >  1:2540112/1‑93 (MQ=255)
                                                                                   aTACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATc      <  2:2287272/93‑1 (MQ=255)
                                                                                              |                                                                                      
CATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCG  >  NC_000913/4315811‑4315989

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: