Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,737,599 G→C 10.6% intergenic (‑154/+48) hutU ← / ← PP_5034 urocanate hydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,737,5990GC10.6% 172.4 / 11.3 65intergenic (‑154/+48)hutU/PP_5034urocanate hydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/36);  new base C (3/4);  total (25/40)
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.63e-01

GGTCCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACC  >  NC_002947/5737452‑5737741
                                                                                                                                                    |                                                                                                                                                 
ggTCCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTAAAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCcgcgc                                                                                                                                                  >  5:383833/1‑149 (MQ=255)
 gTCCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCcgcgcg                                                                                                                                                 <  2:342988/149‑1 (MQ=255)
   ccTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcggg                                                                                                                                               >  6:21016/1‑149 (MQ=255)
         cGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACTCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCgcgggcg                                                                                                                                             <  1:282898/146‑1 (MQ=255)
         aggTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTAGCGTTCCTGGCTTGAGATTTATTTTGGGGTATAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                         <  8:18945/148‑1 (MQ=255)
          ggTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                        >  5:132088/1‑149 (MQ=255)
            tcgtcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTAATTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑g                                                                                                                                    <  8:322807/149‑1 (MQ=255)
              gtcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GAt                                                                                                                                 >  3:200954/1‑149 (MQ=255)
              gtcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCCCGCGCGGCGCTC‑‑GAt                                                                                                                                 <  4:95909/149‑1 (MQ=255)
               tcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GAtc                                                                                                                                >  6:144982/1‑149 (MQ=255)
                 cgcgGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GAtct                                                                                                                               >  2:376461/1‑148 (MQ=255)
                 cgcgGGCGGGTATGGCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTCCCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GAtctc                                                                                                                              <  1:243829/149‑1 (MQ=255)
                        gggTATGTCGTGTTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCCCGCGCGGCGCTC‑‑gg                                                                                                                                   <  6:134363/138‑2 (MQ=255)
                         ggTATGTCGTCGTGAATATACTATCACAACTTTACTTGTATGTAAATGGATGGAAAAACAAAATAATGGGTCGTTAGCGTTACTGGCTTGAGATTTATGTTGGTGTATAGATCGTGCGCCGCGCGGGCTGCTCTC‑‑GATCTCACAGgcgc                                                                                                                      <  4:291258/149‑1 (MQ=255)
                          gTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACCAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCATGAGATTTATGTTTGCGTAAAGATAGAGCCCCGCGCGGGCGGCCCTC‑‑GATCTCACAGgcgc                                                                                                                      >  7:274726/1‑148 (MQ=255)
                               tCGTTGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACGTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCaaaa                                                                                                                 <  5:184256/148‑1 (MQ=255)
                               tCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCaaaaa                                                                                                                >  5:29069/1‑149 (MQ=255)
                                       aaTATACAAACACACCTTTACTTGTATTTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGc                                                                                                        <  3:332236/149‑1 (MQ=255)
                                       aaTATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGAACTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGc                                                                                                        <  3:337674/149‑1 (MQ=255)
                                           taCAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTACTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGaa                                                                                                    >  2:108175/1‑149 (MQ=255)
                                            aCAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                       <  8:109429/116‑1 (MQ=255)
                                            aCAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                       >  7:109429/1‑116 (MQ=255)
                                             cAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAAcc                                                                                                  <  7:209379/149‑1 (MQ=255)
                                             cAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAAcc                                                                                                  >  7:87694/1‑149 (MQ=255)
                                              aaaCACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCt                                                                                                 >  7:30867/1‑149 (MQ=255)
                                               aaCACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCtt                                                                                                <  6:383833/149‑1 (MQ=255)
                                                 cacacCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCa                                                                                              <  4:270182/149‑1 (MQ=255)
                                                   cacCTTTTCTTGTATGTTCAAGCATGGTCAACCAAAATAATTGGAACTTTCCGTTCCTGGCTTGAGATTTATTTTTGGGTATAGATCGAGCGGCGCGCGGGCGGCGCTC‑‑GATTCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAg                                                                                             <  1:253275/149‑1 (MQ=255)
                                                   cacCTTTACTTGTATGCACAAGCATGGACAACCAAAATAATCGTACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAgc                                                                                            <  6:191613/149‑1 (MQ=255)
                                                    acCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAgcg                                                                                           <  4:200954/149‑1 (MQ=255)
                                                     ccTTTACTTGTTTGTACTAGCATGGACATCCAAATTAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGTGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCTCAGGCGCCAAAAATCTACAGCCGAACCTTCAgcgc                                                                                          <  3:120990/149‑1 (MQ=39)
                                                      cTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGcc                                                                                         >  8:225806/1‑149 (MQ=255)
                                                      cTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTACCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAAACTACAGCCGAACATTCAgcgca                                                                                         >  4:176521/1‑148 (MQ=255)
                                                           cTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATc                                                                                     <  8:148721/148‑1 (MQ=255)
                                                           cTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCa                                                                                    <  8:30867/149‑1 (MQ=255)
                                                            ttGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATc                                                                                     >  3:150660/1‑147 (MQ=255)
                                                              gtatgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAg                                                                                  <  8:87694/148‑1 (MQ=255)
                                                               tatgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGtt                                                                                >  2:61700/1‑149 (MQ=255)
                                                                   taCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGa                                                                             <  8:623/148‑1 (MQ=255)
                                                                   taCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGa                                                                             <  1:61700/148‑1 (MQ=255)
                                                                   taCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGa                                                                             <  5:144982/148‑1 (MQ=255)
                                                                     cAAGCATGGACAACCAAAATAATCTTACCTTACCGTTACTGGCTTGAGATTTATGTAGGGGTATAGATCGATCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCa                                                                          <  8:5067/149‑1 (MQ=39)
                                                                     cAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATc                                                                           <  7:84370/148‑1 (MQ=255)
                                                                      aaGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATc                                                                                     <  1:192909/137‑1 (MQ=255)
                                                                      aaGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATc                                                                                     >  2:192909/1‑137 (MQ=255)
                                                                        tcATGGACAACCAAAATAATCGGAACTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATAGAGCGCCGCGCGGGCGGCGCTCTTGATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAg                                                                         <  1:108175/148‑1 (MQ=255)
                                                                         tatgtcgaaCCAAAATACTCGTACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGgc                                                                       <  8:352470/141‑1 (MQ=255)
                                                                           tGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGgcag                                                                     >  4:131289/1‑148 (MQ=255)
                                                                           tGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTAGAGATTTATGTTGGGGTATAGAACGAGCGCCGCGCGGGCGCCGCTC‑‑GATATCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGAGCCCATAAAGTAGGATCAGgc                                                                       >  7:352470/1‑146 (MQ=255)
                                                                                    aaaaTAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCt                                                           <  3:131289/149‑1 (MQ=255)
                                                                                    aaaaTAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCt                                                           >  7:108928/1‑149 (MQ=255)
                                                                                     aaaTAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTg                                                          >  7:185488/1‑149 (MQ=255)
                                                                                       ataatCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGca                                                        <  7:225806/149‑1 (MQ=255)
                                                                                                ccTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGcc                                               <  1:376461/149‑1 (MQ=255)
                                                                                                     ccGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCATAACCAGCCCCAg                                           <  2:4809/148‑1 (MQ=255)
                                                                                                     ccGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAg                                           <  5:21016/148‑1 (MQ=255)
                                                                                                               gCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATtcatca                                 <  8:108928/149‑1 (MQ=255)
                                                                                                                 ttGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCAtt                               >  7:88345/1‑149 (MQ=255)
                                                                                                                 ttGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATAAtt                               >  7:311658/1‑149 (MQ=255)
                                                                                                                  tGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTg                              >  1:331805/1‑149 (MQ=255)
                                                                                                                   gagaTTTATGTTGCGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCTCCAAAAATCTACAGCCGAACCTTCAGCTCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCAtt                               <  8:274726/147‑1 (MQ=255)
                                                                                                                     gaTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGAGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCt                           <  3:176521/149‑1 (MQ=255)
                                                                                                                          aTGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCa                       <  8:185488/148‑1 (MQ=255)
                                                                                                                          aTGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCa                       <  6:150526/148‑1 (MQ=255)
                                                                                                                                       gATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTc          <  6:29069/148‑1 (MQ=255)
                                                                                                                                        aTCGAGCGCCGCCCGCGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCg         >  4:53532/1‑148 (MQ=255)
                                                                                                                                          cGAGCGCCGCCCGCGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGaa       >  4:313751/1‑148 (MQ=255)
                                                                                                                                           gAGCGCCGCCCGCGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGGGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCTTCATTGCCTTCCAGCCGCAAGCAGTCGaaca     >  5:292817/1‑149 (MQ=255)
                                                                                                                                              cgccgcCCGCGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAAcc  <  3:313751/149‑1 (MQ=255)
                                                                                                                                                  gcgcgGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAg               >  8:303666/1‑132 (MQ=255)
                                                                                                                                                  gcgcgGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAg               <  7:303666/132‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGTCCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACC  >  NC_002947/5737452‑5737741

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