Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,586,290 G→A 6.3% V59V (GTC→GTT)  PP_1390 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,586,2900GA 6.3% 189.0 / 2.2 63V59V (GTC→GTT) PP_1390hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base A (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CGTGGTCGGCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCATGACAAGGGCGACTC  >  NC_002947/1586151‑1586438
                                                                                                                                           |                                                                                                                                                     
cGTGGTCGGCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTCCGTTGAATGAACGgcgccc                                                                                                                                              <  2:237413/149‑3 (MQ=255)
 gTGGTCGGCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCa                                                                                                                                             <  5:296738/149‑1 (MQ=255)
       ggCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCTCTGGGCCAAGGAATACATGAGAGTCGAACAGTGATCGGCAGGGCTCCGTTAATACAACGGACGGCTCTAAAG‑Gtag                                                                                                                                      >  4:214822/1‑147 (MQ=255)
       ggCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTg                                                                                                                                      >  7:162899/1‑149 (MQ=255)
        tcataaCAAGTCGTAAAACCTCAGTGCCACCGAGCTCTAGCACAACTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTTATCGGCAGGGCTCCGTTCATTCAACGGACGTCGCTCAAG‑GTTGc                                                                                                                                     <  4:259259/145‑1 (MQ=255)
        gCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTg                                                                                                                                      >  3:259259/1‑148 (MQ=255)
         cAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCa                                                                                                                                    >  7:276235/1‑149 (MQ=255)
         cAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCa                                                                                                                                    <  4:293200/149‑1 (MQ=255)
           gAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACt                                                                                                                                  <  4:328464/149‑1 (MQ=255)
              cAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACtgtg                                                                                                                               <  5:101748/149‑1 (MQ=255)
                    ggAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCAc                                                                                                                         >  5:64298/1‑149 (MQ=255)
                         aCCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTTGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAAt                                                                                                                    <  7:185980/149‑1 (MQ=255)
                           cTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCAc                                                                                                                         <  5:342407/142‑1 (MQ=255)
                           cTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCAc                                                                                                                         >  6:342407/1‑142 (MQ=255)
                           cTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAAtgt                                                                                                                  >  2:356296/1‑149 (MQ=255)
                              aGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAGTGTTGCACTGTGAGGCACGCAATGTGc                                                                                                                <  1:3953/149‑1 (MQ=255)
                                          ctctCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTa                                                                                                   <  5:108223/149‑1 (MQ=255)
                                           tctcGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTAt                                                                                                  >  3:160968/1‑149 (MQ=255)
                                               gCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGAGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTgg                                                                                               <  3:214822/148‑1 (MQ=255)
                                                       acacCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACgctgg                                                                                      >  4:377654/1‑149 (MQ=255)
                                                         acCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACgctgggc                                                                                    >  8:147824/1‑149 (MQ=255)
                                                             agccagAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGTCGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGggctggc                                                                                <  8:238305/149‑1 (MQ=255)
                                                             agccagAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGggctggc                                                                                >  6:182746/1‑149 (MQ=255)
                                                                cagAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GGTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGc                                                                             <  7:216265/149‑1 (MQ=255)
                                                                 agaagaCTACATTGGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCg                                                                            <  8:276235/149‑1 (MQ=255)
                                                                  gaagaCTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAAtgc                                                                                                                  <  6:140248/110‑2 (MQ=255)
                                                                  gaagaCTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAAtgc                                                                                                                  >  5:140248/1‑109 (MQ=255)
                                                                  gaagaCTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGt                                                                           <  8:162899/149‑1 (MQ=255)
                                                                    agaCTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAg                                                                         <  3:270743/149‑1 (MQ=255)
                                                                      aCTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCa                                                                       <  7:147824/149‑1 (MQ=255)
                                                                       cTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCa                                                                       <  6:64298/148‑1 (MQ=255)
                                                                        tACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGt                                                                     >  7:161985/1‑149 (MQ=255)
                                                                        tACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGt                                                                     >  8:232750/1‑149 (MQ=255)
                                                                        tACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGt                                                                     >  8:337331/1‑149 (MQ=255)
                                                                         aCATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGtt                                                                    >  2:220715/1‑149 (MQ=255)
                                                                          cATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAacgga                                                                                                                                                      >  4:191659/1‑67 (MQ=255)
                                                                          cATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAacgga                                                                                                                                                      <  3:191659/67‑1 (MQ=255)
                                                                           aTTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCg                                                                  <  1:356296/149‑1 (MQ=255)
                                                                              cGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGgc                                                                <  5:182746/148‑1 (MQ=255)
                                                                               gCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGgcgc                                                              >  1:79295/1‑149 (MQ=255)
                                                                                     ccAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCcag                                                         >  3:302441/1‑148 (MQ=255)
                                                                                       aaGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCcagcag                                                      >  2:201588/1‑149 (MQ=255)
                                                                                         ggAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGc                                                    >  6:339743/1‑149 (MQ=255)
                                                                                          gAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCt                                                   >  6:265698/1‑149 (MQ=255)
                                                                                                aTGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                                 >  2:22586/1‑145 (MQ=255)
                                                                                                aTGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                                 <  1:22586/145‑1 (MQ=255)
                                                                                                aTGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                                 <  1:16398/145‑1 (MQ=255)
                                                                                                aTGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                                 >  2:16398/1‑145 (MQ=255)
                                                                                                         gAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCt                                    >  3:328554/1‑149 (MQ=255)
                                                                                                         gAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGACTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCt                                    >  1:48452/1‑149 (MQ=255)
                                                                                                           agagTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGc                                  <  4:160968/149‑1 (MQ=255)
                                                                                                                 aTCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                                 <  6:113000/128‑1 (MQ=255)
                                                                                                                 aTCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGa                                                 >  5:113000/1‑128 (MQ=255)
                                                                                                                    ggCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCgg                          >  4:302739/1‑148 (MQ=255)
                                                                                                                     gCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGt                         >  7:112520/1‑148 (MQ=255)
                                                                                                                      ttgagcgcCGTCCGTTGAATGAACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGt                         >  8:122241/5‑147 (MQ=255)
                                                                                                                       aGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCa                                                                       >  6:248905/1‑100 (MQ=255)
                                                                                                                       aGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCa                                                                       <  5:248905/100‑1 (MQ=255)
                                                                                                                        gggCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAgg                     >  5:1726/1‑149 (MQ=255)
                                                                                                                         ggCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAggg                    <  6:189381/149‑1 (MQ=255)
                                                                                                                         agcgcCGTCCGTTGAATGAACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGt                                                                           <  1:147759/93‑1 (MQ=38)
                                                                                                                         agcgcCGTCCGTTGAATGAACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGt                                                                           >  2:147759/2‑94 (MQ=38)
                                                                                                                             ccGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCAt                <  7:337331/149‑1 (MQ=255)
                                                                                                                             ccGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCAt                >  8:291369/1‑149 (MQ=255)
                                                                                                                               gTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCATGa              <  2:79295/149‑1 (MQ=255)
                                                                                                                                        acggacggCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCATGACAAGGGCGa     <  5:339743/149‑1 (MQ=255)
                                                                                                                                           gacggCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCATGACAAGGGCGActc  >  2:266669/1‑149 (MQ=255)
                                                                                                                                           gacggCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCATGACAAGGGCGActc  >  6:190381/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
CGTGGTCGGCAGAACAAGTCGGAAAACCTCAGTGCCACCGAGCTCTCGCACACCTACACCCAGCCAGAAGACTACATTCGCTGGGCCAAGGAATACATGAGCGTCGAAGAGTGATCGGCAGGGCGCCGTTCATTCAACGGACGGCGCTCAAG‑GTTGCACTGTGAGGCACGCAATGTGCCTGCTCGCAGCTATTGGCGACGCTGGGCTGGCTGCGTAGCAGTTCGGCGCGCCAGCAGGCTGAGCCATACAGCCCTGCGGCAGCGGTAAGGGCCATGACAAGGGCGACTC  >  NC_002947/1586151‑1586438

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: ‑