Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,948,597 T→G 6.8% D162A (GAC→GCC) ‡ PP_2582 ← heme oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,948,5970TG 6.8% 215.3 / 5.0 71D162A (GAC→GCC) ‡PP_2582heme oxygenase
Reads supporting (aligned to +/- strand):  ref base T (38/28);  new base G (3/2);  total (41/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.48e-01

GTGCAGCACCCGCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCACTGCAA  >  NC_002947/2948457‑2948743
                                                                                                                                            |                                                                                                                                                  
gtgCAGCACCCGCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGcc                                                                                                                                              >  1:127491/1‑147 (MQ=255)
gtgCAGCACCCGCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTg                                                                                                                                            >  1:272349/1‑149 (MQ=255)
               aaGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGcaac                                                                                                                                       <  6:64182/139‑1 (MQ=255)
               aaGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGcaac                                                                                                                                       >  5:64182/1‑139 (MQ=255)
                 ggTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGATGCAATCGCCCCCGCTACCACATCTCCGTCCTCCTGCTCCCCCCGCACCAGGCTGTGCAGCCAGGCAACAAATTcgcgac                                                                                                                            >  4:189498/1‑146 (MQ=255)
                     ccATTGCGCTCGGTCAACTGCCTGCAAATACCTCAAGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATg                                                                                                                        <  3:189498/148‑1 (MQ=255)
                                     aCTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCg                                                                                                       >  7:167671/1‑149 (MQ=255)
                                      ctgcctgcAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCgg                                                                                                      >  6:222445/1‑149 (MQ=255)
                                         cctgcAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGt                                                                                                    >  2:359065/1‑148 (MQ=255)
                                            gcAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCGTCCTCCTGCTCCCCCAGCACCCGGCTGTCCAGCCTGGCAACACCTTCGCGCCAATGCAGCGCTCGGCCATCCGGCTGCgg                                                                                                >  8:184373/1‑149 (MQ=255)
                                              aaaTACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCgg                                                                                                      >  8:175884/1‑141 (MQ=255)
                                              aaaTACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCgg                                                                                                      <  7:175884/141‑1 (MQ=255)
                                              aaaTACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAgggg                                                                                              >  5:1586/1‑149 (MQ=255)
                                                    cTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGt                                                                                        >  3:37276/1‑149 (MQ=255)
                                                     tCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCggg                                                                                                     <  3:327226/135‑1 (MQ=255)
                                                     tCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCggg                                                                                                     >  4:327226/1‑135 (MQ=255)
                                                     tCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGa                                                                                                  <  3:249972/138‑1 (MQ=255)
                                                     tCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGa                                                                                                  >  4:249972/1‑138 (MQ=255)
                                                      cACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGGACAGCCTggt                                                                                                                                          <  6:156630/97‑2 (MQ=255)
                                                      cACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGGACAGCCTggt                                                                                                                                          >  5:156630/1‑96 (MQ=255)
                                                       aCGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGc                                                                                     >  8:287648/1‑149 (MQ=255)
                                                         gCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGt                                                                                        <  5:198292/144‑1 (MQ=255)
                                                         gCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGt                                                                                        >  6:198292/1‑144 (MQ=255)
                                                         gCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCgg                                                                                   <  2:272349/149‑1 (MQ=255)
                                                           aGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGc                                                                                 <  8:167671/149‑1 (MQ=255)
                                                             aTGCACCCGATAACTGGCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCGCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGccc                                                                               <  7:184373/149‑1 (MQ=255)
                                                              tGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCCTCCGGGTGAGGGGCCAGGTGGCGGGccc                                                                               >  7:229790/1‑148 (MQ=255)
                                                              tGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGGACAGCCTggt                                                                                                                                          <  7:173277/89‑2 (MQ=255)
                                                              tGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGGACAGCCTggt                                                                                                                                          >  8:173277/1‑88 (MQ=255)
                                                                 aCCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCAtt                                                                           >  7:186912/1‑149 (MQ=255)
                                                                   ccGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATc                                                                                                         >  5:76198/1‑117 (MQ=255)
                                                                   ccGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATc                                                                                                         <  6:76198/117‑1 (MQ=255)
                                                                    cGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTc                                                                        <  6:1586/149‑1 (MQ=255)
                                                                     gATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCa                                                                       <  2:127491/149‑1 (MQ=255)
                                                                     gATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCAAGCCTCGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCa                                                                       >  8:323947/1‑149 (MQ=255)
                                                                      aTAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCAc                                                                      <  5:313837/149‑1 (MQ=255)
                                                                         aCTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCAccc                                                                    <  5:222445/148‑1 (MQ=255)
                                                                           tGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTcgc                                                                                                                               >  2:65933/1‑87 (MQ=255)
                                                                           tGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTcgc                                                                                                                               <  1:65933/87‑1 (MQ=255)
                                                                             tCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAg                                                               >  6:132996/1‑149 (MQ=255)
                                                                             tCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAg                                                               <  8:186912/149‑1 (MQ=255)
                                                                             tCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAg                                                               >  4:32871/1‑149 (MQ=255)
                                                                             tCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCATGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTTAGGGGCCAGGTGGCGGGCCCCATTGGCACCCTCCAg                                                               >  4:37469/1‑149 (MQ=255)
                                                                                  cGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCaa                                                          <  8:90693/149‑1 (MQ=255)
                                                                                        aaTCGCCCCCGCTACCACCTCTGCTTCCTCCTAGTGCCCCCGCACCAGGCGGGACAGCCTGGCAACCAATTAGCGCCAATGCAGCGCTAGGCCAACCCGGTTAGGGGCCAGGTGGCGGGGCCCATTGTGACCCGGCAGCCCAAGACGtt                                                    >  1:221633/1‑149 (MQ=255)
                                                                                          tCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGc                                                  <  7:287648/149‑1 (MQ=255)
                                                                                            gCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCgt                                                <  7:323947/149‑1 (MQ=255)
                                                                                                ccGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGc                                            >  2:72012/1‑149 (MQ=255)
                                                                                                      ccaccTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTGGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTag                                      <  3:32871/149‑1 (MQ=255)
                                                                                                      ccaccTATGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGGGAGGGGACAGTGGGCGGGCCCCTTTTTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTag                                      <  3:37469/149‑1 (MQ=255)
                                                                                                        accTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTa                                       >  4:144796/1‑146 (MQ=255)
                                                                                                         ccTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCACCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTagag                                    >  1:117663/1‑148 (MQ=255)
                                                                                                            ctGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAAc                                >  6:24231/1‑149 (MQ=255)
                                                                                                                  cctcctGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAAcgcc                             <  3:144796/146‑1 (MQ=255)
                                                                                                                    tcctGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTaga                                     >  5:13871/1‑136 (MQ=255)
                                                                                                                    tcctGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTaga                                     <  6:13871/136‑1 (MQ=255)
                                                                                                                     cctGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGc                                                                                             >  5:254445/1‑79 (MQ=255)
                                                                                                                     cctGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGc                                                                                             <  6:254445/79‑1 (MQ=255)
                                                                                                                          tCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGtt                  >  8:214814/1‑149 (MQ=255)
                                                                                                                            ccccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGtt                  <  1:72012/147‑1 (MQ=255)
                                                                                                                            ccccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTc                 >  8:198591/1‑148 (MQ=255)
                                                                                                                              ccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCTCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGc                                                  <  8:332236/113‑1 (MQ=255)
                                                                                                                              ccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGc                                                  >  7:332236/1‑113 (MQ=255)
                                                                                                                              ccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGaa              >  1:101650/1‑149 (MQ=255)
                                                                                                                              ccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGCTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTAGCGTGTGCTTGTAGAGAAACGCCGCCCCCACGTACGaa              >  3:317895/1‑149 (MQ=255)
                                                                                                                              ccAGCACAAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGaa              >  4:131327/1‑149 (MQ=255)
                                                                                                                                 gCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCt           >  7:146549/1‑149 (MQ=255)
                                                                                                                                   aCCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCtt          <  5:169274/148‑1 (MQ=255)
                                                                                                                                   aCCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCtt          >  6:169274/1‑148 (MQ=255)
                                                                                                                                    ccAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCa        <  7:214814/149‑1 (MQ=255)
                                                                                                                                    ccAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCa        >  6:183745/1‑149 (MQ=255)
                                                                                                                                            tCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCACTGCaa  <  4:37276/147‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GTGCAGCACCCGCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCACTGCAA  >  NC_002947/2948457‑2948743

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