Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,518,735 T→A 9.6% 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,518,7350TA 9.6% 218.0 / 9.0 82K59I (AAA→ATA) ‡PP_4851phosphate starvation‑inducible protein PsiF
Reads supporting (aligned to +/- strand):  ref base T (40/33);  new base A (2/6);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAG  >  NC_002947/5518588‑5518881
                                                                                                                                                   |                                                                                                                                                  
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                                                                                                                                                 >  5:170460/1‑136 (MQ=255)
               aGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGc                                                                                                                                                 <  6:170460/136‑1 (MQ=255)
               aGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCtt                                                                                                                                    <  5:187931/149‑1 (MQ=255)
                   aaTCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAgg                                                                                                                                <  2:110248/149‑1 (MQ=255)
                   aaTCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATACTTTTTCAGGCAGGTGCTAATGAAAGCTTTGCGTTCATCGACCTTGAgg                                                                                                                                >  4:337927/1‑149 (MQ=255)
                   aaTCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTCTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTc                                                                                                                                             >  2:272536/1‑136 (MQ=255)
                   aaTCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTCTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTc                                                                                                                                             <  1:272536/136‑1 (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                                                                                                                                           >  7:42548/1‑92 (MQ=255)
                                                                 ggtggCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCAAAGCGTTCAt                                                                                                                                           <  8:42548/92‑1 (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                                                                        >  7:46663/1‑121 (MQ=255)
                                                                                                       gggTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTAtt                                                                        <  8:46663/121‑1 (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)
                                                                                                                                                   |                                                                                                                                                  
GGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGCACATCCTTTTTCAGGCAGGTGCTCATGAACGCTTTGCGTTCATCGCCCTTGAGGGCCTTGGCGGTGGCGTCGGCGTTGCAGGTTTTCATCTTTTCCTGCTGCGCGGTATTGGCAGCGAAACCTTGGGCACTGAACAACACCGCCAGCGCCAGTAACGGGATTTGCAGCACTTTCATGGAG  >  NC_002947/5518588‑5518881

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: