Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,632,498 G→T 6.8% intergenic (+103/+56) PP_4096 → / ← PP_t59 hypothetical protein/tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,632,4980GT 6.8% 222.2 / 3.2 71intergenic (+103/+56)PP_4096/PP_t59hypothetical protein/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base G (32/34);  new base T (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGC  >  NC_002947/4632354‑4632643
                                                                                                                                                |                                                                                                                                                 
tttGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAgg                                                                                                                                               >  4:76160/1‑149 (MQ=255)
 ttGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAgg                                                                                                                                               <  8:133403/148‑1 (MQ=255)
  tGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAgg                                                                                                                                               <  8:297075/147‑1 (MQ=255)
  tGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGc                                                                                                                                             >  8:221726/1‑149 (MQ=255)
   gAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGc                                                                                                                                             <  5:195384/148‑1 (MQ=255)
    aaaTACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCt                                                                                                                                           >  2:161278/1‑149 (MQ=255)
        aCCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCttttt                                                                                                                                       >  7:264321/1‑149 (MQ=255)
        aCCCGTTGACGTGATCAGCGATATGTCCAGGAAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCttttt                                                                                                                                       >  1:227488/1‑149 (MQ=255)
         cccGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCtttttt                                                                                                                                      >  6:136435/1‑149 (MQ=255)
                 cGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGAc                                                                                                                               >  5:136320/1‑148 (MQ=255)
                 cGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGAc                                                                                                                               <  6:136320/148‑1 (MQ=255)
                  gTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCTTCGGGGCCttt                                                                                                                                         <  1:76807/137‑1 (MQ=255)
                  gTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCTTCGGGGCCttt                                                                                                                                         >  2:76807/1‑137 (MQ=255)
                   tGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTa                                                                                                                            <  7:141440/149‑1 (MQ=255)
                        aGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGa                                                                                                                        <  3:76160/148‑1 (MQ=255)
                        aGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGAc                                                                                                                       >  4:85072/1‑149 (MQ=255)
                          cGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACtt                                                                                                                     <  2:298735/149‑1 (MQ=255)
                           gATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTc                                                                                                                    >  3:168694/1‑149 (MQ=255)
                            atatGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCg                                                                                                                   <  2:227488/149‑1 (MQ=255)
                             tatGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGa                                                                                                                  <  3:85072/149‑1 (MQ=255)
                               tGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGAc                                                                                                                 <  7:95509/148‑1 (MQ=255)
                               tGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACt                                                                                                                <  1:161278/149‑1 (MQ=255)
                                gTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACt                                                                                                                >  8:301869/1‑148 (MQ=255)
                                 tCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGAc                                                                                                                       <  7:28404/140‑1 (MQ=255)
                                 tCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGGTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGGATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAAGGACTTACAGACTTCGAcgg                                                                                                               >  6:215881/1‑146 (MQ=255)
                                   cAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACa                                                                                                                          <  7:245350/135‑1 (MQ=255)
                                   cAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACa                                                                                                                          >  8:245350/1‑135 (MQ=255)
                                         gCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGt                                                                                                       <  6:258012/148‑1 (MQ=255)
                                           ccACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGa                                                                                                                  <  3:283619/135‑1 (MQ=255)
                                           ccACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGa                                                                                                                  >  4:283619/1‑135 (MQ=255)
                                           ccACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTg                                                                                                        >  2:292702/1‑145 (MQ=255)
                                                  gTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGCCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTAt                                                                                             <  6:345528/149‑1 (MQ=255)
                                                  gTCACGGCCAGGGCGCCGCGGTGGGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTAt                                                                                             <  6:350992/149‑1 (MQ=255)
                                                   tCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTAtt                                                                                            <  5:136435/149‑1 (MQ=255)
                                                   tCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTAtt                                                                                            >  8:84491/1‑149 (MQ=255)
                                                     aCGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGTCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTTGATGTTTATTTg                                                                                          >  5:101813/1‑149 (MQ=255)
                                                             ggCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGa                                                                                  <  7:135700/149‑1 (MQ=255)
                                                               cgccgcGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGaaa                                                                                <  7:301869/149‑1 (MQ=255)
                                                                 ccgcGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAAc                                                                               >  4:36480/1‑148 (MQ=255)
                                                                    ctttgttgTTTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCTAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGTGAAACgaga                                                                           <  4:104834/147‑1 (MQ=255)
                                                                        gtggtgTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTc                                                                        <  7:221726/148‑1 (MQ=255)
                                                                         tggtgTTTGTACTTGAATTTCGTAGAGGGCTGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGa                                                                      <  3:36480/149‑1 (MQ=255)
                                                                           gtgtTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGACCCGAAGGGCCTTTTTTCAACGACTTAAAGACTTCGACTGAAACCTGTAGATATTTATTTGGAGCGGGAAACGAGAATCGaa                                                                     >  3:17148/1‑148 (MQ=255)
                                                                              tttGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCttttt                                                                                                                                       >  8:316658/1‑79 (MQ=255)
                                                                              tttGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCttttt                                                                                                                                       <  7:316658/79‑1 (MQ=255)
                                                                                tGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTcgc                                                                <  7:336700/148‑1 (MQ=255)
                                                                                    cTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGAccc                                                           >  1:356281/1‑149 (MQ=255)
                                                                                    cTTGAATTTCGTAGAGGGAGGGGTATTGAGTCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAAGCTGTAGATCTTTATTTGGAGCGGGAAACGTGACTCGAACTCGCGAccc                                                           >  4:272990/1‑149 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACaaga                                                                           >  5:47022/1‑132 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACaaga                                                                           <  6:47022/132‑1 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGAcccc                                                          >  1:209432/1‑149 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGAcccc                                                          >  7:312259/1‑149 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGAcccc                                                          >  3:278674/1‑149 (MQ=255)
                                                                                               tAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGAATGAAACATGTAGATCTTTATTTGGAGAGGGAAACGAGACTCGCACTCGCGACCCCGACCTTGGc                                                 >  3:276282/1‑148 (MQ=255)
                                                                                                    ggCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTc                                           <  4:278674/149‑1 (MQ=255)
                                                                                                     gCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCg                                          <  2:356281/149‑1 (MQ=255)
                                                                                                       ggggGATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCggg                                        <  5:224998/149‑3 (MQ=255)
                                                                                                        gggTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGc                                       <  1:292702/149‑1 (MQ=255)
                                                                                                        gggTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCACGACCTTGCCAAGGTCgggg                                       <  4:168694/149‑4 (MQ=255)
                                                                                                           tATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctct                                    <  8:264321/149‑1 (MQ=255)
                                                                                                                 ggCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAAc                              <  2:286038/149‑1 (MQ=255)
                                                                                                                          cATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTAtt                     >  3:5016/1‑149 (MQ=255)
                                                                                                                           aTTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTc                    >  3:309262/1‑149 (MQ=255)
                                                                                                                                aacgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCcgcg               >  3:231912/1‑149 (MQ=255)
                                                                                                                                      aaaaGGCCCTTCGGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCggg                                                                                   >  7:24317/1‑75 (MQ=38)
                                                                                                                                      aaaaGGCCCTTCGGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCggg                                                                                   <  8:24317/75‑1 (MQ=38)
                                                                                                                                      aaaaGGCCCTTCGGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGa         >  3:213332/1‑149 (MQ=255)
                                                                                                                                           gCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGt    <  4:5016/149‑1 (MQ=255)
                                                                                                                                           gCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGt    >  3:131654/1‑149 (MQ=255)
                                                                                                                                            ccccGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGATATACCAACTGAGCTATTCCCGCGTCGTGATGGGTg   >  6:299588/1‑149 (MQ=255)
                                                                                                                                             cccGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGc  >  8:123640/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTTGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGC  >  NC_002947/4632354‑4632643

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: