Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,880,869 G→C 6.3% intergenic (+138/‑49) PP_4289 → / → uacT membrane protein/uric acid permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,880,8690GC 6.3% 253.2 / 6.4 79intergenic (+138/‑49)PP_4289/uacTmembrane protein/uric acid permease
Reads supporting (aligned to +/- strand):  ref base G (34/40);  new base C (2/3);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

ACTGACAGCACCTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTGGCCTGCAACACGTTCTGCTGA  >  NC_002947/4880727‑4881017
                                                                                                                                              |                                                                                                                                                    
aCTGACAGCACCTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCt                                                                                                                                                >  2:380876/1‑149 (MQ=255)
 cTGACAGCACCTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTg                                                                                                                                               >  2:252631/1‑149 (MQ=255)
  tGACAGCACCTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGc                                                                                                                                              >  8:192298/1‑149 (MQ=255)
     cAGCACCTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGc                                                                                                                                           <  8:375419/149‑1 (MQ=255)
      aGCACCTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCG‑GCAGGCTTCATGCCTGcgcgct                                                                                                                                         <  4:98185/149‑6 (MQ=255)
      aGCACCTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCG‑GCAGGCTTCATGCCTGcgcgct                                                                                                                                         <  3:171136/149‑6 (MQ=255)
       gCACCTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCtt                                                                                                                                         <  4:228401/149‑1 (MQ=255)
       gCACCTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCtt                                                                                                                                         <  1:286300/149‑1 (MQ=255)
       gCACCTCGCCAGGCAACCCGCCACCCAAGGTAATCACCTTCCCCGGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCtt                                                                                                                                         <  2:341930/149‑1 (MQ=255)
        cACCTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTg                                                                                                                                        >  7:44348/1‑149 (MQ=255)
           cTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGAt                                                                                                                                     >  8:69891/1‑149 (MQ=255)
                 aGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCgag                                                                                                                                >  1:274629/1‑148 (MQ=255)
                   gCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTGCATGTGTTACGGAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGaa                                                                                                                              >  1:294877/1‑148 (MQ=255)
                   gCAACCCGCCCCCCAAGGGAAGCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAAt                                                                                                                             <  6:288529/149‑1 (MQ=255)
                      aCCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATaaa                                                                                                                          <  7:397388/149‑1 (MQ=255)
                      aCCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATaaa                                                                                                                          <  8:345571/149‑1 (MQ=255)
                       cccGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGATGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATaaaa                                                                                                                         <  5:390051/149‑1 (MQ=255)
                       cccGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATaaaa                                                                                                                         <  4:21256/149‑1 (MQ=255)
                       cccGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATaaa                                                                                                                          >  5:83254/1‑148 (MQ=255)
                        ccGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATaaaaa                                                                                                                        >  6:162457/1‑149 (MQ=255)
                              cccAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAaaaacaaaa                                                                                                                   >  2:139459/1‑148 (MQ=255)
                              cccAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCACAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAaaaacaaaac                                                                                                                  >  1:270730/1‑149 (MQ=255)
                                   gggAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAg                                                                                                             <  1:173842/149‑1 (MQ=255)
                                     gAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCGTGGCTCCGAGAATAACAACAAAACTCGAGGt                                                                                                           >  8:161665/1‑149 (MQ=255)
                                      aaTCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGtg                                                                                                          <  4:107962/149‑1 (MQ=255)
                                       aTCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAggcg                                                                                                          <  4:263315/148‑3 (MQ=255)
                                       aTCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGtgc                                                                                                         >  7:357999/1‑149 (MQ=255)
                                          aCCTTCCCCTGCAGGAGGGGGATTGCCCGCGAAGAGTCCTATAGCCCGAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGtgctgc                                                                                                      <  5:296393/149‑1 (MQ=255)
                                          aCCTTCCCCTGCAGGAGCGGGCTTGACCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCCAGGtgctg                                                                                                       >  8:367288/1‑148 (MQ=255)
                                                  cTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGa                                                                                              <  1:139459/149‑1 (MQ=255)
                                                  cTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGa                                                                                              <  7:192298/149‑1 (MQ=255)
                                                     cAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCTCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGt                                                                                            <  1:407425/148‑1 (MQ=255)
                                                     cAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCATGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATCTCCGact                                                                                            >  4:311217/1‑146 (MQ=255)
                                                        gAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACg                                                                                        <  2:227014/149‑1 (MQ=255)
                                                              gCTTGCCCGCGAAGAGGCCTATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGc                                                                                  <  1:252631/149‑1 (MQ=255)
                                                              gCTTGCCCGCGAAGAGGCAGATAGCCCCAGCCTCAAGCTCCCTGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGc                                                                                  <  4:264752/149‑1 (MQ=255)
                                                                   cccGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGtata                                                                             >  7:197236/1‑149 (MQ=255)
                                                                   cccGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCAGCATGTACCAGTCACGCAAGGCGtatc                                                                             >  4:57533/1‑148 (MQ=255)
                                                                    ccGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAg                                                                                    >  8:50266/1‑141 (MQ=255)
                                                                    ccGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAg                                                                                    <  7:50266/141‑1 (MQ=255)
                                                                     cgcgAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATc                                                                           <  6:83254/149‑1 (MQ=255)
                                                                       cgAAGAGGCAGATAGCCCCAGCCTCAATCTCCATGTGTTACGCAACAAGCAATCATCGAGCGCGCAGGCATGAAGACTTCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATccc                                                                         <  3:57533/149‑1 (MQ=255)
                                                                        gAAGAGGCCTATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAAGAATACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATcc                                                                          <  3:311217/147‑1 (MQ=255)
                                                                             ggCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTgc                                                                   <  5:162457/149‑1 (MQ=255)
                                                                              gCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTgcg                                                                  >  6:366776/1‑149 (MQ=255)
                                                                                 gATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGcgccg                                                               <  8:44348/149‑1 (MQ=255)
                                                                                   tAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGcgccgcc                                                             <  5:90181/149‑1 (MQ=255)
                                                                                        ccAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCTGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCcgcgc                                                         >  3:149039/1‑148 (MQ=255)
                                                                                        ccAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCcgcgc                                                         >  3:141004/1‑148 (MQ=255)
                                                                                                aaGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCTGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCtt                                                 <  4:149039/148‑1 (MQ=255)
                                                                                                 aGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCcgcgcgtggtggc                                                 >  3:182220/1‑140 (MQ=255)
                                                                                                  gCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAATCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGcc                                              <  2:270730/149‑1 (MQ=255)
                                                                                                   cTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGccc                                             >  3:206926/1‑149 (MQ=255)
                                                                                                      cATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAg                                                                                    >  2:103465/1‑107 (MQ=255)
                                                                                                      cATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAg                                                                                    <  1:103465/107‑1 (MQ=255)
                                                                                                      cATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTg                                          <  5:366776/149‑1 (MQ=255)
                                                                                                       aTGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGtat                                                                              >  6:287553/1‑112 (MQ=255)
                                                                                                       aTGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGtat                                                                              <  5:287553/112‑1 (MQ=255)
                                                                                                            ttACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCaa                                    <  2:60260/149‑1 (MQ=255)
                                                                                                              aCGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAAc                                   <  4:206926/148‑1 (MQ=255)
                                                                                                                gCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGa                                >  5:287149/1‑149 (MQ=255)
                                                                                                                  aacaaGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCAc                                                                                         <  3:239677/90‑1 (MQ=255)
                                                                                                                  aacaaGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCAc                                                                                         >  4:239677/1‑90 (MQ=255)
                                                                                                                  aacaaGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATc                              >  4:16877/1‑149 (MQ=255)
                                                                                                                    caagcaagCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCt                            >  3:425405/1‑149 (MQ=255)
                                                                                                                      agcaagcaTCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTgg                          <  7:367288/149‑1 (MQ=255)
                                                                                                                         aagcaTCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTg                       <  5:422740/149‑1 (MQ=255)
                                                                                                                         aagcaTCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTg                       >  4:25036/1‑149 (MQ=255)
                                                                                                                         aagcaTCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTg                       >  7:328869/1‑149 (MQ=255)
                                                                                                                         aagcaTCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTg                       >  7:135456/1‑149 (MQ=255)
                                                                                                                         aagcaTCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTg                       >  5:441369/1‑149 (MQ=255)
                                                                                                                           gcaTCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTGGc                     >  4:60278/1‑149 (MQ=255)
                                                                                                                                gAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTGGCCTGc                 <  2:274629/148‑1 (MQ=255)
                                                                                                                                  agcggcAGGCTTCATGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGt                         <  5:2207/134‑1 (MQ=255)
                                                                                                                                  agcggcAGGCTTCATGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGt                         >  6:2207/5‑138 (MQ=255)
                                                                                                                                   gcggcAGGCTTCATGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTGGCCTGCAaca             >  5:264627/4‑149 (MQ=255)
                                                                                                                                     cgcAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTGGCCTGCAACACg           <  2:294877/149‑1 (MQ=255)
                                                                                                                                        aGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTGGCCTGCAACACGTTc        <  3:25036/149‑1 (MQ=255)
                                                                                                                                          gCATGAAGCCTGCCGGTTGGATCCGAGAATAAAAACAAAAGTGGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTGGCCTGCAACACGTTctg      <  1:380876/149‑1 (MQ=255)
                                                                                                                                              gAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTGGCCTGCAACACGTTCTGCTGa  >  5:236447/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
ACTGACAGCACCTCGCCAGGCAACCCGCCCCCCAAGGGAATCACCTTCCCCTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCGATAGCCCCAGCCTCAAGCTCCATGTGTTACGCAACAAGCAAGCATCGAGCGCGCAGGCATGAAGCCTGCCGCTTGGATCCGAGAATAAAAACAAAACTCGAGGTGCTGCATGTCCGAGTCACGCAAGGCGTATATCCCCATTGCGCCGCCGCGCGAACCCTTGCCCCTGTTCCAACTGATCCTGGTTGGCCTGCAACACGTTCTGCTGA  >  NC_002947/4880727‑4881017

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