Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,988,829 A→G 19.5% intergenic (+123/‑423) 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,8290AG19.5% 173.8 / 34.3 82intergenic (+123/‑423)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (27/39);  new base G (10/6);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

GTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGG  >  NC_002947/5988682‑5988973
                                                                                                                                                   |                                                                                                                                                      
gttgttTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCAc                                                                                                                                                       >  3:39384/1‑149 (MQ=255)
   gttTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCAcgc                                                                                                                                                     <  1:5040/148‑1 (MQ=255)
   gttTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCAcgc                                                                                                                                                     >  7:100422/1‑148 (MQ=255)
   gttTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCAcgc                                                                                                                                                     >  2:5040/1‑148 (MQ=255)
       gATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCTCGGATAAACCCGCTCCTACGTTTGTTTCGGGCGAGTTATTCCTGTGGGACGCGCACGCGAGcc                                                                                                                                                >  5:313823/1‑149 (MQ=255)
         ttCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGcccc                                                                                                                                              >  3:288002/1‑149 (MQ=255)
            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                                                                                                                     <  7:303704/149‑1 (MQ=255)
                                  gcgcGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGcgtc                                                                                                                     <  2:239346/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                                                                                                                                >  2:248349/1‑118 (MQ=255)
                                                      aGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTc                                                                                                                                <  1:248349/118‑1 (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                                                                                         >  7:130252/1‑139 (MQ=255)
                                                                  aTTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTcgcgcg                                                                                         <  8:130252/139‑1 (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)
                                                                                                                                                   |                                                                                                                                                      
GTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTC‑‑‑‑‑GTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGG  >  NC_002947/5988682‑5988973

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: