Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,737,602 G→C 10.5% intergenic (‑157/+45) 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,6020GC10.5% 177.7 / 13.1 67intergenic (‑157/+45)hutU/PP_5034urocanate hydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/37);  new base C (3/4);  total (26/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

CCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACC  >  NC_002947/5737455‑5737741
                                                                                                                                                    |                                                                                                                                              
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)
                                                                                                                                                    |                                                                                                                                              
CCTCTGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTC‑‑GATCTCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACC  >  NC_002947/5737455‑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: