Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,518,736 T→A 8.5% K59* (AAA→TAA) ‡ 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,518,7360TA 8.5% 227.9 / 8.8 82K59* (AAA→TAA) ‡PP_4851phosphate starvation‑inducible protein PsiF
Reads supporting (aligned to +/- strand):  ref base T (41/34);  new base A (2/5);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGT  >  NC_002947/5518588‑5518882
                                                                                                                                                    |                                                                                                                                                  
ggTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCttt                                                                                                                                                    <  4:338106/149‑1 (MQ=255)
 gTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTg                                                                                                                                                   >  8:114597/1‑149 (MQ=255)
          tcttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGCTGCTAATGAACGCTTTGCGTTCATCg                                                                                                                                          >  5:26278/1‑149 (MQ=255)
           cttcAGGCAATCGCTCATGAAGGCCTTGCGCTGATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCATGTCTTCATTTTTTCTTGCTTGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGATCGCAAAGCGTTCAtga                                                                                                                                          <  7:56088/148‑3 (MQ=255)
           cttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAt                                                                                                                                            <  8:162405/146‑1 (MQ=255)
            ttcAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGcc                                                                                                                                        >  1:110248/1‑149 (MQ=255)
               aGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGc                                                                                                                                                  <  6:170460/136‑1 (MQ=255)
               aGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGc                                                                                                                                                  >  5:170460/1‑136 (MQ=255)
               aGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCtt                                                                                                                                     <  5:187931/149‑1 (MQ=255)
                   aaTCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAgg                                                                                                                                 <  2:110248/149‑1 (MQ=255)
                   aaTCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATACTTTTTCAGGCAGGTGCTAATGAAAGCTTTGCGTTCATCGACCTTGAgg                                                                                                                                 >  4:337927/1‑149 (MQ=255)
                   aaTCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTCTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTc                                                                                                                                              <  1:272536/136‑1 (MQ=255)
                   aaTCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTCTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTc                                                                                                                                              >  2:272536/1‑136 (MQ=255)
                       gCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCATTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATAGCCCTTGAgggaat                                                                                                                             >  2:215062/1‑146 (MQ=255)
                               aGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTg                                                                                                                     <  7:114597/149‑1 (MQ=255)
                                   cTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTGTTCATTTTTTCATGCTGGGTTTGGGGAACTTGATTTTTAAGGCAGGTGCTAATGAAGGGTTTGCGGTCATCGCCCTTGAGGTCATTGGCGGTGGCGt                                                                                                                 >  2:91113/1‑149 (MQ=255)
                                    ttGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCGTGCTGGGTTGGGGGCTCATCGTTTATCAGGCCTGTGCTCAAGAACGCGGTGCGTTCATAGCCCTTGAGGGCCTTGGCGGTGGCTTc                                                                                                                >  5:189308/1‑149 (MQ=255)
                                        gcTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCg                                                                                                            >  1:334046/1‑149 (MQ=255)
                                            aTCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGc                                                                                                        >  3:58111/1‑149 (MQ=255)
                                             tCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCa                                                                                                       <  2:334046/149‑1 (MQ=255)
                                               ccccTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGTCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATAGCCCTTGAGGGCGTTGGCGGTGGCGTCGGCGTTGCAgg                                                                                                     >  1:168690/1‑149 (MQ=255)
                                               ccccTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAgg                                                                                                     <  8:47834/149‑1 (MQ=255)
                                                     cAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCa                                                                                               >  2:345079/1‑149 (MQ=255)
                                                       gggCTTTGGTGGTTGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATc                                                                                             <  1:215062/149‑1 (MQ=255)
                                                        ggCTTTGGTGGTGGCGTCAGCATTGCAGGTGTTCATTTTTTCTTGCTGGGTTGGGGGCAAATCCTTTTTCAGGCGGGTGATCATGAACGCTTTGCGTTCATAGGCATTGAGGGCGTTGGAGGTTGAGTCGTCGTTGCTCGTTTTCATCt                                                                                            >  8:246729/1‑149 (MQ=255)
                                                         gCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCt                                                                                                                             <  8:80965/115‑1 (MQ=255)
                                                         gCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCt                                                                                                                             >  7:80965/1‑115 (MQ=255)
                                                          cTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCGTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGAGTCGGCGTTGCAGGTTTTCATCttt                                                                                          >  8:126132/1‑149 (MQ=255)
                                                           tttGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACAACAATTTTCAGGCAGGTGATCATGAACGCTTTGATTTCATCGCCATTGAGGGCATTGGCGGTGGCGTCGGCGTTGCAGGTTTTGATCtttt                                                                                         >  1:145328/1‑149 (MQ=255)
                                                            ttGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTc                                                                                        <  1:91113/149‑1 (MQ=255)
                                                             tggtggtggCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTc                                                                                        <  3:335592/148‑1 (MQ=255)
                                                             tgcgggtggCGTCAGCAGTGCAGGTCTTCATTTTTTCTTGCTGGGTCTCCGGCACATCCTTTTTCAGGCAGCCCCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTc                                                                                        <  6:118183/144‑1 (MQ=255)
                                                              ggtggtggCGTCAGCATTTCATTTGTTCATTTTTTCTTGCTGGTTTTTGCGTACATCCTTTTTCATGCAGGTGCTCATGCACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTcc                                                                                       <  4:58111/148‑1 (MQ=255)
                                                               gtggtggCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCctg                                                                                     >  7:161290/1‑149 (MQ=255)
                                                               ggggtggCGTCAGCATTGCAGGTCTTAATTTTTTCTTTCGGGGCTTGGTGCACATCCTTTTTCAGTAATGTGCTCATGATCTCTTTGCGTTCATCGCCCTTGAGGGCATTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCct                                                                                      <  2:145328/146‑1 (MQ=255)
                                                                 ggtggCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAt                                                                                                                                            <  8:42548/92‑1 (MQ=255)
                                                                 ggtggCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAt                                                                                                                                            >  7:42548/1‑92 (MQ=255)
                                                                    ggCGTCAGCATTGCAGGCCTTCATTTTTTCATGCAGCGTTTGGGGCACATCCTTTTTCAGGGAGGTGCGCATGAACGCTATGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTgcg                                                                                <  7:246729/149‑1 (MQ=255)
                                                                     gCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTgcgc                                                                               >  1:17199/1‑149 (MQ=255)
                                                                        tCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCgg                                                                             >  1:231336/1‑148 (MQ=255)
                                                                        tCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGt                                                                            <  3:337927/149‑1 (MQ=255)
                                                                            cATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTg                                                                        >  2:149881/1‑149 (MQ=255)
                                                                            cATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTg                                                                        >  5:331716/1‑149 (MQ=255)
                                                                             aTTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTgg                                                                       >  8:5523/1‑149 (MQ=255)
                                                                              ttGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCGTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTgg                                                                       >  2:233387/1‑148 (MQ=255)
                                                                                gCAGGACTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGtt                                                                                                   >  4:298265/1‑118 (MQ=255)
                                                                                gCAGGACTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGtt                                                                                                   <  3:298265/118‑1 (MQ=255)
                                                                                 cAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCa                                                                                                       >  4:130190/1‑113 (MQ=255)
                                                                                 cAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCa                                                                                                       <  3:130190/113‑1 (MQ=255)
                                                                                    gTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCa                                                                     >  6:139532/1‑144 (MQ=255)
                                                                                    gTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCa                                                                     <  5:139532/144‑1 (MQ=255)
                                                                                       ttCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAAcc                                                             >  8:259617/1‑149 (MQ=255)
                                                                                       ttCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAt                                                                                                                                            <  5:191058/70‑1 (MQ=255)
                                                                                       ttCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAt                                                                                                                                            >  6:191058/1‑70 (MQ=255)
                                                                                          aTTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTg                                                          >  2:98838/1‑149 (MQ=255)
                                                                                           ttttttCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTgg                                                         >  1:185728/1‑149 (MQ=255)
                                                                                                    gCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGaa                                                 <  1:149881/148‑1 (MQ=255)
                                                                                                    gCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGaa                                                 <  6:331716/148‑1 (MQ=255)
                                                                                                       gggTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTAtt                                                                         <  8:46663/121‑1 (MQ=255)
                                                                                                       gggTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTAtt                                                                         >  7:46663/1‑121 (MQ=255)
                                                                                                       gggTTGGGGTCCCATCCTTCATCACTCTGCTGCTCATGTACGTGTTGCGCTCATCGCCCTTGAGGGCGTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGaacaac                                             <  6:189308/149‑1 (MQ=255)
                                                                                                          tttGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACAcc                                          <  8:161290/149‑1 (MQ=255)
                                                                                                                 cacaTCCGGTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACGGAACAACACcgccagcg                                   >  4:189763/1‑149 (MQ=255)
                                                                                                                    aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgcca                                >  8:63855/1‑149 (MQ=255)
                                                                                                                    aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgc                                  >  1:274297/1‑147 (MQ=255)
                                                                                                                    aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgc                                        >  7:7637/1‑141 (MQ=255)
                                                                                                                    aTCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgc                                        <  8:7637/141‑1 (MQ=255)
                                                                                                                     tCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACcgccagcgccag                               >  5:326433/1‑149 (MQ=255)
                                                                                                                       cTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCAc                                                     <  5:176647/125‑1 (MQ=255)
                                                                                                                       cTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCAc                                                     >  6:176647/1‑125 (MQ=255)
                                                                                                                              aggcatgTGCTCATGTACGCTTTTCGTTTATCGCCCTTGAGGGCCTTGGCGGTGGCGTCTGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGa                       <  2:168690/148‑1 (MQ=255)
                                                                                                                                gcaggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGAttt                    <  7:259617/149‑1 (MQ=255)
                                                                                                                                   ggTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCTAGCGCCAGTAACGGGATTTgca                 <  1:113065/149‑1 (MQ=255)
                                                                                                                                    gTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTgcag                <  3:189763/149‑1 (MQ=255)
                                                                                                                                        cgaTGAACGCAAAGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAAcc                                                             <  8:204298/98‑1 (MQ=255)
                                                                                                                                           tGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACttt          >  6:297466/1‑148 (MQ=255)
                                                                                                                                            gAACGCTTTGTGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTgc                  >  1:95620/1‑139 (MQ=255)
                                                                                                                                            gAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTgc                  <  2:95620/139‑1 (MQ=255)
                                                                                                                                             aaCGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTTCAGGTTTTCATCTTTTCCTGCTGCGCAGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCAt       <  3:301763/149‑1 (MQ=255)
                                                                                                                                             aaCGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCa        <  2:17199/148‑1 (MQ=255)
                                                                                                                                             aaCGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCAt       >  6:220945/1‑149 (MQ=255)
                                                                                                                                                 cTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACt            >  1:138329/1‑140 (MQ=255)
                                                                                                                                                 cTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAg   <  2:231336/149‑1 (MQ=255)
                                                                                                                                                   ttGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACgg                         <  2:291180/125‑1 (MQ=255)
                                                                                                                                                    tGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGt  <  7:5523/147‑1 (MQ=255)
                                                                                                                                                    tGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGt  <  5:297466/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAGT  >  NC_002947/5518588‑5518882

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