Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,998 T→A 9.2% K246* (AAG→TAG) ‡ PP_4563 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,9980TA 9.2% 182.9 / 8.4 65K246* (AAG→TAG) ‡PP_4563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base A (3/3);  total (32/33)
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.72e-01

GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACT  >  NC_002947/5182851‑5183146
                                                                                                                                                   |                                                                                                                                                      
gTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTg                                                                                                                                                       >  1:132602/1‑149 (MQ=255)
 tGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGc                                                                                                                                                      <  7:56992/149‑1 (MQ=255)
  ggCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAAGAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCt                                                                                                                                                     >  1:282435/1‑149 (MQ=255)
    cGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                                   <  5:345910/149‑1 (MQ=255)
    cGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                                   <  3:252362/149‑1 (MQ=255)
         agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCaa                                                                                                                                               <  7:57472/148‑1 (MQ=255)
         agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAAc                                                                                                                                              <  2:272370/149‑1 (MQ=255)
          gcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACg                                                                                                                                             >  7:21781/1‑149 (MQ=255)
                cTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATTGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGt                                                                                                                                       <  5:2517/149‑1 (MQ=255)
                   ccGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAgg                                                                                                                                         >  7:169151/1‑144 (MQ=255)
                   ccGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAgg                                                                                                                                         <  8:169151/144‑1 (MQ=255)
                   ccGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCAc                                                                                                                                    >  1:12718/1‑149 (MQ=255)
                     gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActc                                                                                                                                  <  6:332457/149‑1 (MQ=255)
                       aaCTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCtt                                                                                                                                <  2:282435/149‑1 (MQ=255)
                            cGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                                                                                                             <  1:4652/147‑1 (MQ=255)
                                 aGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACAcc                                                                                                                     >  6:17491/1‑149 (MQ=255)
                                 aGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACAcc                                                                                                                     <  8:21781/149‑1 (MQ=255)
                                      tAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCaa                                                                                                                <  7:135729/149‑1 (MQ=255)
                                           ccGTCAGGCTTTCTCGATGCGATCGTCATGGATTACGATCAGGCCCTGCTTGAACAACCAGCCAATGGCCTTCTTGAAGTTGCACTTGCTGACGTTGAACAGCTTGCTGATCAACGAAGGGTCACTCTTGTCGC‑ACAACCCCAacacac                                                                                                           >  3:157906/1‑147 (MQ=255)
                                             gTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCg                                                                                                         >  8:35066/1‑149 (MQ=255)
                                                aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCt                                                                                                      >  7:313002/1‑149 (MQ=255)
                                                   cTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCat                                                                                                                >  7:1188/1‑135 (MQ=255)
                                                   cTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCat                                                                                                                <  8:1188/136‑2 (MQ=255)
                                                          gATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCtt                                                                                               <  8:313002/146‑1 (MQ=255)
                                                          gATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCtt                                                                                               <  5:17491/146‑1 (MQ=255)
                                                           aTGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑Ag                                                                                          <  2:132602/149‑1 (MQ=255)
                                                           aTGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑Ag                                                                                          <  5:45611/149‑1 (MQ=255)
                                                             gCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCACCACCCCGCCTTCGGCTTCC‑AGgc                                                                                        >  8:41317/1‑149 (MQ=255)
                                                               gATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGgcgc                                                                                      <  3:313289/149‑1 (MQ=255)
                                                                   gTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGcc                                                                                   <  5:81872/148‑1 (MQ=255)
                                                                    tCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGc                                                                                                                                                      >  2:70071/1‑82 (MQ=255)
                                                                    tCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGc                                                                                                                                                      <  1:70071/82‑1 (MQ=255)
                                                                    tCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCAt                                                                                 >  7:22412/1‑149 (MQ=255)
                                                                         gATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATc                                                                             <  8:22412/148‑1 (MQ=255)
                                                                         gATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCt                                                                            >  4:35744/1‑149 (MQ=255)
                                                                               gATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                              >  8:227845/1‑79 (MQ=38)
                                                                               gATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                              <  7:227845/79‑1 (MQ=38)
                                                                                       ccTGCTTGAACAACCCGCCAATGGCCTTCTTGATGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGAAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGgtg                                                              >  4:53250/1‑149 (MQ=255)
                                                                                       ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                                 >  4:290153/1‑146 (MQ=255)
                                                                                       ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                                 <  3:290153/146‑1 (MQ=255)
                                                                                       ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGgtg                                                              >  3:163669/1‑149 (MQ=255)
                                                                                              gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGTTCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGGCTTCGGCTTCCAAGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                           <  1:41846/146‑1 (MQ=255)
                                                                                              gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                           >  2:41846/1‑145 (MQ=255)
                                                                                              gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                           >  4:164412/1‑145 (MQ=255)
                                                                                              gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                           <  3:164412/145‑1 (MQ=255)
                                                                                                 caacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAggg                                                    <  7:152127/149‑1 (MQ=255)
                                                                                                     ccGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGc                                                         >  4:37135/1‑140 (MQ=255)
                                                                                                     ccGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGc                                                         <  3:37135/140‑1 (MQ=255)
                                                                                                               tcttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGc                                      >  3:115663/2‑149 (MQ=255)
                                                                                                                 ttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTg                                    >  1:323281/1‑149 (MQ=255)
                                                                                                                 ttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACAACCCGCCTTCGGATTAC‑AGGCTCGCCATGATATGCTCTTGCAGGGTGTCCGCATCGGCTGCGCCTACAGGCt                                     >  1:227391/1‑148 (MQ=255)
                                                                                                                     tGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTAc                                          <  2:54764/139‑1 (MQ=255)
                                                                                                                     tGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTAc                                          >  1:54764/1‑139 (MQ=255)
                                                                                                                         gTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                                 <  2:264821/112‑1 (MQ=255)
                                                                                                                         gTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                                 >  1:264821/1‑112 (MQ=255)
                                                                                                                          ttGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAg                           >  5:278332/1‑149 (MQ=255)
                                                                                                                                   cctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCtt                  >  4:239915/5‑149 (MQ=255)
                                                                                                                                    tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg                 >  7:122502/1‑149 (MQ=255)
                                                                                                                                    ctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAgg                                >  3:210065/4‑134 (MQ=255)
                                                                                                                                    ctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAgg                                <  4:210065/131‑1 (MQ=255)
                                                                                                                                      aCGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGcc               <  3:35744/149‑1 (MQ=255)
                                                                                                                                         ttGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTc            <  2:12718/149‑1 (MQ=255)
                                                                                                                                            aaCAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAgg                                                     >  8:6805/1‑105 (MQ=255)
                                                                                                                                            aaCAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAgg                                                     <  7:6805/105‑1 (MQ=255)
                                                                                                                                              cAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCAACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTATTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTgg        >  1:51125/1‑149 (MQ=255)
                                                                                                                                                   tGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACt  >  6:126397/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACT  >  NC_002947/5182851‑5183146

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