Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,963,299 Δ1 bp 5.6% coding (931/1485 nt) gpp ← guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,963,2990G. 5.6% 400.6 / 19.0 124coding (931/1485 nt)gppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (66/50);  new base . (3/4);  total (69/55)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-02

Only 100 of 124 total aligned reads displayed.

CAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAAC  >  NC_000913/3963154‑3963444
                                                                                                                                                  |                                                                                                                                                  
cAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCt                                                                                                                                                 >  2:504684/1‑149 (MQ=255)
  ggCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGc                                                                                                                                               <  2:553599/149‑1 (MQ=255)
  ggCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGc                                                                                                                                               <  3:245201/149‑1 (MQ=255)
      cTGATGAGCAAATCGCGGCTTATTGCTTCTAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAATTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgtaa                                                                                                                                           <  1:794559/149‑4 (MQ=255)
      cTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑Aat                                                                                                                                          <  3:1039099/149‑1 (MQ=255)
         atgaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATc                                                                                                                                       <  3:256695/149‑1 (MQ=255)
           gaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCt                                                                                                                                     >  1:176140/1‑149 (MQ=255)
             gCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGc                                                                                                                                   >  1:406667/1‑149 (MQ=255)
               aaaTCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTc                                                                                                                                 <  3:175427/149‑1 (MQ=255)
                 aTCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgtaat                                                                                                                                          <  2:898370/139‑5 (MQ=255)
                 aTCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgtaat                                                                                                                                          >  1:898370/1‑135 (MQ=255)
                  tCGCGGCTTATTGCTTCAAGATGCCATTAGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTTTGCGTAATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGAc                                                                                                                              >  2:787149/1‑149 (MQ=255)
                         ttATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGa                                                                                                                       >  4:326187/1‑149 (MQ=255)
                         ttATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGa                                                                                                                       <  4:126553/149‑1 (MQ=255)
                         ttATTGCTTCAAAATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATTTCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGa                                                                                                                       >  4:913768/1‑149 (MQ=255)
                           aTTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATg                                                                                                                     >  2:1012696/1‑149 (MQ=255)
                            ttGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGc                                                                                                                    >  3:448225/1‑149 (MQ=255)
                             tGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTATGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCa                                                                                                                   >  2:272995/1‑149 (MQ=255)
                              gCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCaa                                                                                                                  <  2:762988/149‑1 (MQ=255)
                               cTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAAc                                                                                                                 <  1:1075835/149‑1 (MQ=255)
                                ttCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgt                                                                                                                                             >  3:475615/1‑120 (MQ=255)
                                ttCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgt                                                                                                                                             <  4:475615/121‑2 (MQ=255)
                                ttCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATT‑CGCAGCCGTAC‑GCTGCGTAatattcaac                                                                                                                                   <  1:377091/130‑6 (MQ=38)
                                ttCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATT‑CGCAGCCGTAC‑GCTGCGTAatattcaac                                                                                                                                   >  2:377091/1‑125 (MQ=38)
                                 tCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCa                                                                                                                         <  2:838932/139‑1 (MQ=255)
                                 tCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCa                                                                                                                         >  1:838932/1‑139 (MQ=255)
                                 tCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATACTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACAt                                                                                                               <  4:43573/149‑1 (MQ=255)
                                 tCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGAGCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACAt                                                                                                               >  3:1080316/1‑149 (MQ=255)
                                  cAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATg                                                                                                              >  2:86834/1‑149 (MQ=255)
                                   aaGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGc                                                                                                             <  2:359864/149‑1 (MQ=255)
                                    aGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGcc                                                                                                            <  3:224321/149‑1 (MQ=255)
                                     gATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATcc                                                                                                                                      <  1:108106/122‑1 (MQ=255)
                                     gATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATcc                                                                                                                                      >  2:108106/1‑122 (MQ=255)
                                       tGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTa                                                                                                         >  3:130284/1‑149 (MQ=255)
                                          cATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGAc                                                                                                      >  1:86891/1‑149 (MQ=255)
                                           aTTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGAcc                                                                                                     >  2:60342/1‑149 (MQ=255)
                                           aTTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGAcc                                                                                                     >  2:820033/1‑149 (MQ=255)
                                           aTTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGAcc                                                                                                     >  2:437495/1‑149 (MQ=255)
                                            ttCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                                    >  3:348987/1‑149 (MQ=255)
                                            ttCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                                    >  1:936207/1‑149 (MQ=255)
                                            ttCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                                    >  3:884786/1‑149 (MQ=255)
                                            ttCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                                    >  1:708572/1‑149 (MQ=255)
                                            atCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGGGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                                    >  2:936150/2‑149 (MQ=255)
                                             tCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAg                                                                                                   <  4:130284/149‑1 (MQ=255)
                                             tCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAg                                                                                                   <  3:893603/149‑1 (MQ=255)
                                               gTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGc                                                                                                 >  2:738359/1‑149 (MQ=255)
                                                      aCCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACg                                                                                          >  3:1000383/1‑149 (MQ=255)
                                                      aCCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACg                                                                                          >  3:505266/1‑149 (MQ=255)
                                                       ccTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGc                                                                                         <  1:1069103/149‑1 (MQ=255)
                                                       ccTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGc                                                                                         >  2:134118/1‑149 (MQ=255)
                                                        cTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCa                                                                                        >  4:441727/1‑149 (MQ=255)
                                                        cTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCa                                                                                        >  1:831421/1‑149 (MQ=255)
                                                              gaagaaGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAc                                                                                  >  2:271107/1‑149 (MQ=255)
                                                              gaagaaGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAc                                                                                  >  3:539177/1‑149 (MQ=255)
                                                               aagaagTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCACGCCTTCACGCAGCGCAcc                                                                                 >  3:336829/1‑149 (MQ=255)
                                                                 gaagTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAccgc                                                                               <  1:1012696/149‑1 (MQ=255)
                                                                             aaCTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCGGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCa                                                                   <  2:182713/149‑1 (MQ=255)
                                                                                   ggCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTTAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACt                                                             <  4:947886/149‑1 (MQ=255)
                                                                                   ggCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACt                                                             <  2:62508/149‑1 (MQ=255)
                                                                                     cTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGa                                                           >  2:1063903/1‑149 (MQ=255)
                                                                                      tACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGaa                                                          >  1:474142/1‑149 (MQ=255)
                                                                                      tACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGaa                                                          <  2:1077136/149‑1 (MQ=255)
                                                                                        cgcgCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAata                                                        <  3:1063534/149‑1 (MQ=255)
                                                                                            cTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCa                                                    <  1:958612/149‑1 (MQ=255)
                                                                                              gtgCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGt                                                  <  3:1723/149‑1 (MQ=255)
                                                                                                gCCTGGTCAATATCGATCATAAAGCGGCGGTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGGCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTc                                                <  2:842715/149‑1 (MQ=255)
                                                                                                gCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTc                                                <  1:504684/149‑1 (MQ=255)
                                                                                                 ccTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCg                                               >  1:494667/1‑149 (MQ=255)
                                                                                                  cTGATCAATCTCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTAATACACTGAATATTCAGTTCgg                                              >  1:947251/1‑149 (MQ=255)
                                                                                                  cTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCgg                                              >  4:881687/1‑149 (MQ=255)
                                                                                                   tGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGt                                             >  2:664156/1‑149 (MQ=255)
                                                                                                      tCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCa                                                                         >  2:878750/1‑118 (MQ=255)
                                                                                                         atatCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAAt                                       >  2:430243/1‑149 (MQ=255)
                                                                                                           atCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATgg                                     <  1:738359/149‑1 (MQ=255)
                                                                                                            tCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                                    <  3:441727/149‑1 (MQ=255)
                                                                                                              gATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGa                                  >  2:520902/1‑149 (MQ=255)
                                                                                                               aTCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                                    <  4:824127/146‑1 (MQ=255)
                                                                                                               aTCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                                    >  3:824127/1‑146 (MQ=255)
                                                                                                               aTCATAAAGCGGCGTTGAATATTACGCAG‑CGTAAGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGAt                                 >  3:230259/1‑149 (MQ=255)
                                                                                                                 cATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCa                               <  1:437495/149‑1 (MQ=255)
                                                                                                                 cATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACCGCCAGATACAACATGCCGTAAACCAGGCCTTCGCGCAGCGAACCGCCCGCCAGGGTAATACACTGAATATTCAGTTCGGTAAAAATGGCGATCa                               <  2:406667/149‑1 (MQ=255)
                                                                                                                  aTAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACcgcccgcc                                                                          >  2:982963/1‑105 (MQ=255)
                                                                                                                  aTAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                              >  4:167830/1‑149 (MQ=255)
                                                                                                                     aaGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGAt                           <  1:430243/149‑1 (MQ=255)
                                                                                                                          gcgTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCa                      <  1:60342/149‑1 (MQ=255)
                                                                                                                          gcgTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCa                      >  2:865175/1‑149 (MQ=255)
                                                                                                                           cgTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCat                     <  4:1000383/149‑2 (MQ=255)
                                                                                                                            gTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGa                    <  4:306904/149‑1 (MQ=255)
                                                                                                                             ttGAATATTACGCAG‑CGTATGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                   >  4:796610/1‑149 (MQ=255)
                                                                                                                             ttGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                   >  2:829543/1‑149 (MQ=255)
                                                                                                                             ttGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                   >  2:919019/1‑149 (MQ=255)
                                                                                                                             ttGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                   >  3:399684/1‑149 (MQ=255)
                                                                                                                                  tatTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTc              >  1:459658/1‑149 (MQ=255)
                                                                                                                                  tatTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTc              >  3:97806/1‑149 (MQ=255)
                                                                                                                                  tatTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTc              >  3:871131/1‑149 (MQ=255)
                                                                                                                                       cGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGaa         <  2:535250/149‑1 (MQ=255)
                                                                                                                                         cAG‑CGTACGGCTGCG‑ACTATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATGCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAaca       <  4:879912/149‑1 (MQ=255)
                                                                                                                                               tACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAAc  <  4:693131/149‑1 (MQ=255)
                                                                                                                                               tACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAAc  <  1:85930/149‑1 (MQ=255)
                                                                                                                                               tACAGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAAc  <  2:896635/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCG‑AATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAAC  >  NC_000913/3963154‑3963444

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑