Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,150,078 A→T 5.8% intergenic (+33/‑364) PP_3654 → / → PP_3655 leucine‑responsive regulatory protein/cytosine/purine/uracil/thiamine/allantoin permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,150,0780AT 5.8% 233.7 / 4.6 85intergenic (+33/‑364)PP_3654/PP_3655leucine‑responsive regulatory protein/cytosine/purine/uracil/thiamine/allantoin permease family protein
Reads supporting (aligned to +/- strand):  ref base A (37/42);  new base T (2/3);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

CTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCT  >  NC_002947/4149929‑4150227
                                                                                                                                                       |                                                                                                                                                     
cTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTTCCTATGTGAACCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTA                                                                                                                                                       >  2:1263913‑M1/1‑150 (MQ=255)
   ccGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAC                                                                                                                                                    >  1:247803‑M1/1‑150 (MQ=255)
   ccGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAC                                                                                                                                                    <  1:1874613‑M1/150‑1 (MQ=255)
       cGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCATGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGC                                                                                                                                                <  1:707145‑M1/150‑1 (MQ=255)
          ttCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaa                                                                                                                                             <  1:2307048/150‑3 (MQ=255)
              gAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCa                                                                                                                                         >  2:1567533/1‑150 (MQ=255)
               aaGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCaa                                                                                                                                        >  2:2108885/1‑150 (MQ=255)
                 gCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAg                                                                                                                                      <  2:1300523/150‑1 (MQ=255)
                      aTGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCCGCCGCCGCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGAc                                                                                                                                 <  1:436119/150‑1 (MQ=255)
                       tGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACg                                                                                                                                <  1:170820/150‑1 (MQ=255)
                           gAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCCCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGc                                                                                                                            >  2:121141/1‑150 (MQ=255)
                                ccTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcaccttt                                                                                                                       >  1:693575‑M1/1‑125 (MQ=255)
                                ccTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaagg                                                                                                                                           >  1:139966/1‑126 (MQ=255)
                                ccTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaagg                                                                                                                                           <  2:139966/130‑5 (MQ=255)
                                ccTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTc                                                                                                                       <  1:2108885/150‑1 (MQ=255)
                                 cTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCg                                                                                                                      >  2:319632/1‑150 (MQ=255)
                                  tGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgc                                                                                                                     <  2:693575‑M1/150‑28 (MQ=255)
                                    aCCAGTATTCCTAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAg                                                                                                                   <  1:1165362/150‑1 (MQ=255)
                                    aCCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAg                                                                                                                   <  2:1491883/150‑1 (MQ=255)
                                    aCCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCCGCAAAGGTGACGTAGCCCGgcgccg                                                                                                                   >  2:1741210/1‑148 (MQ=255)
                                      cAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAg                                                                                                                                      <  2:863763/129‑1 (MQ=255)
                                      cAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAg                                                                                                                                      >  1:863763/1‑129 (MQ=255)
                                      cAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGa                                                                                                                 >  2:419410/1‑150 (MQ=255)
                                       aGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaa                                                                                                                                             <  2:65516/121‑3 (MQ=255)
                                       aGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaa                                                                                                                                             >  1:65516/1‑119 (MQ=255)
                                         tATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcg                                                                                                              >  1:2021610‑M1/1‑116 (MQ=255)
                                         tATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATAGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgc                                                                                                               >  1:733333‑M1/1‑116 (MQ=255)
                                         tATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATAGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgc                                                                                                               <  2:733333‑M1/149‑34 (MQ=255)
                                          aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcccctttgcaggcgcgg                                                                                                             <  2:90971‑M1/150‑36 (MQ=255)
                                           ttCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcggg                                                                                                            <  2:2021610‑M1/150‑37 (MQ=255)
                                           ttCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcggg                                                                                                            >  1:1481330‑M2/1‑145 (MQ=255)
                                              ccAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGTAAAGGTGACGTAGCCCGTCGCAGGAgcgggttt                                                                                                         >  1:2305502‑M2/1‑142 (MQ=255)
                                                       aaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgggtttacccgcgaa                                                                                                <  1:1503717‑M2/150‑18 (MQ=255)
                                                         cGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTACACTTTTCACCTTCTCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAaga                                                                                              <  2:566010/150‑1 (MQ=255)
                                                          gCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgc                                                                                                               <  2:1135549‑M1/132‑34 (MQ=255)
                                                          gCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgc                                                                                                               >  1:1135549‑M1/1‑99 (MQ=255)
                                                          gCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAagag                                                                                             <  2:2431477/150‑1 (MQ=255)
                                                          gCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAagag                                                                                             >  2:741943/1‑150 (MQ=255)
                                                            ttGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCaaa                                                                                                                                       >  1:1094481/1‑106 (MQ=255)
                                                            ttGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCaaa                                                                                                                                       <  2:1094481/106‑1 (MQ=255)
                                                                  tCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGccc                                                                                     >  2:2152926/1‑150 (MQ=255)
                                                                   cGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCa                                                                                    <  1:741943/150‑1 (MQ=255)
                                                                   cGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCa                                                                                    <  1:99277/150‑1 (MQ=255)
                                                                     agcTTGGTGATGAAGGAAGTGATCCAGCATCGGCTGCAAAAACCCGGCAGC‑‑TTTGGTCTTGTGCCGGCCCCTTCACAGTTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAgg                                                                                  <  1:1741210/148‑1 (MQ=19)
                                                                     aCCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAgg                                                                                  <  2:247803/150‑1 (MQ=255)
                                                                     aCCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAgg                                                                                  <  1:2201834/150‑1 (MQ=255)
                                                                     aCCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAgg                                                                                  <  1:319632/150‑1 (MQ=255)
                                                                         tGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGNAAc                                                                              <  2:1503717/150‑1 (MQ=255)
                                                                              atgaAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTc                                                                         <  2:778378/150‑1 (MQ=255)
                                                                               tgaAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCa                                                                        <  2:581972/150‑1 (MQ=255)
                                                                                 aaggaagTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATc                                                                      <  1:1567533/150‑1 (MQ=255)
                                                                                     aagTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGcc                                                                  <  1:2152926/150‑1 (MQ=255)
                                                                                           tCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGaaa                                                            >  1:2256820/1‑150 (MQ=255)
                                                                                           tCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGACTAGAGGCCGGCCCAGGCAACCTTCCATCTGCCCTCaaa                                                            >  1:1840378/1‑150 (MQ=39)
                                                                                           tCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCATGGAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGaaa                                                            >  1:1593064/1‑150 (MQ=255)
                                                                                             cAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAAAGGCCGGCCTAGGCAACCTTTCATCTGCCCTGAAAAc                                                          >  2:1198807/1‑150 (MQ=255)
                                                                                              aGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGGTTACCCGCGAAGAGGCCGGCCCAGCCAACCTTTCATCTGCCCTGAAAAcg                                                         >  2:741763/1‑150 (MQ=255)
                                                                                               ggATCGGGTGC‑‑TCCTGCGGTGGCTGTTTCTTCTGGTGCCGACCCAGTCGCGGGTAAATCTGCGCCTGCAAAGGGGACGTAGCTCGTCGCAGGAGCGTGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAAcgc                                                        <  1:1998470/150‑1 (MQ=255)
                                                                                                 aTCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAggccggcc                                                                                      <  2:128299/118‑1 (MQ=255)
                                                                                                 aTCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAggccggcc                                                                                      >  1:128299/1‑118 (MQ=255)
                                                                                                     gTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGt                                                                                                                                    <  1:473779/68‑1 (MQ=255)
                                                                                                     gTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGt                                                                                                                                    >  2:473779/1‑68 (MQ=255)
                                                                                                     gTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGcacac                                                  >  2:1469269/1‑149 (MQ=255)
                                                                                                      ttGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATc                                                 >  1:1800148/1‑150 (MQ=255)
                                                                                                         c‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCa                                              <  2:1800148/150‑1 (MQ=255)
                                                                                                           ccccccATGGAGCTCTTTCCTCATGTGACGGACCCTTAGCGGGTCAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGTCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCa                                              <  2:1593064/148‑1 (MQ=255)
                                                                                                            tCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcag                                                                                                                   >  1:2379902‑M1/1‑51 (MQ=255)
                                                                                                            tCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcag                                                                                                                   <  2:2379902‑M1/80‑30 (MQ=255)
                                                                                                              cccATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                             >  1:686137‑M1/1‑49 (MQ=255)
                                                                                                              cccATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                             <  2:686137‑M1/84‑36 (MQ=255)
                                                                                                              cccATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGt                                                                                                                                    >  1:1080821/1‑61 (MQ=255)
                                                                                                              cccATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGt                                                                                                                                    <  2:1080821/61‑1 (MQ=255)
                                                                                                              cccATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCa                                           >  1:1325412/1‑150 (MQ=255)
                                                                                                                                  tgtgCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCg                       <  2:2256820/150‑1 (MQ=255)
                                                                                                                                     gCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGa                                                              <  2:2328576/108‑1 (MQ=255)
                                                                                                                                     gCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGa                                                              >  1:2328576/1‑108 (MQ=255)
                                                                                                                                     gCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGtt                    <  2:1827425/150‑1 (MQ=255)
                                                                                                                                      ccGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAAGTGACGTAGCCCGGCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCCACCTTTCCTCTGCCCTGAAAACGCGCCCCCCGCACCAACCCCCCCCACCCTGGTCCCGTTg                   >  2:2419645/1‑150 (MQ=25)
                                                                                                                                       cGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCaa                                                                               <  2:629092/89‑1 (MQ=255)
                                                                                                                                       cGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCaa                                                                               >  1:629092/1‑89 (MQ=255)
                                                                                                                                           cccTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCAcc              >  1:1340813/1‑150 (MQ=255)
                                                                                                                                           ccattCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGCGCCACCCAGGCCACCTTTCCTCTGCCCTGAAAACGCGAATAGCGCACCAACCCAGCCCACCGTGATCCGGTTGCcccc              >  1:1983153/4‑150 (MQ=255)
                                                                                                                                              ttCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtt           <  1:1503717‑M1/150‑134 (MQ=255)
                                                                                                                                                   gggTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCcgccgc      >  1:1253896/1‑150 (MQ=255)
                                                                                                                                                    ggTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCgccgcc     <  1:419410/150‑1 (MQ=255)
                                                                                                                                                       aaaCCCGCtcctgcgacgggctacgtcacctttacaggcgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttcacattgg  >  1:2305502‑M1/1‑8 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGC‑‑TCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCT  >  NC_002947/4149929‑4150227

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