Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,024,313 G→C 7.4% S279R (AGC→AGG opdP ← glycine‑glutamate dipeptide porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,024,3130GC 7.4% 577.5 / 27.4 176S279R (AGC→AGGopdPglycine‑glutamate dipeptide porin
Reads supporting (aligned to +/- strand):  ref base G (90/73);  new base C (7/6);  total (97/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 177 total aligned reads displayed.

TCATTGCCGTTGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGTTTCC  >  NC_002947/1024164‑1024460
                                                                                                                                                     |                                                                                                                                                   
      ccGTTGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGc                                                                                                                                               >  1:662004/1‑150 (MQ=255)
           gACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTgc                                                                                                                                          <  1:3333512/150‑1 (MQ=255)
              ctgctgCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCt                                                                                                                                       >  2:390798/1‑150 (MQ=255)
              ctgctgCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCt                                                                                                                                       >  2:4084646/1‑150 (MQ=255)
              ctgctgCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCt                                                                                                                                       <  1:1404386/150‑1 (MQ=255)
                gctgcCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGt                                                                                                                                     >  1:3168690/1‑150 (MQ=255)
                gctgcCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGt                                                                                                                                     >  2:2550622/1‑150 (MQ=255)
                gctgcCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGt                                                                                                                                     <  2:4628739/150‑1 (MQ=255)
                 ctgcCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTc                                                                                                                                    <  1:35055/150‑1 (MQ=255)
                   gcCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGc                                                                                                                                                    >  1:3765744/1‑132 (MQ=255)
                   gcCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGc                                                                                                                                                    <  2:3765744/132‑1 (MQ=255)
                    ccAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGcc                                                                                                                                 >  1:1135867/1‑150 (MQ=255)
                      aggaggCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGCTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTTCGCTGTCTCCGa                                                                                                                               <  2:3295226/150‑1 (MQ=255)
                         aggCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGca                                                                                                                            >  1:1017437/1‑150 (MQ=255)
                          ggCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGcac                                                                                                                           >  2:3627387/1‑150 (MQ=255)
                             gcTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCAGTGCTCAGCCTGAGCAccgg                                                                                                                                          <  1:2497742/132‑4 (MQ=255)
                                  gggTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTgg                                                                                                                   <  2:2122179/150‑1 (MQ=255)
                                    gtgtGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGgcg                                                                                                                 >  2:4958507/1‑150 (MQ=255)
                                     tgtgGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCAGTGCTCAGCCTGAGCAccggc                                                                                                                                         >  2:599040/1‑121 (MQ=255)
                                         ggCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGaa                                                                                                            <  2:5196597/150‑1 (MQ=255)
                                          gCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAag                                                                                                           <  1:4099777/150‑1 (MQ=255)
                                           cTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCCGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGcac                                                                                                                           <  1:5319973/133‑1 (MQ=255)
                                           cTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGcac                                                                                                                           >  2:5319973/1‑133 (MQ=255)
                                            ttGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAagta                                                                                                         <  1:1235119/150‑1 (MQ=255)
                                                     cAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGtt                                                                                                <  2:3451464/150‑1 (MQ=255)
                                                      aGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACCTGGTTGGCGCCGAAGTAGTACTGGTTc                                                                                               >  1:1317074/1‑150 (MQ=255)
                                                      aGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTc                                                                                               <  2:4113095/150‑1 (MQ=255)
                                                      aGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTc                                                                                               >  1:1744936/1‑150 (MQ=255)
                                                        gCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGc                                                                                                                                                    <  1:3309213/95‑1 (MQ=255)
                                                        gCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGc                                                                                                                                                    >  2:3309213/1‑95 (MQ=255)
                                                        gCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCa                                                                                             >  2:4867151/1‑150 (MQ=255)
                                                          ccAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAagta                                                                                                         >  2:1164018/1‑136 (MQ=255)
                                                          ccAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGa                                                                                           >  1:208203/1‑150 (MQ=255)
                                                          ccAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGa                                                                                           >  1:3488977/1‑150 (MQ=255)
                                                           cAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGaa                                                                                          <  2:5150315/150‑1 (MQ=255)
                                                           cAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGaa                                                                                          <  2:5394776/150‑1 (MQ=255)
                                                           cAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGaa                                                                                          <  1:663569/150‑1 (MQ=255)
                                                           cAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGaa                                                                                          <  1:1659966/150‑1 (MQ=255)
                                                            aGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGACTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAg                                                                                         >  2:5334951/1‑150 (MQ=255)
                                                             gTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGt                                                                                        >  2:4986609/1‑150 (MQ=255)
                                                              tGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTc                                                                                       <  2:208203/150‑1 (MQ=255)
                                                                 cAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTc                                                                                    <  2:686621/150‑1 (MQ=255)
                                                                  aGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCg                                                                                   >  2:5182825/1‑150 (MQ=255)
                                                                   ggCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGa                                                                                  <  1:4773299/150‑1 (MQ=255)
                                                                     cTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGAcc                                                                                <  1:3702862/150‑1 (MQ=255)
                                                                     cTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGAcc                                                                                <  2:1744936/150‑1 (MQ=255)
                                                                      tGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCt                                                                               >  1:3862285/1‑150 (MQ=255)
                                                                           gtgtCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTtggtg                                                                          <  1:2758593/150‑1 (MQ=255)
                                                                                 ttGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTAc                                                                    >  1:4980938/1‑150 (MQ=255)
                                                                                      gATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCt                                                               >  1:1369502/1‑150 (MQ=255)
                                                                                       aTCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTg                                                              <  2:2560723/150‑1 (MQ=255)
                                                                                            gCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACt                                                         >  1:3073673/1‑150 (MQ=255)
                                                                                            gCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACt                                                         >  1:3073667/1‑150 (MQ=255)
                                                                                              ccATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCa                                                       <  1:2659784/150‑1 (MQ=255)
                                                                                              ccATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCa                                                       >  1:734181/1‑150 (MQ=255)
                                                                                              ccATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCa                                                       <  1:552581/150‑1 (MQ=255)
                                                                                               cATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAg                                                      >  1:4941126/1‑150 (MQ=255)
                                                                                                aTTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAgg                                                     <  1:4084646/150‑1 (MQ=255)
                                                                                                       ggCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGa                                              >  1:2872033/1‑150 (MQ=255)
                                                                                                          tGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACg                                           >  2:2093831/1‑150 (MQ=255)
                                                                                                            ccGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGc                                         <  2:4203995/150‑1 (MQ=255)
                                                                                                            ccGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGc                                         <  2:859127/150‑1 (MQ=255)
                                                                                                            ccGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGc                                         >  2:4564554/1‑150 (MQ=255)
                                                                                                             cGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGaa                                                                                                            <  2:2822103/82‑1 (MQ=255)
                                                                                                                 gTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTa                                    <  1:167529/150‑1 (MQ=255)
                                                                                                                 gTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTa                                    <  2:3862285/150‑1 (MQ=255)
                                                                                                                 gTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTa                                    <  1:5182825/150‑1 (MQ=255)
                                                                                                                     aCGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTg                                >  1:3633589/1‑150 (MQ=255)
                                                                                                                       ggTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGAt                              <  2:4236671/150‑1 (MQ=255)
                                                                                                                       ggTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGAt                              >  2:1004215/1‑150 (MQ=255)
                                                                                                                         tCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGc                            >  2:1687864/1‑150 (MQ=255)
                                                                                                                         tCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGc                            >  2:2706160/1‑150 (MQ=255)
                                                                                                                         tCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACAGGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGc                            >  2:1405362/1‑150 (MQ=255)
                                                                                                                            tGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCcgg                         >  1:4451100/1‑150 (MQ=255)
                                                                                                                            tGTAGGAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCcgg                         <  2:3084236/150‑1 (MQ=255)
                                                                                                                             gtagtaGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGct                                                   >  2:1121897/1‑123 (MQ=255)
                                                                                                                             gtagtaGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTa                                    >  1:198799/1‑138 (MQ=255)
                                                                                                                             gtagtaGTTCAGGCCGGTGCTCAGGCTGACCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCcggc                        >  1:282377/1‑150 (MQ=255)
                                                                                                                              tagtagTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGaa                                                                                          <  1:4916454/83‑1 (MQ=255)
                                                                                                                               agtagtTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCATACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCcggcgg                      >  2:206826/1‑150 (MQ=255)
                                                                                                                                gtagtTCTGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCAGATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCTGGATCTTGAACGGCTGGTAGTTGATGCCGGCGGt                     >  2:3052049/1‑150 (MQ=255)
                                                                                                                                   gtTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCt                  >  1:3706874/1‑150 (MQ=255)
                                                                                                                                    ttCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCtt                 <  1:4867151/150‑2 (MQ=255)
                                                                                                                                    ttCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTg                 >  1:884792/1‑150 (MQ=255)
                                                                                                                                    ttCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTg                 <  1:1298461/150‑1 (MQ=255)
                                                                                                                                    ttCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGGTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCGTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGGTGATGCCGGCGGTGCTg                 >  2:55594/1‑150 (MQ=255)
                                                                                                                                      cAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGc               <  1:3706985/150‑1 (MQ=255)
                                                                                                                                       aGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTCGTTGATGCCGGCGGTGCTGGCa              >  2:2375686/1‑150 (MQ=255)
                                                                                                                                       aGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCa              >  2:1809240/1‑150 (MQ=255)
                                                                                                                                         cgcagTGCTCAGCCTGAGCACCGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGGAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCAcg            >  1:5174164/5‑150 (MQ=255)
                                                                                                                                          gcagTGCTCAGCCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCGGAAGTAGTACTGGTTCCAGAAGTCCTCGACCGTGGTGGCGTACAGGCTGGGACTCAGGCTCTTGAACGGCTGGTCGTTGATGCCGGCGGTGCTGGCAcgg           >  1:11101/4‑150 (MQ=255)
                                                                                                                                          ccGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCAcgg           >  2:216863/1‑150 (MQ=255)
                                                                                                                                            agTGCTCAGCCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCAcggtc         >  2:4906769/2‑150 (MQ=255)
                                                                                                                                                cTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGtt     >  2:3107594/1‑150 (MQ=255)
                                                                                                                                                 tCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGttt    >  2:1935827/1‑150 (MQ=255)
                                                                                                                                                 tCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGttt    >  1:2735048/1‑150 (MQ=255)
                                                                                                                                                 tCAGCCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGAt                              <  2:2550134/124‑1 (MQ=255)
                                                                                                                                                   aGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGTTTcc  <  1:4224659/150‑1 (MQ=255)
                                                                                                                                                   aGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGTTTcc  >  1:2275527/1‑150 (MQ=255)
                                                                                                                                                   aGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGTTTcc  >  1:2275063/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCATTGCCGTTGACCTGCTGCCAGGAGGCGCTGAGGGTGTGGGCTTGATGGGTCAGGCCCAGTGACAGGCTGTAAGTGTCGTTGTCGATCTTGCCCATTTTCTGGCTGCCGGCGTCGACGGTCTTGTAGTAGTTCAGGCCGGTGCTCAGGCTGAGCACTGCGCTGTCGCCGAGCACATGGTTGGCGCCGAAGTAGTACTGGTTCCAGAAGTCCTCGACCTTGGTGGCGTACAGGCTGGTACTCAGGCTCTTGAACGGCTGGTAGTTGATGCCGGCGGTGCTGGCACGGTCGGTTTCC  >  NC_002947/1024164‑1024460

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