Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,889,832 T→C 5.4% L368S (TTG→TCG)  tnaA → tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,889,8320TC 5.4% 221.1 / 3.1 76L368S (TTG→TCG) tnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (48/23);  new base C (2/2);  total (50/26)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

GGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTG  >  NC_000913/3889683‑3889975
                                                                                                                                                     |                                                                                                                                               
ggCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGtt                                                                                                                                                 >  1:496265/1‑150 (MQ=255)
  cGGTAGGTCTGTATTCCGGAATAAAAATCGACTGTCTGGCTTATCGTATCGCGCAGGTACACTATCTTGTCCAAGGGCTGGAAGAGATTGGCGTTGTCTGCCAGAAGGCGGGCGGTCACGCTGCATTCGTTGATGCCGGTACACTGTTGc                                                                                                                                               >  2:62350/1‑150 (MQ=255)
            gTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCCGCAGGCGGGCGGACACCCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCAAATccc                                                                                                                                     >  2:229339/1‑150 (MQ=255)
             tATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTCGCCGCATATCCCg                                                                                                                                    >  2:127693/1‑150 (MQ=255)
             tATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGAATTCGTTGATGCCGGTAAACTGTCGCCGCATATCCCg                                                                                                                                    >  2:131887/1‑150 (MQ=255)
              aTGACGGCCTGAATCTCGACTGGCCGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCgg                                                                                                                                   <  2:75880/150‑1 (MQ=255)
              aTGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCATTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCgg                                                                                                                                   <  2:943640/150‑1 (MQ=255)
               tgccGGCATGAAGCGCGACTGGCTGGCGTATCGTCTCGCGCAGGTACAGGAGCTGGTCGATGGGCTGGAAGAGATGGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGc                                                                                                                                  <  1:190303/147‑1 (MQ=255)
                   ggCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAgcc                                                                                                                              >  1:583269/1‑148 (MQ=255)
                    gCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACATTATCTGGTAGATGGTCTGGAAGAGATTTGCGTTGTCTGACAGCAGGCGTGCGGTCACGCGGCATTCGTTTATGACGGTAAACTGTTACCGCATATCCCGGCATAca                                                                                                                             >  2:982760/1‑149 (MQ=255)
                     cATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCa                                                                                                                            >  2:891371/1‑150 (MQ=255)
                        gAATCTCGACTTGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGCGAGGTCACGCGGCATTCGTTGATGACTGTAAACTGTTGCCGCATATCCAGGCAGACGAGtt                                                                                                                         >  1:106914/1‑150 (MQ=255)
                         aaTCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTc                                                                                                                        >  1:1250447/1‑150 (MQ=255)
                          aTCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGCCCAGTTcc                                                                                                                       >  1:1323261/1‑150 (MQ=255)
                          aTCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTcc                                                                                                                       >  2:956642/1‑150 (MQ=255)
                           tctcGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTccc                                                                                                                      >  1:1072212/1‑150 (MQ=255)
                             tcGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCgg                                                                                                                    >  1:1186338/1‑150 (MQ=255)
                             tcGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCgg                                                                                                                    <  2:1239267/150‑1 (MQ=255)
                                     cttgctTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGACACGCGGCATTCGTTGATGCCGGTAACCTGTTGCCGCATATCCCGGCAGACCAGTTCCAGGCACAGgcg                                                                                                            >  2:1176838/4‑150 (MQ=255)
                                         ctTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTgg                                                                                                        >  1:1023674/1‑150 (MQ=255)
                                            atcgtatcgCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTACCGGCACAGGCGCTGGCCt                                                                                                     >  2:925349/1‑150 (MQ=255)
                                                tatcgCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGa                                                                                                 <  1:84949/150‑1 (MQ=255)
                                                tatcgCGCAGGTACAGAATATGGTCGATGGTCTGGAAGAGATAGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGAGGCATTCGTTGATTCCGGAAAACTGTTGCAGCATATCCCGCAAGACCAGTTCCCGGAGCAGGCGCTGTCCTGCGa                                                                                                 >  1:81929/1‑150 (MQ=255)
                                                  tcgCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGc                                                                                               <  2:1072212/150‑1 (MQ=255)
                                                  tcgCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGc                                                                                               <  1:453780/150‑1 (MQ=255)
                                                    gcgcAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTg                                                                                             <  1:991951/150‑1 (MQ=255)
                                                    gcgcAGGTACAGTATCTGGGCGATGGTCTGGAAGAGATTGGCGGTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTg                                                                                             <  2:786822/150‑1 (MQ=255)
                                                     cgcAGGTCCAGTATCTGGTCGCTGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGt                                                                                            <  1:661432/150‑1 (MQ=255)
                                                      gcAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGta                                                                                           >  1:1159420/1‑150 (MQ=255)
                                                                tATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGcc                                                                                 <  1:1043167/150‑1 (MQ=255)
                                                                   cTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGTCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGCGCTGTATTAAGTCGCCGGt                                                                              >  2:19888/1‑150 (MQ=255)
                                                                    tGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTa                                                                             >  2:905336/1‑150 (MQ=255)
                                                                     ggTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTAt                                                                            >  1:827980/1‑150 (MQ=255)
                                                                         gATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGAGCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAATGTCGCGGGTATCCgg                                                                        >  2:249117/1‑149 (MQ=255)
                                                                         gATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGt                                                                        >  1:1107705/1‑150 (MQ=255)
                                                                         gATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGt                                                                        >  2:716018/1‑150 (MQ=255)
                                                                              tcgggAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGt                                                                   <  2:868020/147‑1 (MQ=255)
                                                                                tGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAg                                                                 >  2:24769/1‑150 (MQ=255)
                                                                                 ggAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGa                                                                >  2:935792/1‑150 (MQ=255)
                                                                                    agagATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAAt                                                             >  2:1005225/1‑150 (MQ=255)
                                                                                    agagATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAAt                                                             >  2:701202/1‑150 (MQ=255)
                                                                                      agaTTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTg                                                           >  2:351008/1‑150 (MQ=255)
                                                                                       gaTTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTgg                                                          <  1:438829/150‑1 (MQ=255)
                                                                                        agtGGAGTTGTCTGCCAGCAGGCTGGCGGTCACGCGGCATTCGTTGTTGGCGGTAAACTGTTGCCGCAAATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGc                                                         <  2:877317/148‑1 (MQ=255)
                                                                                           ggCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGctct                                                      <  2:1272291/150‑1 (MQ=255)
                                                                                           ggCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGcgct                                                      >  2:518482/1‑150 (MQ=255)
                                                                                           ggCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTAGAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGctcg                                                      >  1:511209/1‑149 (MQ=255)
                                                                                            gCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCtt                                                     >  1:2355/1‑150 (MQ=255)
                                                                                              gTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTc                                                   >  1:743404/1‑150 (MQ=255)
                                                                                               ttGGCTGCCAGAAGGCGGTCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTAGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTcc                                                  <  2:836687/150‑1 (MQ=255)
                                                                                                tGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGCGCTGTATAAAGTCGCCGGTATCCGTGCGGGAGAAATTGGCTCTTTCCt                                                 >  2:230108/1‑150 (MQ=255)
                                                                                                     gCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAg                                            >  1:887844/1‑150 (MQ=255)
                                                                                                       cagcagGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCGGGTATCCGTGGCGTAGAAATGGGCTGTTTCCTGTTGGGc                                          >  1:955257/1‑150 (MQ=255)
                                                                                                       cagcagGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGc                                          <  1:1314975/150‑1 (MQ=255)
                                                                                                        agcagGGGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTCGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGcc                                         <  1:131887/150‑1 (MQ=255)
                                                                                                        agcagGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTCGCCGCATATCCCGGCAGTCCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGcc                                         <  1:127693/150‑1 (MQ=255)
                                                                                                          cagGCGGGCGGGCACGCGGCATTCGTTGATGCCGGTAAACTGTTGACGCATATCCCGGCAGACCAGTTCCCGGCTCAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATAGGCTCTTTCCTGTTAGGCcgc                                       <  1:833754/150‑1 (MQ=255)
                                                                                                            ggcgggcgGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGa                                     >  2:1209505/1‑150 (MQ=255)
                                                                                                             gcgggcggTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAATTTGGCTCTTTCCTGTTAGGCCGCGAt                                    >  2:566406/1‑150 (MQ=255)
                                                                                                             gcgggcggTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGAt                                    <  2:364891/150‑1 (MQ=255)
                                                                                                             gcgggcggTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGGATAAAGTCGCCGGGATCCGTGCGGTAGATATTGGCTCTTTCCTGTTAGGCcgcgct                                    >  2:267974/1‑148 (MQ=255)
                                                                                                             gcgggcggTCACGCGGCATTCGGTGATGCCGGTAAACTGGTGCCGCATATCCCGGCAGACCAGTTCCCGGCACTGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGAt                                    <  2:46654/150‑1 (MQ=255)
                                                                                                                      cACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAAcc                           >  2:1134341/1‑150 (MQ=255)
                                                                                                                      cACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAAcc                           <  1:476470/150‑1 (MQ=255)
                                                                                                                       aCGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCg                          >  1:971689/1‑150 (MQ=255)
                                                                                                                                  gTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTg               >  1:39310/1‑150 (MQ=255)
                                                                                                                                  gTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTg               >  1:965557/1‑150 (MQ=255)
                                                                                                                                     gATGCCGGTAAACTGTTGCCGCATATCCCGGCAGCCCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGGGGCCGTGCGGTAGCAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTCACCCACTGCCa            >  1:374480/1‑150 (MQ=255)
                                                                                                                                     gATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCa            >  2:735347/1‑150 (MQ=255)
                                                                                                                                     gATGCCGGTAAACTGTTGCCGCATATCACGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCa            >  2:61936/1‑150 (MQ=255)
                                                                                                                                        gCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGc         >  2:53071/1‑150 (MQ=255)
                                                                                                                                         ccGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGcc        >  2:1181538/1‑150 (MQ=255)
                                                                                                                                         ccGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAAACGGTTAAAACCTGCCATGcc        >  2:1181540/1‑150 (MQ=255)
                                                                                                                                          cGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGccc       <  1:1028364/150‑1 (MQ=255)
                                                                                                                                           ggTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAGACAACTGCCATGCCCg      >  2:1121177/1‑150 (MQ=255)
                                                                                                                                             tAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGc    <  1:891371/150‑1 (MQ=255)
                                                                                                                                               aaCTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGGATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTg  <  1:1176838/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTCGCCGGTATCCGTGCGGTAGAAATTGGCTCTTTCCTGTTAGGCCGCGATCCGAAAACCGGTAAACAACTGCCATGCCCGGCTG  >  NC_000913/3889683‑3889975

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: