Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,520,705 T→A 5.0% K59I (AAA→ATA) ‡ PP_4851 ← phosphate starvation‑inducible protein PsiF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,520,7050TA 5.0% 757.3 / 19.4 218K59I (AAA→ATA) ‡PP_4851phosphate starvation‑inducible protein PsiF
Reads supporting (aligned to +/- strand):  ref base T (118/89);  new base A (7/4);  total (125/93)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 230 total aligned reads displayed.

AGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTG  >  NC_002947/5520557‑5520853
                                                                                                                                                    |                                                                                                                                                    
 ggTCATTTCTTCTTCAGGCAATCGCTCATGAAGGACTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTg                                                                                                                                                    <  2:928227/150‑1 (MQ=255)
  gTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGc                                                                                                                                                   <  1:1010330/150‑1 (MQ=255)
  gTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGc                                                                                                                                                   <  2:3394899/150‑1 (MQ=255)
  gTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGc                                                                                                                                                   <  2:343900/150‑1 (MQ=255)
   tCATTTCTTCTACAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCg                                                                                                                                                  <  2:3913879/150‑1 (MQ=255)
    cATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGt                                                                                                                                                 >  1:3894970/1‑150 (MQ=255)
     aTTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGtt                                                                                                                                                >  2:33809/1‑150 (MQ=255)
       ttcttcttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCa                                                                                                                                              <  1:4029164/150‑1 (MQ=255)
       ttcttcttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCa                                                                                                                                              <  2:2070050/150‑1 (MQ=255)
        tcttcttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAt                                                                                                                                             <  1:1507246/150‑1 (MQ=255)
            cttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGcc                                                                                                                                         <  2:1816653/150‑1 (MQ=255)
                 ggCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGa                                                                                                                                    <  2:3974036/150‑1 (MQ=255)
                  gCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAg                                                                                                                                   >  2:1271361/1‑150 (MQ=255)
                   cAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCGTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAgg                                                                                                                                  <  2:2781826/150‑1 (MQ=255)
                     aTCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGc                                                                                                                                <  1:2460572/150‑1 (MQ=255)
                      tCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGcc                                                                                                                               >  1:455062/1‑150 (MQ=255)
                              gAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGt                                                                                                                       >  1:4237913/1‑150 (MQ=255)
                               aaGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCAGTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTg                                                                                                                      >  1:792870/1‑150 (MQ=255)
                                     ttGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCg                                                                                                                >  1:3119641/1‑150 (MQ=255)
                                     ttGCGCTCATCCCCCGTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCg                                                                                                                >  2:14029/1‑150 (MQ=255)
                                      tGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCgg                                                                                                               <  2:4167063/150‑1 (MQ=255)
                                           tCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTg                                                                                                          >  2:3000793/1‑150 (MQ=255)
                                             aTCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCa                                                                                                        <  2:2039889/150‑1 (MQ=255)
                                              tCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAg                                                                                                       >  2:2576444/1‑150 (MQ=255)
                                                ccccTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGc                                                                                                                    >  1:3160108/1‑135 (MQ=255)
                                                 cccTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGtt                                                                                                    >  1:2158381/1‑150 (MQ=255)
                                                    ttCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTc                                                                                                 >  2:3330402/1‑150 (MQ=255)
                                                      cAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTACTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCAt                                                                                               <  2:2870400/150‑1 (MQ=255)
                                                        gggCTTTGGTGGTGGCGTCAGCATTGCTGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCt                                                                                             <  2:229836/150‑1 (MQ=255)
                                                         ggCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCtt                                                                                            >  1:4080102/1‑150 (MQ=255)
                                                            tttGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTc                                                                                         <  2:1112854/150‑1 (MQ=255)
                                                            tttGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTc                                                                                         <  1:955392/150‑1 (MQ=255)
                                                            tttGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTc                                                                                         >  1:520258/1‑150 (MQ=255)
                                                                gtggtggCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCctgc                                                                                     >  1:3210409/1‑150 (MQ=255)
                                                                 tggtggCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCctgct                                                                                    >  1:1284634/1‑150 (MQ=255)
                                                                        gTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGt                                                                             >  1:1889010/1‑150 (MQ=255)
                                                                         tCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTa                                                                            >  1:2061560/1‑150 (MQ=255)
                                                                           aGCCTTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTAtt                                                                          <  1:3435240/150‑1 (MQ=255)
                                                                             cATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTgg                                                                        <  1:3268053/150‑1 (MQ=255)
                                                                             cATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTgg                                                                        <  1:3068189/150‑1 (MQ=255)
                                                                              aTTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGc                                                                       <  1:1437774/150‑1 (MQ=255)
                                                                                  cAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTAGTGGCAGCg                                                                   >  2:4249156/1‑150 (MQ=255)
                                                                                   aGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGa                                                                  <  1:2580557/150‑1 (MQ=255)
                                                                                      tCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGTTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCATGCTGCGCGGTATTGGCAGCGAAAc                                                               >  1:4516279/1‑150 (MQ=255)
                                                                                      tCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAAc                                                               >  2:387476/1‑150 (MQ=255)
                                                                                        ttCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCt                                                             >  1:2391652/1‑150 (MQ=255)
                                                                                          cATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTg                                                           <  2:3136935/150‑1 (MQ=255)
                                                                                           aTTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGTAACCTTgg                                                          <  1:3623450/150‑1 (MQ=255)
                                                                                           aTTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTgg                                                          <  1:3475256/150‑1 (MQ=255)
                                                                                            ttttttCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTggg                                                         >  2:1781038/1‑150 (MQ=255)
                                                                                             tttttCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGc                                                        >  2:3639920/1‑150 (MQ=255)
                                                                                             tttttCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGc                                                        >  1:1782164/1‑150 (MQ=255)
                                                                                                ttCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACt                                                     <  1:185289/150‑1 (MQ=255)
                                                                                                ttCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACt                                                     >  1:3403810/1‑150 (MQ=255)
                                                                                                   ttGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGaa                                                  <  2:56547/150‑1 (MQ=255)
                                                                                                   ttGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGaa                                                  >  1:2227453/1‑150 (MQ=255)
                                                                                                    tGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGaac                                                 <  1:683789/150‑1 (MQ=255)
                                                                                                    tGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGaac                                                 <  2:259345/150‑1 (MQ=255)
                                                                                                       tGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGaacaac                                              >  2:695304/1‑150 (MQ=255)
                                                                                                          gTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACAcc                                           >  2:1428653/1‑150 (MQ=255)
                                                                                                          gTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACAcc                                           >  1:1684763/1‑150 (MQ=255)
                                                                                                          gTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACAcc                                           <  1:118627/150‑1 (MQ=255)
                                                                                                               gggCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccag                                      >  2:1488875/1‑150 (MQ=255)
                                                                                                                    caTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgcca                                 <  1:4295936/150‑1 (MQ=255)
                                                                                                                     aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgccag                                >  2:348681/1‑150 (MQ=255)
                                                                                                                     aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgccag                                >  1:2622576/1‑150 (MQ=255)
                                                                                                                     aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgccag                                >  1:624204/1‑150 (MQ=255)
                                                                                                                      tCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGt                               >  2:1577153/1‑150 (MQ=255)
                                                                                                                      tCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGt                               >  2:2199481/1‑150 (MQ=255)
                                                                                                                      tCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGt                               >  1:3348058/1‑150 (MQ=255)
                                                                                                                         tttttCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAAc                            >  1:1084361/1‑150 (MQ=255)
                                                                                                                         tttttCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCAATGAACAACACCGCCAGCGCCAGTAAc                            <  1:84559/150‑1 (MQ=255)
                                                                                                                                  caggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTgc                   <  1:301254/150‑1 (MQ=255)
                                                                                                                                    ggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTgcag                 >  2:108407/1‑150 (MQ=255)
                                                                                                                                     gTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGTGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTgcagc                <  1:2800896/150‑1 (MQ=255)
                                                                                                                                        cTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACt             >  2:3593692/1‑150 (MQ=255)
                                                                                                                                         tCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACtt            >  2:1297586/1‑150 (MQ=255)
                                                                                                                                         cgaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACtt            >  2:1925675/3‑150 (MQ=255)
                                                                                                                                          gaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTgc                   >  2:2820748/2‑142 (MQ=255)
                                                                                                                                          gaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACttt           >  2:3209573/2‑150 (MQ=255)
                                                                                                                                          gaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACttt           >  1:82489/2‑150 (MQ=255)
                                                                                                                                           aTGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTCTTGGCAGCGAAACCGTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGAACTTTc          >  1:3022419/1‑150 (MQ=255)
                                                                                                                                           aTGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTc          >  2:1707051/1‑150 (MQ=255)
                                                                                                                                            tGAACGCTTTGCGTTCATCGCCCTTGAGGTCCTTGGCGGTGGCGTCGGCGTTGCAGGTATTCATCTTTTACTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTACCGGCATTTGCAGCACTATCa         >  2:4168514/1‑150 (MQ=255)
                                                                                                                                            tGAACGCAAAGCGTTCTTCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCa         >  2:3462570/1‑150 (MQ=255)
                                                                                                                                             gAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCAt        >  2:875157/1‑150 (MQ=255)
                                                                                                                                              aaCGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATg       >  2:262393/1‑150 (MQ=255)
                                                                                                                                              aaCGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATg       >  1:3528766/1‑150 (MQ=255)
                                                                                                                                              aaCGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATg       <  2:1211243/150‑1 (MQ=255)
                                                                                                                                              aaCGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATg       >  2:2007082/1‑150 (MQ=255)
                                                                                                                                                cGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGa     <  2:3754387/150‑1 (MQ=255)
                                                                                                                                                cGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGa     <  2:4237913/150‑1 (MQ=255)
                                                                                                                                                cGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGa     <  2:3241640/150‑1 (MQ=255)
                                                                                                                                                   tttGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgcca                                       <  2:2227185/113‑1 (MQ=255)
                                                                                                                                                   tttGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgcca                                       >  1:2227185/1‑113 (MQ=255)
                                                                                                                                                   tttGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTg  <  2:346989/150‑1 (MQ=255)
                                                                                                                                                   tttGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTg  <  2:455041/150‑1 (MQ=255)
                                                                                                                                                   tttGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTg  <  2:455062/150‑1 (MQ=255)
                                                                                                                                                   tttGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTg  >  1:2513827/1‑150 (MQ=255)
                                                                                                                                                   tttGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGGAACGGGATTTGCAGCACTTTCATGGAGTg  >  1:355467/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGTG  >  NC_002947/5520557‑5520853

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: ‑