Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,829,972 T→G 5.6% H192Q (CAT→CAG PP_2482 → molybdenum cofactor biosynthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,829,9720TG 5.6% 257.0 / 4.4 82H192Q (CAT→CAGPP_2482molybdenum cofactor biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base G (2/3);  total (41/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 = 9.95e-01

TTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCGGGGCGT  >  NC_002947/2829824‑2830120
                                                                                                                                                      |                                                                                                                                                    
ttGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCAt                                                                                                                                                      <  6:104789/149‑1 (MQ=255)
ttGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCAt                                                                                                                                                      <  6:355912/149‑1 (MQ=255)
 tgtgCAGGCGCGCCGTGCCGGCATGAGCATCAAGATCTACATGGTGCCGTTGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCAt                                                                                                                                                      <  6:93479/145‑1 (MQ=255)
   aGCACGTGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCAGGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTg                                                                                                                                                   <  5:274735/149‑1 (MQ=255)
        gcgcgcCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGcc                                                                                                                                                >  2:387555/1‑147 (MQ=255)
            gcCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCAt                                                                                                                                          >  4:138974/1‑149 (MQ=255)
            gcCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCAt                                                                                                                                          >  4:7637/1‑149 (MQ=255)
              cGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCAGATGGCCCatgcgca                                                                                                                                        <  3:194252/149‑7 (MQ=255)
                 gCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGcaca                                                                                                                                     <  1:387555/149‑1 (MQ=255)
                  ccGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAg                                                                                                                                    >  1:292698/1‑149 (MQ=255)
                   cGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAgg                                                                                                                                   <  2:334958/149‑1 (MQ=255)
                     gCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCt                                                                                                                                 >  6:122936/1‑149 (MQ=255)
                         gAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGc                                                                                                                                        >  4:370824/1‑138 (MQ=255)
                         gAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGc                                                                                                                                        <  3:370824/138‑1 (MQ=255)
                          aGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTc                                                                                                                             <  1:298312/148‑1 (MQ=255)
                             aTCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTAGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTa                                                                                                                          >  1:290388/1‑149 (MQ=255)
                             aTCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTAc                                                                                                                         <  5:122936/149‑1 (MQ=255)
                             aTCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTAc                                                                                                                         <  3:7637/149‑1 (MQ=255)
                             aTCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTAc                                                                                                                         <  3:138974/149‑1 (MQ=255)
                              tCAAGATCAACATGGTGCCGATGCGTGGTAGCGAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGTGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTAc                                                                                                                         >  1:46652/1‑149 (MQ=255)
                                 aGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCttt                                                                                                                               <  6:334579/139‑1 (MQ=255)
                                 aGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCttt                                                                                                                               >  5:334579/1‑139 (MQ=255)
                                    tCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAGCTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTtc                                                                                                                  <  4:240603/149‑1 (MQ=255)
                                            gTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGc                                                                                                            >  6:245897/1‑147 (MQ=255)
                                               ccGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAgcg                                                                                                       <  6:307937/149‑1 (MQ=255)
                                                cGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAgcgc                                                                                                      <  2:263978/149‑1 (MQ=255)
                                                cGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAgcgc                                                                                                      <  8:135253/149‑1 (MQ=255)
                                                 gATGCGTGGTAGCAACCTCGCTCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCAACCAGTAGCGTATTATCGAACTTATGAGCATGGGCCATCTGGCCCAGCATGCACAGGATTTTCCTACCCAGTTCGTCGGCAAgcgc                                                                                                      >  7:160283/1‑148 (MQ=255)
                                                 gATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAgcgc                                                                                                      >  8:195327/1‑148 (MQ=255)
                                                         gTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGctgctg                                                                                             <  3:244999/149‑1 (MQ=255)
                                                         gTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGctgctg                                                                                             <  7:276766/149‑1 (MQ=255)
                                                         ctAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGctgctg                                                                                             <  5:379198/148‑1 (MQ=255)
                                                          tAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTgg                                                                                            <  4:125785/149‑1 (MQ=255)
                                                            gCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGCGctgctg                                                                                             >  6:379198/1‑146 (MQ=255)
                                                                     gATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCg                                                                                   <  2:292698/148‑1 (MQ=255)
                                                                      aTCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCACCGCGAGCTGCTGGAGCTTATCGGc                                                                                 <  4:381289/149‑1 (MQ=255)
                                                                            ttttGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGc                                                                                 >  4:119266/1‑143 (MQ=255)
                                                                            ttttGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGc                                                                                 <  3:119266/143‑1 (MQ=255)
                                                                            ttttGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCAAGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGt                                                                           >  1:314511/1‑149 (MQ=255)
                                                                                  cccATTTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGc                                                                      <  5:44472/149‑1 (MQ=255)
                                                                                  cccATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTAcc                                                                                                                        >  8:314046/1‑99 (MQ=255)
                                                                                  cccATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTAcc                                                                                                                        <  7:314046/99‑1 (MQ=255)
                                                                                   gcATGTTGGAGTTCTGCCCTACTCCGGGCAACGCTCTGCGTTTTATCGAACTTATGCGCCTGGGGCAGCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGGGAGCGTTATGcc                                                                     <  8:252907/148‑1 (MQ=255)
                                                                                   ccATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGcc                                                                     >  3:186409/1‑149 (MQ=255)
                                                                                       gTTGGCGTTCTGCCTTCAGTAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAgct                                                                                                 <  6:179261/117‑1 (MQ=255)
                                                                                       gTTGGCGTTCTGCCTTCAGTAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAgct                                                                                                 >  5:179261/1‑117 (MQ=255)
                                                                                             gTTCTGCCTTGAGCAGGGCTTCGAGTTGGGTTTTATGGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTAGCGAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAgg                                                           <  8:160283/149‑1 (MQ=255)
                                                                                             gTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTAt                                                                        >  2:366111/1‑136 (MQ=255)
                                                                                             gTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTAt                                                                        <  1:366111/136‑1 (MQ=255)
                                                                                             gTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAgg                                                           >  2:266600/1‑149 (MQ=255)
                                                                                                 tGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCt                                                         >  4:163245/1‑147 (MQ=255)
                                                                                                  gCCTTCAGCAGGGCTACGAGAAGCGATTTATCGAACTTATGCGCAAGGGCCAACTCGCCCAGCATACACAGGCTGATCGTACCCAGTGCGTAGCCCACCGCCAGCTGCTGGAGCTAATCGGCGAGAGGACAGCCTACACCCAGCCTCAg                                                      >  8:327979/1‑149 (MQ=255)
                                                                                                         gcagGGCTACGAGTTGCGTTTTATCGAACTTCTGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCTGCCAGCGCGAGCTGCTGGAGCTTATCGTCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCgg                                               <  6:263025/149‑1 (MQ=255)
                                                                                                           agGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTa                                             <  2:290388/149‑1 (MQ=255)
                                                                                                             ggCTAGGAGTTGCGTTTTATCGAACTTCTGCGCATGGTCCATCTGGCCGAGCATGCACAGGCTTTTCGTACCCTGGTCGGCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTTTGCGTTCACCCAGGCTGAGGCGCCGGTAGa                                           <  7:327979/149‑1 (MQ=255)
                                                                                                               cTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACt                                         >  7:127303/1‑149 (MQ=255)
                                                                                                                tACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTc                                        >  3:187480/1‑149 (MQ=255)
                                                                                                                  cGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTc                                        <  4:187480/147‑1 (MQ=255)
                                                                                                                   gAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGc                                                                                                                                        >  5:116111/1‑50 (MQ=255)
                                                                                                                   gAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGc                                                                                                                                        <  6:116111/50‑1 (MQ=255)
                                                                                                                           ttttATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGc                             >  3:199685/1‑149 (MQ=255)
                                                                                                                           ttttATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGc                             >  4:338087/1‑149 (MQ=255)
                                                                                                                              tATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGt                           <  5:245897/148‑1 (MQ=255)
                                                                                                                               aTCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTta                         >  5:229162/1‑149 (MQ=255)
                                                                                                                                tCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTta                         >  3:357411/1‑148 (MQ=255)
                                                                                                                                tCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTta                         >  3:106443/1‑148 (MQ=255)
                                                                                                                                tCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTta                         >  3:163733/1‑148 (MQ=255)
                                                                                                                                  gAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAg                      >  5:345280/1‑149 (MQ=255)
                                                                                                                                  gAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAg                      >  7:98017/1‑149 (MQ=255)
                                                                                                                                   aaCTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGc                     >  2:56132/1‑149 (MQ=255)
                                                                                                                                     cTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGt                   <  4:186409/149‑1 (MQ=255)
                                                                                                                                      ttATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTg                              >  7:20885/1‑137 (MQ=255)
                                                                                                                                      ttATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTg                              <  8:20885/137‑1 (MQ=255)
                                                                                                                                      ttATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGtt                  >  2:86703/1‑149 (MQ=255)
                                                                                                                                      ttATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGtt                  >  5:282462/1‑149 (MQ=255)
                                                                                                                                        aTGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGtt                  <  7:285499/147‑1 (MQ=255)
                                                                                                                                           atgctGGGCCAGATGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGa                                                       >  8:51579/5‑107 (MQ=255)
                                                                                                                                           atgctGGGCCAGATGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGa                                                       <  7:51579/103‑1 (MQ=255)
                                                                                                                                            tgctGGGCCAGATGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAg            >  7:31726/4‑149 (MQ=255)
                                                                                                                                              aTGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAgcg          >  8:242328/1‑149 (MQ=255)
                                                                                                                                                gggCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGtt                                                                          <  3:33679/83‑1 (MQ=255)
                                                                                                                                                gggCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGtt                                                                          >  4:33679/1‑83 (MQ=255)
                                                                                                                                                  gCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCggg      <  3:338087/149‑1 (MQ=255)
                                                                                                                                                  gCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCggg      <  6:229162/149‑1 (MQ=255)
                                                                                                                                                  gCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCggg      <  4:232973/149‑1 (MQ=255)
                                                                                                                                                      tCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCGGGGCGt  >  5:317050/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TTGAGCAGGCGCGCCGAGCCGGCATGAGCATCAAGATCAACATGGTGCCGATGCGTGGTAGCAACCTCGATCAAGTTTTGCCCATGTTGGCGTTCTGCCTTCAGCAGGGCTACGAGTTGCGTTTTATCGAACTTATGCGCATGGGCCATCTGGCCCAGCATGCACAGGCTTTTCGTACCCAGTTCGTCGGCCAGCGCGAGCTGCTGGAGCTTATCGGCGAGCGTTATGCCTTCACCCAGGCTGAGGCGCCGGTAGACTCCACCGCGCTGCGTTATAGCGTTGGCGAGCGCGGGGCGT  >  NC_002947/2829824‑2830120

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