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.A96.2% 441.4 / 7.1 130coding (440/759 nt)phnKATP‑binding component of phosphonate transport
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base A (67/58);  total (71/59)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

CATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCG  >  NC_000913/4315811‑4315989
                                                                                             |                                                                                        
cttCCAGCCCGCCGGTCGGTTCATCCATAAACCCCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCgg                                                                                           <  1:1874527/91‑1 (MQ=255)
  tCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaa                                                                                       >  2:3237155/1‑93 (MQ=255)
    cAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaa                                                                                    <  1:1757653/93‑1 (MQ=255)
      gCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaag                                                                                  <  1:1641452/93‑3 (MQ=255)
        ccgccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaaaggg                                                                                >  1:191021/1‑89 (MQ=255)
          gccgGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAggtgg                                                                               >  1:4495738/1‑93 (MQ=255)
            cggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTc                                                                             <  1:328812/93‑1 (MQ=255)
            cggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTc                                                                             <  1:4595130/93‑1 (MQ=255)
            cggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTc                                                                             <  2:2863958/93‑1 (MQ=255)
             ggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCg                                                                            <  1:882637/93‑1 (MQ=255)
             ggtcggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCg                                                                            >  2:1487792/1‑93 (MQ=255)
              gtcgttACATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCgg                                                                           <  2:3186016/93‑1 (MQ=255)
               tcggtTCATCCATAAACCCCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGc                                                                          <  2:1816609/93‑1 (MQ=255)
                 ggtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAg                                                                        >  2:1501497/1‑93 (MQ=255)
                  gtTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAGAAAGGTGGTCGGCAgg                                                                       >  1:3202754/1‑93 (MQ=255)
                   ttCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCTGGCAATCTGCAAACGCTGCTGAATCCCGCGGGAAAAGGGGGTCGGCCggcg                                                                     >  1:3913958/1‑91 (MQ=255)
                   ttCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgt                                                                      <  1:4581093/93‑1 (MQ=255)
                   ttCATCCATAAACACCAGCTTCGGATGCGTCACCAGGGTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGGCGGCAGgt                                                                      >  2:1601899/1‑93 (MQ=255)
                    tCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtc                                                                     <  2:1624370/93‑1 (MQ=255)
                    tCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtc                                                                     <  2:2310075/93‑1 (MQ=255)
                     cATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcg                                                                    <  2:3246590/93‑1 (MQ=255)
                     cATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcg                                                                    <  1:2069885/93‑1 (MQ=255)
                      aTCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgt                                                                   >  1:3924557/1‑93 (MQ=255)
                      aTCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgt                                                                   >  2:3448979/1‑93 (MQ=255)
                      aTCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgt                                                                   <  2:1461768/93‑1 (MQ=255)
                       tCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgtc                                                                  <  2:3377322/93‑1 (MQ=255)
                       accATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGgtcgtc                                                                  <  2:3940608/92‑1 (MQ=255)
                        ccATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGtcgtcg                                                                 >  1:2515816/1‑93 (MQ=255)
                        ccATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGtcgtcg                                                                 >  1:1769439/1‑93 (MQ=255)
                         cATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGa                                                                >  2:2122466/1‑93 (MQ=255)
                         cATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGa                                                                <  1:2121277/93‑1 (MQ=255)
                          aTAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGAt                                                               <  1:4019577/93‑1 (MQ=255)
                           tAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATc                                                              >  2:3964713/1‑93 (MQ=255)
                            aaaCACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATcc                                                             >  2:223438/1‑93 (MQ=255)
                            aaaCACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGGGGTCGGCAGGTCGTCGATcc                                                             >  2:3830164/1‑93 (MQ=255)
                                cccAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGtt                                                         <  1:506558/92‑1 (MQ=255)
                                acCAGCTACGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGtt                                                         <  2:3186451/93‑1 (MQ=255)
                                  cAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTTCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTgg                                                       >  1:2912713/1‑93 (MQ=255)
                                  cAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTgg                                                       >  1:1046622/1‑93 (MQ=255)
                                  cAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTgg                                                       >  1:3598966/1‑93 (MQ=255)
                                   aGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGc                                                      >  2:2628282/1‑93 (MQ=255)
                                    gCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCg                                                     >  1:243932/1‑93 (MQ=255)
                                      ttCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCggg                                                   >  1:1210285/1‑93 (MQ=255)
                                      ttCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCggg                                                   >  1:3202855/1‑93 (MQ=255)
                                      ttCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCggg                                                   <  1:2274033/93‑1 (MQ=255)
                                      ttCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCggg                                                   >  2:482944/1‑93 (MQ=255)
                                       tCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGa                                                  <  2:688071/93‑1 (MQ=255)
                                       tCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGa                                                  <  1:2679006/93‑1 (MQ=255)
                                        cGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGaa                                                 <  1:4444220/93‑1 (MQ=255)
                                          gATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtc                                               >  1:2181011/1‑93 (MQ=255)
                                          gATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACAGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtc                                               >  1:1124756/1‑93 (MQ=255)
                                            tGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtctc                                             >  2:3724025/1‑93 (MQ=255)
                                             gCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTcc                                            <  2:1025127/93‑1 (MQ=255)
                                               gTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCAc                                          >  2:2888252/1‑93 (MQ=255)
                                               gTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCAc                                          <  1:4347476/93‑1 (MQ=255)
                                                 cACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACct                                        >  1:1061177/1‑93 (MQ=255)
                                                  aCCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctc                                       >  1:1605258/1‑93 (MQ=255)
                                                  aCCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctc                                       >  1:4746524/1‑93 (MQ=255)
                                                   ccAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctct                                      <  2:1911538/93‑1 (MQ=255)
                                                   ccAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctct                                      <  1:743374/93‑1 (MQ=255)
                                                     aGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTc                                    <  2:932331/93‑1 (MQ=255)
                                                      ggTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTcc                                   >  1:3977388/1‑93 (MQ=255)
                                                        ttGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAg                                 >  1:3131187/1‑93 (MQ=255)
                                                         tGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCAATCCGGTTGGCGGGAATCTCCACCTCTTCCAGc                                >  1:4583355/1‑93 (MQ=255)
                                                          gCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGcc                               <  2:287754/93‑1 (MQ=255)
                                                          gCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGcc                               <  1:979126/93‑1 (MQ=255)
                                                            gggCAATCTGCAAACGCTGCTTCATACCGCCGGAAGAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATTTCCACCTCTTCCCGCCaa                             >  1:851001/1‑92 (MQ=255)
                                                            gggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCAc                             <  1:930716/93‑1 (MQ=255)
                                                            gggCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCAc                             >  1:3803620/1‑93 (MQ=255)
                                                             ggCAATCTGCAAACGCTGCTGCATGCCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACt                            <  1:1208647/93‑1 (MQ=255)
                                                                aaTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCt                         >  2:2650593/1‑93 (MQ=255)
                                                                aaTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCt                         <  1:3953216/93‑1 (MQ=255)
                                                                aaTCTGCAAACGCTGCTCCATCCCGCCGGAAAAAGGTGGGCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCt                         <  1:3953217/93‑1 (MQ=255)
                                                                 aTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTg                        >  2:3295919/1‑93 (MQ=255)
                                                                 aTCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTg                        >  2:1306290/1‑93 (MQ=255)
                                                                  tCTGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCATCCACTTCTgc                       >  1:4581338/1‑93 (MQ=255)
                                                                    tGCAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgcgc                     >  2:1726711/1‑93 (MQ=255)
                                                                      cAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCgg                   <  2:1058197/93‑1 (MQ=255)
                                                                      cAAACGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCgg                   <  2:4565739/93‑1 (MQ=255)
                                                                       aaaCGCTGCTGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGt                  >  1:2808571/1‑93 (MQ=255)
                                                                           gctgctGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGTTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCa              >  1:1722728/1‑93 (MQ=255)
                                                                           gctgctGCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCCCTTCTGCGCCGTggta              >  2:4070643/1‑91 (MQ=255)
                                                                            gtgctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCAc             <  1:2009222/92‑1 (MQ=255)
                                                                            ctgctgCATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCAc             <  1:2009223/93‑1 (MQ=255)
                                                                            ctgctgCATACCGCCGGAAAAAGGGGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCAACCTCTTCCAGCCACTTATGCGCGGTGGCAc             >  2:3619990/1‑93 (MQ=255)
                                                                             tgctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCGCCTCTTCCAGCCACTTCTGCGCGGTGGCACg            >  1:2614681/1‑93 (MQ=255)
                                                                             tgctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACg            >  2:1059583/1‑93 (MQ=255)
                                                                              gctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGa           >  2:1492231/1‑93 (MQ=255)
                                                                               ctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGaa          <  2:4124634/93‑1 (MQ=255)
                                                                               ctgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCCGTGCGACGaa          >  1:3088585/1‑93 (MQ=255)
                                                                                tgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAat         >  2:1226975/1‑93 (MQ=255)
                                                                                tgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAat         >  2:1693585/1‑93 (MQ=255)
                                                                                tgcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAat         >  2:1276199/1‑93 (MQ=255)
                                                                                 gcATCCCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAata        <  1:2695167/93‑1 (MQ=255)
                                                                                 gcATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAata        >  1:4614003/1‑93 (MQ=255)
                                                                                  cATACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAatat       <  1:4629987/93‑1 (MQ=255)
                                                                                   aTACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCGCCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATc      <  1:2102426/93‑1 (MQ=255)
                                                                                   aTACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATc      >  1:4615404/1‑93 (MQ=255)
                                                                                    tACCGCCGGAAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCg     >  2:4729270/1‑93 (MQ=255)
                                                                                      ccgccgGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCg  >  2:3993813/1‑93 (MQ=255)
                                                                                             |                                                                                        
CATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCG  >  NC_000913/4315811‑4315989

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