Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,988,844 C→G 22.1% intergenic (+138/‑408) 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,8440CG22.1% 166.3 / 42.9 77intergenic (+138/‑408)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (22/38);  new base G (11/6);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 5.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTAT  >  NC_002947/5988700‑5988988
                                                                                                                                                |                                                                                                                                                      
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                                                                                                         <  5:306244/149‑1 (MQ=255)
                                     gTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTcgcgc                                                                                                         >  1:109951/1‑149 (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                  >  1:190334/1‑149 (MQ=255)
                                                                                                                            cgcgCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATg                  >  8:69166/1‑149 (MQ=255)
                                                                                                                             gcgcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATgg                 >  2:306078/1‑149 (MQ=255)
                                                                                                                             gcgcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATgg                 <  1:306078/149‑1 (MQ=255)
                                                                                                                               gcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACaggag                                                                                         <  2:59869/75‑1 (MQ=255)
                                                                                                                               gcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACaggag                                                                                         >  1:59869/1‑75 (MQ=255)
                                                                                                                                     gaccgccttggggCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaa         <  2:69558‑M2/143‑47 (MQ=255)
                                                                                                                                           cttggGGCATGTCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctc  >  1:69558‑M2/3‑96 (MQ=255)
                                                                                                                                           ccTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaa                                                                            <  8:206819/76‑1 (MQ=255)
                                                                                                                                           ccTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaa                                                                            >  7:206819/1‑76 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTAT  >  NC_002947/5988700‑5988988

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: