Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,988,840 C→T 21.0% intergenic (+134/‑412) kefB‑III → / → PP_mr64 glutathione‑regulated potassium/H+ antiporter/C4 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,988,8400CT21.0% 160.2 / 32.7 76intergenic (+134/‑412)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (22/38);  new base T (10/6);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 8.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAA  >  NC_002947/5988694‑5988981
                                                                                                                                                  |                                                                                                                                                   
ggATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTg                                                                                                                                                   >  8:303704/1‑149 (MQ=255)
ggATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTg                                                                                                                                                   <  4:189042/149‑1 (MQ=255)
      ggACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCAc                                                                                                                                             >  2:194407/1‑149 (MQ=255)
      ggACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCAc                                                                                                                                             <  1:97158/149‑1 (MQ=255)
          tGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCt                                                                                                                                         >  6:19425/1‑149 (MQ=255)
              gCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGa                                                                                                                                      >  8:345142/1‑148 (MQ=255)
               ctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGata                                                                                                                                    <  2:61200‑M2/93‑1 (MQ=255)
                 ttttGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACACGCTCCTACGATTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATc                                                                                                                                  >  4:349281/1‑149 (MQ=255)
                      gtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGcgtc                                                                                                                             <  2:153853‑M2/100‑1 (MQ=255)
                      gcgcGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGcgtc                                                                                                                             <  2:239346/149‑1 (MQ=255)
                      gcgcGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGcgtc                                                                                                                             <  7:303704/149‑1 (MQ=255)
                      gcgcGTGCCGTGAGTGGTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGATAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTTTTCCTGTGGGAGGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGcgtc                                                                                                                             <  6:313823/149‑1 (MQ=37)
                             gcgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTacaa                                                                                                                  <  2:200611‑M2/106‑1 (MQ=255)
                             ccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTacaaa                                                                                                                 >  3:151871‑M2/43‑149 (MQ=255)
                              cTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTacaaa                                                                                                                 >  4:151152/1‑148 (MQ=255)
                               ttGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTacaaac                                                                                                                >  1:271777/1‑148 (MQ=255)
                                  aGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAg                                                                                                            <  7:345142/148‑1 (MQ=255)
                                   gtgtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGc                                                                                                          <  7:98905‑M2/113‑1 (MQ=255)
                                   gggTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGc                                                                                                          <  6:299194/149‑1 (MQ=255)
                                   gggTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGc                                                                                                          <  3:151152/149‑1 (MQ=255)
                                     gTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCgg                                                                                                        <  8:100422/149‑1 (MQ=255)
                                      tttcagcgcctgggagatcgagcgccgcccgcgCTGCGCTTCGCGGGTAACCCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGt                                                                                                       <  5:326934‑M2/116‑1 (MQ=255)
                                          aGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTc                                                                                                                                        <  1:248349/118‑1 (MQ=255)
                                          aGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTc                                                                                                                                        >  2:248349/1‑118 (MQ=255)
                                           gcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGcgtc                                                                                                                             <  3:347552‑M2/100‑1 (MQ=255)
                                           gcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGcgtc                                                                                                                             >  4:347552‑M2/29‑128 (MQ=255)
                                           gTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTcgcgc                                                                                                  <  5:311377/149‑1 (MQ=255)
                                           gTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTcgcgc                                                                                                  >  1:109951/1‑149 (MQ=255)
                                           gTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTcgcgc                                                                                                  <  5:306244/149‑1 (MQ=255)
                                             gCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTAATCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTcgcgcgc                                                                                                <  2:116677/149‑1 (MQ=255)
                                                  tgAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGcc                                                                                             <  2:109951/147‑1 (MQ=255)
                                                  tgAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCt                                                                                            <  4:339551/148‑1 (MQ=255)
                                                   gagaTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTCAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTg                                                                                           <  8:296136/148‑1 (MQ=255)
                                                     gaTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGca                                                                                        <  6:333502/149‑1 (MQ=255)
                                                     gaTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGca                                                                                        >  8:281137/1‑149 (MQ=255)
                                                     cgcgaagcgccgcgcgggCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGca                                                                                        >  3:106075‑M2/19‑149 (MQ=255)
                                                      aTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTcgcgcg                                                                                                 <  8:130252/139‑1 (MQ=255)
                                                      aTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTcgcgcg                                                                                                 >  7:130252/1‑139 (MQ=255)
                                                                 cccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTg                                                                            <  1:194407‑M2/143‑1 (MQ=255)
                                                                  ccgcgCGGCGCTTCGCGGGTAAACAAGCTCCTACGTTTGTTTCGGGCCAGTTATTACTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAAAAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTgg                                                                           >  1:249682‑M2/6‑149 (MQ=255)
                                                                   cgggCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGc                                                                          <  4:106075‑M2/145‑1 (MQ=255)
                                                                   cgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGc                                                                          <  3:191662‑M2/145‑1 (MQ=255)
                                                                          cgcTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaac                                                                   <  2:122283/149‑1 (MQ=255)
                                                                            cTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTgg                                                                                          >  3:6406/1‑124 (MQ=255)
                                                                            cTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTgg                                                                                          <  4:6406/124‑1 (MQ=255)
                                                                             ttCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGAGTGGCCAGaaacaaa                                                                >  8:248139/1‑149 (MQ=255)
                                                                               cgcgGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGGCGCTCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaac                                                               <  7:226569/148‑1 (MQ=255)
                                                                                gcgGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACg                                                              <  4:288002/148‑1 (MQ=255)
                                                                                gcgGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACg                                                              <  3:349281/148‑1 (MQ=255)
                                                                                 cgGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTa                                                            <  2:271777/149‑1 (MQ=255)
                                                                                  gggTACTCCCGCTCCTACGTGTGTTGCGGGCCAGTTATTCCTTTGTGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGGCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAg                                                           <  7:146494/149‑1 (MQ=255)
                                                                                  gggTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAg                                                           <  8:166407/149‑1 (MQ=255)
                                                                                  gggTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGACTggt                                                                                         >  1:296650/1‑118 (MQ=11)
                                                                                    gTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCgg                                                                                                        >  4:175912/1‑102 (MQ=18)
                                                                                        aCCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCgg                                                     <  3:185311/149‑1 (MQ=255)
                                                                                        aCCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCgg                                                     <  7:281137/149‑1 (MQ=255)
                                                                                         cccGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCggg                                                    <  5:19425/149‑1 (MQ=255)
                                                                                          ccGCTCCTAAGTTTGTTTCTGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACTCCTCGATATCGCGTC‑‑‑‑‑GTACAAACCAGGCGGTCGCGCGCGCCTTGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGt                                                   <  2:249682/149‑1 (MQ=255)
                                                                                              tCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCTAAACAAACGTAGGAGCGGGTTTAcc                                              >  4:59253/1‑149 (MQ=11)
                                                                                                cTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGAGGGGGAAAAAGACAGGCGGACGCGCGCGCCTGGCACAGGAGGGATAGGCAAGAAACAAACGTAGGAGCGGGt                                                   >  1:122395/1‑149 (MQ=25)
                                                                                                cTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAccc                                             <  7:145261‑M2/149‑4 (MQ=255)
                                                                                                cTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGGCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAc                                               >  2:154003‑M2/1‑146 (MQ=255)
                                                                                                    gtttgtttCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCcgcg                                          >  7:167978/1‑147 (MQ=11)
                                                                                                    gtttgtttCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCcgcg                                          >  1:219909/1‑147 (MQ=11)
                                                                                                         tttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgc                                    <  2:139377‑M2/149‑13 (MQ=255)
                                                                                                         tttCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGcgcc                                   >  7:149804/1‑149 (MQ=14)
                                                                                                              ggCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCcgcgcg                              >  7:342565/1‑149 (MQ=14)
                                                                                                                ccAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgg                             >  1:139377‑M2/1‑130 (MQ=255)
                                                                                                                   attATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCgcgggcgg                         >  8:324378/2‑149 (MQ=14)
                                                                                                                         ccTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctc                    <  4:191491‑M2/149‑29 (MQ=255)
                                                                                                                                  cgcgCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATg           >  8:69166/1‑149 (MQ=255)
                                                                                                                                  cgcgCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATg           >  1:190334/1‑149 (MQ=255)
                                                                                                                                   gcgcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATgg          <  1:306078/149‑1 (MQ=255)
                                                                                                                                   gcgcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATgg          >  2:306078/1‑149 (MQ=255)
                                                                                                                                     gcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACaggag                                                                                  >  1:59869/1‑75 (MQ=255)
                                                                                                                                     gcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACaggag                                                                                  <  2:59869/75‑1 (MQ=255)
                                                                                                                                           gaccgccttggggCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaa  <  2:69558‑M2/143‑47 (MQ=255)
                                                                                                                                                 ccTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaa                                                                     <  8:206819/76‑1 (MQ=255)
                                                                                                                                                 ccTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaa                                                                     >  7:206819/1‑76 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAA  >  NC_002947/5988694‑5988981

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