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.A99.6% 915.7 / ‑3.8 246coding (440/759 nt)phnKATP‑binding component of phosphonate transport
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (123/122);  total (123/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

Only 100 of 279 total aligned reads displayed.

CATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGT  >  NC_000913/4315811‑4315990
                                                                                             |                                                                                       
 aTCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCAGCAGGTTGCCGGCAATCTGCAAACGCCGCTGCATACCGCCGGa                                                                                        >  2:5266185/1‑93 (MQ=255)
     aGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaa                                                                                    <  1:7098185/93‑2 (MQ=255)
     aGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaa                                                                                    >  2:159608/1‑92 (MQ=255)
      gCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaag                                                                                   >  2:2499073/1‑91 (MQ=255)
      gCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaag                                                                                   <  2:1777780/93‑3 (MQ=255)
      gCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaag                                                                                   >  1:3248278/1‑91 (MQ=255)
      gCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaag                                                                                   >  2:3663029/1‑91 (MQ=255)
       cccGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaagg                                                                                  <  1:8449905/93‑4 (MQ=255)
        ccgccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaaggt                                                                                 <  1:6425338/93‑5 (MQ=255)
        ccgccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaaggt                                                                                 <  1:8146432/93‑5 (MQ=255)
        ccgccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTCCATACCGCCGGaaaaaggg                                                                                 >  1:382823/1‑89 (MQ=255)
          gccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAggtgg                                                                                >  2:3857402/1‑93 (MQ=255)
            cggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTc                                                                              <  1:8865283/93‑1 (MQ=255)
            cggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTc                                                                              <  2:3303896/93‑1 (MQ=255)
             ggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCg                                                                             >  2:3953824/1‑93 (MQ=255)
               tcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGc                                                                           <  2:7138109/93‑1 (MQ=255)
               tcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACACTGCTGCATACCGCCGGAAAAAGGTGGCCGGc                                                                           >  1:1412147/1‑93 (MQ=255)
                cggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCa                                                                          <  2:3007987/93‑1 (MQ=255)
                 ggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAg                                                                         >  2:4766877/1‑93 (MQ=255)
                 ggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAg                                                                         <  1:9811104/93‑1 (MQ=255)
                  gtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAgg                                                                        >  1:2045354/1‑93 (MQ=255)
                   ttCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgt                                                                       >  2:6113890/1‑93 (MQ=255)
                   ttCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgt                                                                       >  1:4139822/1‑93 (MQ=255)
                     cATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcg                                                                     <  2:3574242/93‑1 (MQ=255)
                     cATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcg                                                                     >  1:1881618/1‑93 (MQ=255)
                     cATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcg                                                                     >  1:4550306/1‑93 (MQ=255)
                       tCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCACTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgtg                                                                   >  2:8539841/1‑92 (MQ=255)
                       tCCATAAACACCAGCTTCGGATGCATCACCAGGTTGCGGGCAATCTGCAAACGCTGTTGCATACCGCCGGAAAAAGCTGGTCGGCAGctcgtc                                                                   >  1:445185/1‑93 (MQ=255)
                         cATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGa                                                                 >  1:4332503/1‑93 (MQ=255)
                         cATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGa                                                                 <  1:3282378/93‑1 (MQ=255)
                           tAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATc                                                               <  1:5501414/93‑1 (MQ=255)
                           aaaaCACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATc                                                               <  2:8715284/92‑1 (MQ=255)
                            aaaCACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATGCCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATcc                                                              <  2:4574501/93‑1 (MQ=255)
                             aaCACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCg                                                             >  2:2688169/1‑93 (MQ=255)
                                acCATCCTCGGATGCGTCACCAGGTTGCGCGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTg                                                         <  1:9530501/93‑1 (MQ=255)
                                acCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGtt                                                          <  1:6805691/93‑1 (MQ=255)
                                   aGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGc                                                       >  1:4314832/1‑93 (MQ=255)
                                        cGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGaa                                                  <  1:5175371/93‑1 (MQ=255)
                                          gATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtc                                                <  1:5206131/93‑1 (MQ=255)
                                          gATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtc                                                >  1:517959/1‑93 (MQ=255)
                                          gATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtc                                                >  1:832581/1‑93 (MQ=255)
                                             gCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTcc                                             >  2:1406063/1‑93 (MQ=255)
                                              cGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCGGGTCGTCGATCCGGTTGGCGGGAATCTCCa                                            <  1:1073900/93‑1 (MQ=255)
                                              cGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCa                                            <  2:6003087/93‑1 (MQ=255)
                                               gTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCAc                                           >  1:1022422/1‑93 (MQ=255)
                                               gTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCAc                                           >  1:1618958/1‑93 (MQ=255)
                                                tCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCAcc                                          >  2:4200725/1‑93 (MQ=255)
                                                 cACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACct                                         >  1:9922398/1‑93 (MQ=255)
                                                 cACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACct                                         >  2:8271485/1‑93 (MQ=255)
                                                  aCCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctc                                        >  2:2960814/1‑93 (MQ=255)
                                                  aCCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctc                                        <  1:5948199/93‑1 (MQ=255)
                                                   ccAGGGTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctct                                       <  1:1548625/93‑1 (MQ=255)
                                                     aGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTc                                     <  1:9996809/93‑1 (MQ=255)
                                                      ggTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTcc                                    <  2:8848801/93‑1 (MQ=255)
                                                      ggTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTcc                                    <  2:2103534/93‑1 (MQ=255)
                                                       gtggCGGGCAATCTGCAAACGCTGCTCCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                   <  2:3761355/90‑1 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGGATCTCCACCTCTTCCa                                   >  1:6183969/1‑93 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                   >  2:4335675/1‑93 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                   >  2:8156574/1‑93 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                   >  1:7629136/1‑93 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                   >  1:2077664/1‑93 (MQ=255)
                                                           cGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCa                               <  1:7506938/93‑1 (MQ=255)
                                                            gggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCAc                              >  2:7915634/1‑93 (MQ=255)
                                                               cAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTc                           <  2:96676/93‑1 (MQ=255)
                                                               cAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTc                           <  1:4276484/93‑1 (MQ=255)
                                                               cAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGGTGGCGGGAATCTCCACCTCTTCCAGCCACTTc                           >  2:7756986/1‑93 (MQ=255)
                                                                aaTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCt                          <  2:6280831/93‑1 (MQ=255)
                                                                aaTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCt                          >  1:8018352/1‑93 (MQ=255)
                                                                  tCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgc                        >  1:3176541/1‑93 (MQ=255)
                                                                  tCTGCAAACGCGGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgc                        <  2:4280302/93‑1 (MQ=255)
                                                                    tGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgcgc                      >  2:6848402/1‑93 (MQ=255)
                                                                     gCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGcgcg                     >  1:658603/1‑93 (MQ=255)
                                                                     gCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGcgcg                     <  1:8693476/93‑1 (MQ=255)
                                                                      cAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCgg                    <  1:4765802/93‑1 (MQ=255)
                                                                       aaaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCggg                   >  2:9279306/1‑92 (MQ=255)
                                                                       aaaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGt                   >  2:9085029/1‑93 (MQ=255)
                                                                       aaaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGt                   >  2:5142965/1‑93 (MQ=255)
                                                                       aaaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGt                   <  2:1214557/93‑1 (MQ=255)
                                                                        aaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTg                  >  1:2690989/1‑93 (MQ=255)
                                                                        aaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTg                  >  2:7397513/1‑93 (MQ=255)
                                                                        aaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTCCTGCGCGGTg                  <  2:6548287/93‑1 (MQ=255)
                                                                            ctgctgCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCAc              >  2:7685947/1‑93 (MQ=255)
                                                                            ctgctgCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCAc              >  2:4664397/1‑93 (MQ=255)
                                                                             tgctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGACACTTCTGCGCGGTGGCACg             >  1:254856/1‑93 (MQ=255)
                                                                              gctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGa            >  2:4922476/1‑93 (MQ=255)
                                                                               ctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGaa           <  1:2859978/93‑1 (MQ=255)
                                                                                 gcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAata         <  2:3471270/93‑1 (MQ=255)
                                                                                 gcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAata         >  1:2956718/1‑93 (MQ=255)
                                                                                 gcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAata         >  2:8052832/1‑93 (MQ=255)
                                                                                  cATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAatat        <  1:9442273/93‑1 (MQ=255)
                                                                                  cATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAatat        >  2:3279208/1‑93 (MQ=255)
                                                                                   aTACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATc       <  2:6993314/93‑1 (MQ=255)
                                                                                   aTACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATc       >  2:349193/1‑93 (MQ=255)
                                                                                   aTACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATc       <  1:5796214/93‑1 (MQ=255)
                                                                                   aTACCGCCGGAAAAAGGTGGTCGGCAGGGCGTCGATCCGGGTGGCGGGAATCTGCACCTCTTCCAGCCACTTCTGCGCGGCGGCACCCATATc       >  2:550046/1‑93 (MQ=255)
                                                                                    tACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCg      <  1:2289895/93‑1 (MQ=255)
                                                                                     aCCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGc     <  1:9114343/93‑1 (MQ=255)
                                                                                     aCCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGc     <  2:522968/93‑1 (MQ=255)
                                                                                     aCCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGc     >  2:9402303/1‑93 (MQ=255)
                                                                                        gccgGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGGATATCGCCGt  >  2:200520/1‑93 (MQ=255)
                                                                                             |                                                                                       
CATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGT  >  NC_000913/4315811‑4315990

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