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.

CAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAAC  >  NC_000913/3963154‑3963444
                                                                                                                                                  |                                                                                                                                                  
cAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCt                                                                                                                                                 >  2:504684/1‑149 (MQ=255)
  ggCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGc                                                                                                                                               <  2:553599/149‑1 (MQ=255)
  ggCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGc                                                                                                                                               <  3:245201/149‑1 (MQ=255)
      cTGATGAGCAAATCGCGGCTTATTGCTTCTAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAATTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTACGCTGCgtaa                                                                                                                                           <  1:794559/149‑4 (MQ=255)
      cTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAat                                                                                                                                          <  3:1039099/149‑1 (MQ=255)
         atgaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATc                                                                                                                                       <  3:256695/149‑1 (MQ=255)
           gaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCt                                                                                                                                     >  1:176140/1‑149 (MQ=255)
             gCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGc                                                                                                                                   >  1:406667/1‑149 (MQ=255)
               aaaTCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTc                                                                                                                                 <  3:175427/149‑1 (MQ=255)
                 aTCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTACGCTGCgtaat                                                                                                                                          <  2:898370/139‑5 (MQ=255)
                 aTCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTACGCTGCgtaat                                                                                                                                          >  1:898370/1‑135 (MQ=255)
                  tCGCGGCTTATTGCTTCAAGATGCCATTAGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTTTGCGTAATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGAc                                                                                                                              >  2:787149/1‑149 (MQ=255)
                         ttATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGa                                                                                                                       >  4:326187/1‑149 (MQ=255)
                         ttATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGa                                                                                                                       <  4:126553/149‑1 (MQ=255)
                         ttATTGCTTCAAAATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATTTCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGa                                                                                                                       >  4:913768/1‑149 (MQ=255)
                           aTTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATg                                                                                                                     >  2:1012696/1‑149 (MQ=255)
                            ttGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGc                                                                                                                    >  3:448225/1‑149 (MQ=255)
                             tGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTATGGCTGCGAATATCCTGCTCGACGGCCAGATGCa                                                                                                                   >  2:272995/1‑149 (MQ=255)
                              gCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCaa                                                                                                                  <  2:762988/149‑1 (MQ=255)
                               cTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAAc                                                                                                                 <  1:1075835/149‑1 (MQ=255)
                                ttCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTACGCTGCgt                                                                                                                                             >  3:475615/1‑120 (MQ=255)
                                ttCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTACGCTGCgt                                                                                                                                             <  4:475615/121‑2 (MQ=255)
                                ttCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTCGCAGCCGTACGCTGCGTAatattcaac                                                                                                                                   <  1:377091/130‑6 (MQ=38)
                                ttCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTCGCAGCCGTACGCTGCGTAatattcaac                                                                                                                                   >  2:377091/1‑125 (MQ=38)
                                 tCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCa                                                                                                                         <  2:838932/139‑1 (MQ=255)
                                 tCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCa                                                                                                                         >  1:838932/1‑139 (MQ=255)
                                 tCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATACTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACAt                                                                                                               <  4:43573/149‑1 (MQ=255)
                                 tCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGAGCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACAt                                                                                                               >  3:1080316/1‑149 (MQ=255)
                                  cAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATg                                                                                                              >  2:86834/1‑149 (MQ=255)
                                   aaGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGc                                                                                                             <  2:359864/149‑1 (MQ=255)
                                    aGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGcc                                                                                                            <  3:224321/149‑1 (MQ=255)
                                     gATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATcc                                                                                                                                      <  1:108106/122‑1 (MQ=255)
                                     gATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATcc                                                                                                                                      >  2:108106/1‑122 (MQ=255)
                                       tGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTa                                                                                                         >  3:130284/1‑149 (MQ=255)
                                          cATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGAc                                                                                                      >  1:86891/1‑149 (MQ=255)
                                           aTTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGAcc                                                                                                     >  2:60342/1‑149 (MQ=255)
                                           aTTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGAcc                                                                                                     >  2:820033/1‑149 (MQ=255)
                                           aTTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGAcc                                                                                                     >  2:437495/1‑149 (MQ=255)
                                            ttCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                                    >  3:348987/1‑149 (MQ=255)
                                            ttCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                                    >  1:936207/1‑149 (MQ=255)
                                            ttCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                                    >  3:884786/1‑149 (MQ=255)
                                            ttCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                                    >  1:708572/1‑149 (MQ=255)
                                            atCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGGGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                                    >  2:936150/2‑149 (MQ=255)
                                             tCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAg                                                                                                   <  4:130284/149‑1 (MQ=255)
                                             tCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAg                                                                                                   <  3:893603/149‑1 (MQ=255)
                                               gTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGc                                                                                                 >  2:738359/1‑149 (MQ=255)
                                                      aCCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACg                                                                                          >  3:1000383/1‑149 (MQ=255)
                                                      aCCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACg                                                                                          >  3:505266/1‑149 (MQ=255)
                                                       ccTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGc                                                                                         <  1:1069103/149‑1 (MQ=255)
                                                       ccTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGc                                                                                         >  2:134118/1‑149 (MQ=255)
                                                        cTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCa                                                                                        >  4:441727/1‑149 (MQ=255)
                                                        cTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCa                                                                                        >  1:831421/1‑149 (MQ=255)
                                                              gaagaaGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAc                                                                                  >  2:271107/1‑149 (MQ=255)
                                                              gaagaaGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAc                                                                                  >  3:539177/1‑149 (MQ=255)
                                                               aagaagTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCACGCCTTCACGCAGCGCAcc                                                                                 >  3:336829/1‑149 (MQ=255)
                                                                 gaagTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAccgc                                                                               <  1:1012696/149‑1 (MQ=255)
                                                                             aaCTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCGGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCa                                                                   <  2:182713/149‑1 (MQ=255)
                                                                                   ggCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTTAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACt                                                             <  4:947886/149‑1 (MQ=255)
                                                                                   ggCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACt                                                             <  2:62508/149‑1 (MQ=255)
                                                                                     cTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGa                                                           >  2:1063903/1‑149 (MQ=255)
                                                                                      tACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGaa                                                          >  1:474142/1‑149 (MQ=255)
                                                                                      tACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGaa                                                          <  2:1077136/149‑1 (MQ=255)
                                                                                        cgcgCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAata                                                        <  3:1063534/149‑1 (MQ=255)
                                                                                            cTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCa                                                    <  1:958612/149‑1 (MQ=255)
                                                                                              gtgCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGt                                                  <  3:1723/149‑1 (MQ=255)
                                                                                                gCCTGGTCAATATCGATCATAAAGCGGCGGTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGGCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTc                                                <  2:842715/149‑1 (MQ=255)
                                                                                                gCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTc                                                <  1:504684/149‑1 (MQ=255)
                                                                                                 ccTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCg                                               >  1:494667/1‑149 (MQ=255)
                                                                                                  cTGATCAATCTCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTAATACACTGAATATTCAGTTCgg                                              >  1:947251/1‑149 (MQ=255)
                                                                                                  cTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCgg                                              >  4:881687/1‑149 (MQ=255)
                                                                                                   tGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGt                                             >  2:664156/1‑149 (MQ=255)
                                                                                                      tCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCa                                                                         >  2:878750/1‑118 (MQ=255)
                                                                                                         atatCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAAt                                       >  2:430243/1‑149 (MQ=255)
                                                                                                           atCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATgg                                     <  1:738359/149‑1 (MQ=255)
                                                                                                            tCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                                    <  3:441727/149‑1 (MQ=255)
                                                                                                              gATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGa                                  >  2:520902/1‑149 (MQ=255)
                                                                                                               aTCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                                    <  4:824127/146‑1 (MQ=255)
                                                                                                               aTCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                                    >  3:824127/1‑146 (MQ=255)
                                                                                                               aTCATAAAGCGGCGTTGAATATTACGCAGCGTAAGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGAt                                 >  3:230259/1‑149 (MQ=255)
                                                                                                                 cATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCa                               <  1:437495/149‑1 (MQ=255)
                                                                                                                 cATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACCGCCAGATACAACATGCCGTAAACCAGGCCTTCGCGCAGCGAACCGCCCGCCAGGGTAATACACTGAATATTCAGTTCGGTAAAAATGGCGATCa                               <  2:406667/149‑1 (MQ=255)
                                                                                                                  aTAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACcgcccgcc                                                                          >  2:982963/1‑105 (MQ=255)
                                                                                                                  aTAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                              >  4:167830/1‑149 (MQ=255)
                                                                                                                     aaGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGAt                           <  1:430243/149‑1 (MQ=255)
                                                                                                                          gcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCa                      <  1:60342/149‑1 (MQ=255)
                                                                                                                          gcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCa                      >  2:865175/1‑149 (MQ=255)
                                                                                                                           cgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCat                     <  4:1000383/149‑2 (MQ=255)
                                                                                                                            gTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGa                    <  4:306904/149‑1 (MQ=255)
                                                                                                                             ttGAATATTACGCAGCGTATGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                   >  4:796610/1‑149 (MQ=255)
                                                                                                                             ttGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                   >  2:829543/1‑149 (MQ=255)
                                                                                                                             ttGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                   >  2:919019/1‑149 (MQ=255)
                                                                                                                             ttGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc                   >  3:399684/1‑149 (MQ=255)
                                                                                                                                  tatTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTc              >  1:459658/1‑149 (MQ=255)
                                                                                                                                  tatTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTc              >  3:97806/1‑149 (MQ=255)
                                                                                                                                  tatTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTc              >  3:871131/1‑149 (MQ=255)
                                                                                                                                       cGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGaa         <  2:535250/149‑1 (MQ=255)
                                                                                                                                         cAGCGTACGGCTGCGACTATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATGCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAaca       <  4:879912/149‑1 (MQ=255)
                                                                                                                                               tACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAAc  <  4:693131/149‑1 (MQ=255)
                                                                                                                                               tACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAAc  <  1:85930/149‑1 (MQ=255)
                                                                                                                                               tACAGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAAc  <  2:896635/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAAC  >  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: