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.A97.8% 852.9 / 7.6 233coding (440/759 nt)phnKATP‑binding component of phosphonate transport
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base A (121/106);  total (124/109)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 263 total aligned reads displayed.

CATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGT  >  NC_000913/4315811‑4315990
                                                                                              |                                                                                       
   ccAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaa                                                                                      >  1:5897521/1‑93 (MQ=255)
      gCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaag                                                                                   <  1:4599257/93‑3 (MQ=255)
      gCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaag                                                                                   >  2:3576526/1‑91 (MQ=255)
      ccccGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaag                                                                                   <  2:4478467/92‑3 (MQ=255)
        ccgcggGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCATACGCTGCTGCATACCGCCGGaaaaaggt                                                                                 <  2:4016947/93‑5 (MQ=255)
        ccgccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaaggt                                                                                 >  1:1356786/1‑89 (MQ=255)
        ccgccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaaggt                                                                                 <  2:4553020/93‑5 (MQ=255)
          gccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAggtgg                                                                                <  2:3852899/93‑1 (MQ=255)
          gccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAggtgg                                                                                >  2:2165300/1‑93 (MQ=255)
          gccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAggtgg                                                                                >  1:2487432/1‑93 (MQ=255)
            cggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTc                                                                              <  1:1196694/93‑1 (MQ=255)
            cggtcggtTCATCCATAAACACCAGCTACGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTc                                                                              <  2:2958203/93‑1 (MQ=255)
             ggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCg                                                                             <  2:5717367/93‑1 (MQ=255)
              gtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCgg                                                                            >  2:5495748/1‑93 (MQ=255)
               tcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGc                                                                           <  2:4121723/93‑1 (MQ=255)
               tcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGc                                                                           <  2:4265324/93‑1 (MQ=255)
                cggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCa                                                                          >  1:5190046/1‑93 (MQ=255)
                  gtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAgg                                                                        >  2:72104/1‑93 (MQ=255)
                  gtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAgg                                                                        <  2:1776600/93‑1 (MQ=255)
                   ttCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAATGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgt                                                                       >  2:1453066/1‑93 (MQ=255)
                   ttCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGGGGTCGGCAGgt                                                                       >  2:4394414/1‑93 (MQ=255)
                      aTCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgt                                                                    <  2:3760471/93‑1 (MQ=255)
                      aTCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgt                                                                    <  1:5799148/93‑1 (MQ=255)
                       tCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgtc                                                                   <  2:860995/93‑1 (MQ=255)
                       tCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgtc                                                                   >  1:3585640/1‑93 (MQ=255)
                       tCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgtc                                                                   <  2:2563403/93‑1 (MQ=255)
                         cATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGTAAAAGGTGGTCGGCAGGTCGTCGa                                                                 >  2:328477/1‑93 (MQ=255)
                           tAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATc                                                               <  2:9384/93‑1 (MQ=255)
                           tAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATc                                                               >  1:4528794/1‑93 (MQ=255)
                            aaaCACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATcc                                                              <  2:3191438/93‑1 (MQ=255)
                                acCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGtt                                                          >  2:1838499/1‑93 (MQ=255)
                                  cAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTgg                                                        <  2:5620518/93‑1 (MQ=255)
                                  cAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTgg                                                        >  1:913801/1‑93 (MQ=255)
                                    gCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCg                                                      <  1:1741284/93‑1 (MQ=255)
                                         ggATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAt                                                 <  2:3217556/93‑1 (MQ=255)
                                          gATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtc                                                >  2:816477/1‑93 (MQ=255)
                                          gATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtc                                                >  2:3377592/1‑93 (MQ=255)
                                            tGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtctc                                              <  1:4407030/93‑1 (MQ=255)
                                            tGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtctc                                              <  2:1410564/93‑1 (MQ=255)
                                             gCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTcc                                             <  2:2745641/93‑1 (MQ=255)
                                              cGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCa                                            <  2:4839293/93‑1 (MQ=255)
                                              cGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCa                                            <  1:2335142/93‑1 (MQ=255)
                                                tCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCAcc                                          >  2:5942019/1‑93 (MQ=255)
                                                tCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCAcc                                          >  1:3441642/1‑93 (MQ=255)
                                                  aCCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctc                                        >  2:1139005/1‑93 (MQ=255)
                                                  aCCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctc                                        <  2:6274325/93‑1 (MQ=255)
                                                   ccAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctct                                       <  2:3474166/93‑1 (MQ=255)
                                                      ggTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTcc                                    >  1:2543116/1‑93 (MQ=255)
                                                      ggTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTcc                                    >  1:5633417/1‑93 (MQ=255)
                                                      ggTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTcc                                    >  1:5564038/1‑93 (MQ=255)
                                                      ggTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTcc                                    >  2:5656523/1‑93 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                   <  1:4172658/93‑1 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                   >  1:1641103/1‑93 (MQ=255)
                                                       gTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCa                                   >  1:3340238/1‑93 (MQ=255)
                                                        ttGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAg                                  >  1:1067447/1‑93 (MQ=255)
                                                        ttGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAg                                  >  2:1063704/1‑93 (MQ=255)
                                                           cGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCa                               >  2:2379680/1‑93 (MQ=255)
                                                           cGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCa                               <  1:1977057/93‑1 (MQ=255)
                                                            gggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCCCCTCTTCCAGCCAc                              <  1:2862543/93‑1 (MQ=255)
                                                            gggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCAc                              >  2:2983157/1‑93 (MQ=255)
                                                             ggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACt                             >  1:2236076/1‑93 (MQ=255)
                                                             ggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACt                             >  1:4312009/1‑93 (MQ=255)
                                                             ggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGGTGGCGGGAATCTCCACCTCTTCCAGCCACt                             >  1:5061168/1‑93 (MQ=255)
                                                               cAATCTGCAAACGCTGCTGCATCCCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTc                           <  1:3918484/93‑1 (MQ=255)
                                                               cAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTc                           >  1:109394/1‑93 (MQ=255)
                                                               cAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTc                           <  2:2491755/93‑1 (MQ=255)
                                                                aaTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCt                          <  1:953653/93‑1 (MQ=255)
                                                                 aTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTg                         <  2:1294620/93‑1 (MQ=255)
                                                                  tCTGCAACCGCTGCTGCATCCCGCGGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGAAATCCCCCCCTCTCCCAGCCACTTCTgc                        <  2:4992982/93‑1 (MQ=255)
                                                                  tCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgc                        >  2:5332427/1‑93 (MQ=255)
                                                                  tCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgc                        <  1:6109469/93‑1 (MQ=255)
                                                                  tCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgc                        <  2:2455268/93‑1 (MQ=255)
                                                                    tGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgcgc                      >  1:500183/1‑93 (MQ=255)
                                                                      cAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCgg                    <  2:2841388/93‑1 (MQ=255)
                                                                      cAAACGCGGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCgg                    <  1:449684/93‑1 (MQ=255)
                                                                       aaaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGGCACTTCTGCGCggg                   >  2:5050203/1‑92 (MQ=255)
                                                                         aCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTgg                 <  1:4090771/93‑1 (MQ=255)
                                                                           gctgctGCATACCGCCGAAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCa               <  2:6041036/93‑1 (MQ=255)
                                                                           gctgctGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCa               <  2:6178620/93‑1 (MQ=255)
                                                                           gctgctGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCa               <  2:5151395/93‑1 (MQ=255)
                                                                           gctgctGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCa               <  2:4072515/93‑1 (MQ=255)
                                                                            ctgctgCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCAc              <  1:3958614/93‑1 (MQ=255)
                                                                            ctgctgCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCAc              >  1:2628800/1‑93 (MQ=255)
                                                                             tgctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACg             <  1:6193535/93‑1 (MQ=255)
                                                                             tgctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACg             >  1:3252707/1‑93 (MQ=255)
                                                                             tgctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCAACCTCTTCCAGCCACTTCTGCGCGGTGGCACg             >  1:1928581/1‑93 (MQ=255)
                                                                               ctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGaa           <  2:5825377/93‑1 (MQ=255)
                                                                                tgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAat          >  2:1580291/1‑93 (MQ=255)
                                                                                tgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGCGGCACGAat          >  1:4923576/1‑93 (MQ=255)
                                                                                 gcATGCCGCCGGAAAAAGGTGGTCGGCAGGTCGTCTATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAata         <  1:5532919/93‑1 (MQ=255)
                                                                                 gcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAata         <  2:2070110/93‑1 (MQ=255)
                                                                                  cATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAatat        <  2:1250131/93‑1 (MQ=255)
                                                                                   aTACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATc       >  1:2994039/1‑93 (MQ=255)
                                                                                   aTACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATc       >  1:4645889/1‑93 (MQ=255)
                                                                                   aTACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATc       <  1:759944/93‑1 (MQ=255)
                                                                                    tACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCg      >  2:169869/1‑93 (MQ=255)
                                                                                     aCCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGc     >  1:1627/1‑93 (MQ=255)
                                                                                     aCCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGc     <  1:1532433/93‑1 (MQ=255)
                                                                                      ccgccgGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGcc    >  2:4887792/1‑93 (MQ=255)
                                                                                        gccgGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGt  <  2:1636942/93‑1 (MQ=255)
                                                                                              |                                                                                       
CATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGT  >  NC_000913/4315811‑4315990

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