Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,004,526 G→T 6.1% G152V (GGC→GTC) ‡ PP_5259 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,004,5260GT 6.1% 210.8 / 2.4 65G152V (GGC→GTC) ‡PP_5259LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base T (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CGAGTTTTGGCGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCGGAG  >  NC_002947/6004382‑6004671
                                                                                                                                                 |                                                                                                                                                 
cGAGTTTTGGCGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCAcc                                                                                                                                               >  6:52549/1‑149 (MQ=255)
   gTTATGGCGTGCCCAGGATGATCCCGCACCTCAAGGGGTATGGAAAACGGCACGGCAATATCCACCTGGACATTCGAAACGAGATCGAGCCGTTCGACCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTAAGGCCAGGGCACCTgg                                                                                                                                            >  8:17626/1‑149 (MQ=255)
       tGGCGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGTGCCCTGGcca                                                                                                                                         <  5:140144/148‑2 (MQ=255)
         gCGTGCGGTTGTTTATCCAGCACCTCAAGGGGGTTGGGAAACGGCAGGGCAATATCAAGCTGGAGATTCGCAACGATATGGAGCCGTTCGCCCTGCTGCAAGGCTGGGCTGACGTGTTGTTCTTCTTTGGCCATGGCACCTGGCCggg                                                                                                                                       <  7:17626/148‑1 (MQ=255)
         gCGTGCGCTGGTTGATCCCGCACCTCAGGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCgggg                                                                                                                                      <  1:305573/149‑1 (MQ=255)
         gCGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCgggg                                                                                                                                      <  5:52549/149‑1 (MQ=255)
          cGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCggggg                                                                                                                                     <  3:74385/149‑1 (MQ=255)
             gcgcTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGCCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCa                                                                                                                                  <  2:154458/149‑1 (MQ=255)
              cgcTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCAc                                                                                                                                 <  7:106058/149‑1 (MQ=255)
                    ttGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGt                                                                                                                           <  8:70851/149‑1 (MQ=255)
                      gATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTg                                                                                                                          >  6:108819/1‑148 (MQ=255)
                          ccGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAgtt                                                                                                                     >  2:22949/1‑149 (MQ=255)
                            gCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAgttgt                                                                                                                   <  1:22949/149‑1 (MQ=255)
                               ccTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTc                                                                                                                <  8:150503/149‑1 (MQ=255)
                               ccTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGAGTGGAGTTGTTTc                                                                                                                >  1:265505/1‑149 (MQ=255)
                                   aaGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTcgc                                                                                                              >  5:326914/1‑147 (MQ=255)
                                   aaGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCga                                                                                                            <  1:230128/149‑1 (MQ=255)
                                   aaGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCga                                                                                                            <  2:276102/149‑1 (MQ=255)
                                           tGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAggtg                                                                                                     <  6:80797/148‑1 (MQ=255)
                                              cAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCGCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGc                                                                                                 <  4:103905/149‑1 (MQ=255)
                                                aaCGGCACGGCAATATCCACCTGGATAATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGa                                                                                                                        >  1:164185/1‑125 (MQ=255)
                                                      acggcaATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGc                                                                                                                                                  <  7:120007/92‑1 (MQ=255)
                                                      acggcaATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGc                                                                                                                                                  >  8:120007/1‑92 (MQ=255)
                                                        ggcaATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTgc                                                                                       <  8:144163/149‑1 (MQ=255)
                                                        ggcaATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTgc                                                                                       <  5:210237/149‑1 (MQ=255)
                                                                          cATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCg                                                                      >  3:61892/1‑149 (MQ=255)
                                                                           aTTCGCATCGAGATGGAGCCGTTCGCCCTGCTGCAAGGCGCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCgg                                                                     <  5:183134/149‑1 (MQ=255)
                                                                               gCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGgcg                                                                  >  4:307733/1‑148 (MQ=255)
                                                                                     agaTGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTg                                                           >  8:119232/1‑149 (MQ=255)
                                                                                        tGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGcc                                                        >  2:57759/1‑149 (MQ=255)
                                                                                         ggAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGgt                                                                                             >  3:160470/1‑111 (MQ=255)
                                                                                         ggAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGgt                                                                                             <  4:160470/111‑1 (MQ=255)
                                                                                         ggAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCg                                                       >  8:210938/1‑149 (MQ=255)
                                                                                             ccGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGc                                                   <  6:326914/149‑1 (MQ=255)
                                                                                                     ccTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGc                                           >  7:289906/1‑149 (MQ=255)
                                                                                                      ctgctgCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCt                                          <  7:119232/149‑1 (MQ=255)
                                                                                                        gctgcAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTgg                                        >  8:96545/1‑149 (MQ=255)
                                                                                                                gCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGccg                                <  6:241774/149‑3 (MQ=255)
                                                                                                                  tCGGCTTACTTGGGGTTCTTCTATGGCCTGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGTTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGt                              <  2:109145/149‑1 (MQ=255)
                                                                                                                  tCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGcc                                                                       <  1:299307/108‑1 (MQ=255)
                                                                                                                  tCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGcc                                                                       >  2:299307/1‑108 (MQ=255)
                                                                                                                  tCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGcc                                                        >  2:244004/1‑123 (MQ=255)
                                                                                                                    ggCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGc                                  >  1:170201/1‑143 (MQ=255)
                                                                                                                    ggCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGc                                  <  2:170201/143‑1 (MQ=255)
                                                                                                                    ggCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGtgc                            >  4:148742/1‑149 (MQ=255)
                                                                                                                      cTGCCGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTGGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGtgctg                          >  6:229290/1‑149 (MQ=255)
                                                                                                                         aCGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCgagga                                             >  5:293657/1‑127 (MQ=255)
                                                                                                                         aCGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCgagga                                             <  6:293657/127‑1 (MQ=255)
                                                                                                                            tggtgTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAgctg                                                                             <  1:124344/92‑1 (MQ=255)
                                                                                                                            tggtgTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAgctg                                                                             >  2:124344/1‑92 (MQ=255)
                                                                                                                              gtgtTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCaaa                   <  3:148742/148‑1 (MQ=255)
                                                                                                                                   cttctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCAcc              <  3:307733/148‑1 (MQ=255)
                                                                                                                                   cttctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCt             >  3:208803/1‑149 (MQ=255)
                                                                                                                                     tctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTcg           >  5:313701/1‑149 (MQ=255)
                                                                                                                                     tctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTcg           >  5:145170/1‑149 (MQ=255)
                                                                                                                                      ctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTcgc          >  3:239861/1‑149 (MQ=255)
                                                                                                                                      ctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTc            <  7:210938/147‑1 (MQ=255)
                                                                                                                                      ccccggCCAGGTGCCCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTg                                                           <  2:24900/96‑1 (MQ=255)
                                                                                                                                      ccccggCCAGGTGCCCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTg                                                           >  1:24900/5‑100 (MQ=255)
                                                                                                                                        ccggCCAGGTGCCCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCgcgc        >  5:101079/3‑149 (MQ=255)
                                                                                                                                           gCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCg     >  1:227194/1‑149 (MQ=255)
                                                                                                                                             cAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCg     >  7:180031/1‑147 (MQ=255)
                                                                                                                                             cAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCg     >  3:198651/1‑147 (MQ=255)
                                                                                                                                             cAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCGGa   <  4:61892/149‑1 (MQ=255)
                                                                                                                                              aGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCGGAg  >  7:212964/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGAGTTTTGGCGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGACATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCGGAG  >  NC_002947/6004382‑6004671

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