Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,506,923 A→G 12.2% intergenic (+44/‑73) petC → / → sspA ubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,506,9230AG12.2% 182.2 / 14.5 74intergenic (+44/‑73)petC/sspAubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base G (4/5);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

AAACTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGG  >  NC_002947/1506775‑1507065
                                                                                                                                                    |                                                                                                                                               
aaaCTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTa                                                                                                                                                 <  5:176453/149‑1 (MQ=255)
      gAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATa                                                                                                                                           >  7:19382/1‑149 (MQ=255)
       aaaGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAt                                                                                                                                          >  6:26303/1‑149 (MQ=255)
       aaaGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAt                                                                                                                                          >  4:286771/1‑149 (MQ=255)
           ccAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGc                                                                                                                                       >  5:287675/1‑148 (MQ=255)
            cAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGACGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCa                                                                                                                                     <  5:26303/149‑1 (MQ=255)
            cAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGcta                                                                                                                                     <  5:54193/149‑3 (MQ=255)
            cAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGcta                                                                                                                                     <  5:48851/149‑3 (MQ=255)
            cAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCa                                                                                                                                     <  3:286771/149‑1 (MQ=255)
              gcgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAgg                                                                                                                                   >  1:266801/1‑149 (MQ=255)
               cgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc                                                                                                                                  <  2:122630/149‑1 (MQ=255)
                gcATCGGCACCTATTTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGCACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGcc                                                                                                                                 <  4:264867/149‑1 (MQ=255)
                gcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGc                                                                                                                                       <  8:28675/143‑1 (MQ=255)
                gcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGc                                                                                                                                       >  7:28675/1‑143 (MQ=255)
                 cATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCa                                                                                                                                <  6:287675/149‑1 (MQ=255)
                  aTCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg                                                                                                                                    <  5:140710/144‑1 (MQ=255)
                  aTCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg                                                                                                                                    >  6:140710/1‑144 (MQ=255)
                   tCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCAcc                                                                                                                              >  2:122779/1‑149 (MQ=255)
                   tCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCAcc                                                                                                                              <  3:308093/149‑1 (MQ=255)
                            aTGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGcgcg                                                                                                                    >  8:223926/1‑149 (MQ=255)
                             tgtgTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGATTACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGt                                                                                                                   >  5:198107/1‑149 (MQ=255)
                                 ttGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTc                                                                                                               <  2:307786/149‑1 (MQ=255)
                                   gCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGAAGTATTGCCAGGCCACCCCGGGCGCGTTTTCAt                                                                                                             >  3:201740/1‑149 (MQ=255)
                                     tGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCAtt                                                                                                            >  5:203381/1‑148 (MQ=255)
                                       tATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCAttt                                                                                                           <  2:35452/147‑1 (MQ=255)
                                       tATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTa                                                                                                          >  7:171257/1‑148 (MQ=255)
                                         tCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTTAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTaaa                                                                                                        <  7:202427/148‑1 (MQ=255)
                                         tCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTaaa                                                                                                        <  2:266801/148‑1 (MQ=255)
                                              atttCTTCTTCGTATTCGCCTACTTGCTCAAGCGTTAATACTGGAAGGACGTGCACTCATCACGCAGTAGTTTCTGGTTGTTCAACGCGCCCGGGGGGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGca                                                                                                  <  3:3620/148‑1 (MQ=255)
                                               tttCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAg                                                                                                                                                <  7:152853/103‑1 (MQ=255)
                                               tttCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAg                                                                                                                                                >  8:152853/1‑103 (MQ=255)
                                               tttCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCCGGCGCGGTTTCATTTAAAGTTGca                                                                                                  >  6:41358/1‑148 (MQ=255)
                                                  cttcttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACaa                                                                                               <  1:112921/148‑1 (MQ=255)
                                                  cttcttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACaa                                                                                               <  3:215588/148‑1 (MQ=255)
                                                    tcttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTATGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCa                                                                                            >  8:227911/1‑149 (MQ=255)
                                                     cttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCCGCGTTTTCATTTAAAGTTGCACaa                                                                                               <  7:223926/146‑1 (MQ=255)
                                                      ttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTAATTTAAAGTTTCACAAGCAta                                                                                          >  8:226159/1‑148 (MQ=255)
                                                           aTTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGc                                                                                     >  7:217961/1‑149 (MQ=255)
                                                               gCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCgag                                                                                  <  6:203381/148‑1 (MQ=255)
                                                                ccTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCATTATTGCCAGGCAACCCCGGGCGCGGTTTCATTTAAAGTTGCACAAGCATACCTGCgatga                                                                                >  6:26242/1‑149 (MQ=255)
                                                                    cTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCg                                                                            <  4:201740/149‑1 (MQ=255)
                                                                         tCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGc                                                                           <  2:292203/145‑1 (MQ=255)
                                                                         tCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGc                                                                           >  1:292203/1‑145 (MQ=255)
                                                                         tCAAGCGTGAATACTCGAAGGACGTGCACTGATCACGAAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGCCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTAATTTAAAGTTGCACAAGCATCCTTACGAGGAGCCGCTACa                                                                       >  7:164926/1‑149 (MQ=255)
                                                                            aGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATgg                                                                    <  6:198107/149‑1 (MQ=255)
                                                                             gCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAg                                                                                                                                                <  7:113225/73‑1 (MQ=255)
                                                                             gCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAg                                                                                                                                                >  8:113225/1‑73 (MQ=255)
                                                                                tGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGc                                                                                     <  2:224267/128‑1 (MQ=255)
                                                                                tGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGc                                                                                     >  1:224267/1‑128 (MQ=255)
                                                                                             gACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTa                                                   <  5:26242/149‑1 (MQ=255)
                                                                                               cGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCCCCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGAGCGACGAGCCGCTACATGGGCGCACCCAACAGGATAGc                                                 >  7:238528/1‑149 (MQ=255)
                                                                                                   cACTGATCACGCATTAGTTGCTGGTTGTTGAACGCGCCCGGGGTGGCCTCGCAGTATGGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTg                                              <  5:41358/148‑1 (MQ=255)
                                                                                                      tGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTAt                                          >  8:135876/1‑149 (MQ=255)
                                                                                                        tgcactgagtagtTTCTGGTTGTTTAAGGGGCCCGGGGTGGCCTAGCAGTTTTGCCAGGCCACCCGGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCGGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTAtt                                         <  8:164926/141‑1 (MQ=255)
                                                                                                        ctCACGCAGTAGTTTCTGGTTGTGGACCGCGCCCGGGGTGGCCTAGCAGTATGGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTGGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTc                                        <  8:238528/148‑1 (MQ=255)
                                                                                                        aTCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCt                                                                          >  8:203600/1‑115 (MQ=255)
                                                                                                        aTCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCt                                                                          <  7:203600/115‑1 (MQ=255)
                                                                                                         tCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTc                                        <  8:171257/148‑1 (MQ=255)
                                                                                                          cACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGc                                                                           <  3:93635/112‑1 (MQ=255)
                                                                                                          cACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGc                                                                           >  4:93635/1‑112 (MQ=255)
                                                                                                          cACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTcc                                       >  4:246666/1‑148 (MQ=255)
                                                                                                            cGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGAt                                    >  2:116295/1‑149 (MQ=255)
                                                                                                                 tagtTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGc                               <  1:122779/149‑1 (MQ=255)
                                                                                                                        tGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCATGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGATGATCCGCTACATGGGCGCAACCAACATGTTATCCTGCTATTCCGATCCCGCTGATCa                         <  7:227911/148‑1 (MQ=39)
                                                                                                                          gttgttGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTa                      <  7:135876/149‑1 (MQ=255)
                                                                                                                             gttgAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTcc                                       >  1:113229/1‑129 (MQ=255)
                                                                                                                             gttgAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTcc                                       <  2:113229/129‑1 (MQ=255)
                                                                                                                               tgAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATtct                  >  1:212136/1‑148 (MQ=255)
                                                                                                                                  aCGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATc              <  7:226159/149‑1 (MQ=255)
                                                                                                                                      gcCCGGGGTGGCCTGGCAATACTGCTAGGCCCCCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGATGATCACTATTATCATCGGGt          >  5:23125/1‑149 (MQ=38)
                                                                                                                                        ccGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTc        >  2:129356/1‑149 (MQ=255)
                                                                                                                                          ggggTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGc      <  3:246666/149‑1 (MQ=255)
                                                                                                                                              tGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTgg  >  6:313129/1‑149 (MQ=255)
                                                                                                                                               ggCCTGGCAATTCTGCTAGGCCACCCCGGGCGCGTTTTGATTTAAAGTTGCACAAGCATCCGTTCGAGGAGCCGCTAAATGGGCGCAACCAACAGGTTAGCCTGCTATTACGATACCGCGGATAACTATTCTAAACGGGttgagatag  >  4:205539/1‑141 (MQ=255)
                                                                                                                                               ggCCTGGCAATTCTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGtt         <  3:205539/141‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AAACTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGG  >  NC_002947/1506775‑1507065

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: