Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,741,967 T→C 5.9% intergenic (‑35/+378) PP_3305 ← / ← PP_3306 TerC family membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,741,9670TC 5.9% 204.1 / 2.1 68intergenic (‑35/+378)PP_3305/PP_3306TerC family membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/27);  new base C (2/2);  total (39/29)
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.67e-01

GCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAAACTAACTTCGCA  >  NC_002947/3741821‑3742113
                                                                                                                                                  |                                                                                                                                                   
gCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATg                                                                                                                                                    <  5:261357/148‑1 (MQ=255)
      cTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGa                                                                                                                                             >  7:79196/1‑149 (MQ=255)
      cTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGa                                                                                                                                             >  5:57131/1‑149 (MQ=255)
         aGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGa                                                                                                                                                   <  3:358222/140‑1 (MQ=255)
         aGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGa                                                                                                                                                   >  4:358222/1‑140 (MQ=255)
          gTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATg                                                                                                                                         >  6:257692/1‑149 (MQ=255)
           tCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGCCATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATgg                                                                                                                                        >  6:81594/1‑149 (MQ=255)
            cAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGc                                                                                                                                       >  2:112525/1‑149 (MQ=255)
                   cAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATaa                                                                                                                               >  4:195481/1‑149 (MQ=255)
                   cAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATaa                                                                                                                               >  4:201910/1‑149 (MQ=255)
                    aaCAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGACtatac                                                                                                                                             <  5:68827/135‑5 (MQ=255)
                    aaCAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGACtatac                                                                                                                                             >  6:68827/1‑131 (MQ=255)
                       aGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGt                                                                                                                           <  5:257692/149‑1 (MQ=255)
                             aGCAGAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCagag                                                                                                                      <  5:69661/148‑1 (MQ=255)
                              gCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGACAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAggt                                                                                                                    >  4:22545/1‑148 (MQ=255)
                               cacaACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAgggg                                                                                                                   >  6:293789/1‑149 (MQ=255)
                               cacaACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAgggg                                                                                                                   >  3:85239/1‑149 (MQ=255)
                                  aacaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGa                                                                                                                                             >  3:195164/1‑121 (MQ=255)
                                  aacaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGa                                                                                                                                             <  4:195164/121‑1 (MQ=255)
                                  aacaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCg                                                                                                                                     >  5:185901/1‑128 (MQ=255)
                                  aacaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCg                                                                                                                                     <  6:185901/128‑1 (MQ=255)
                                       cGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGt                                                                                                           >  2:398986/1‑149 (MQ=255)
                                       cGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGt                                                                                                           >  4:250023/1‑149 (MQ=255)
                                         aaaaaGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGACtataca                                                                                                                                            >  8:259262/1‑110 (MQ=39)
                                         aaaaaGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGACtataca                                                                                                                                            <  7:259262/115‑6 (MQ=39)
                                                  agccaaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAg                                                                                                <  1:398986/149‑1 (MQ=255)
                                                    ccaaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGaa                                                                                              <  3:201910/149‑1 (MQ=255)
                                                    ccaaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGaa                                                                                              <  3:195481/149‑1 (MQ=255)
                                                      aaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAAtt                                                                                            >  4:347104/1‑149 (MQ=255)
                                                      aaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCACTCCTGTATAGTCGTAGTTGAACAACTATGACTCCTACATGGCGGTATAACAGTCAGAGGGGCCCGGCGTCCCGAATAAAGAAtt                                                                                            >  6:317293/1‑149 (MQ=255)
                                                         gccaGAGGTCAGTAACAATAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTAt                                                                                         <  4:275375/149‑1 (MQ=255)
                                                                gTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAgcgc                                                                                  <  3:347104/149‑1 (MQ=255)
                                                                gTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAgcgc                                                                                  >  2:194446/1‑149 (MQ=255)
                                                                      aCCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTa                                                                             <  5:317293/148‑1 (MQ=255)
                                                                          aGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTg                                                                        <  1:194446/149‑1 (MQ=255)
                                                                          aGAACGTCTTGCTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTTTTTAGTCGTAGTTGAACAACTTTGACTCCGACTTGGCGGTTTAACCGTCAGAGGGCCCCGGCGACACGAATAAAGAATTTTTGTGGCGGTAGTAGGCt                                                                         >  6:319885/1‑148 (MQ=255)
                                                                           gAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTg                                                                        >  7:338559/1‑148 (MQ=255)
                                                                                tCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAgcgc                                                                                  >  4:282233/1‑133 (MQ=255)
                                                                                tCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAgcgc                                                                                  <  3:282233/133‑1 (MQ=255)
                                                                                tCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTa                                                                  <  5:293789/149‑1 (MQ=255)
                                                                                  ttGGTAGAGATACGTGCTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAg                                                                 >  7:89080/1‑148 (MQ=255)
                                                                                    ggTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGaa                                                               >  7:317204/1‑148 (MQ=255)
                                                                                     gTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGaaaa                                                             >  1:146458/1‑149 (MQ=255)
                                                                                            tACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTAc                                                      <  8:338559/149‑1 (MQ=255)
                                                                                                tACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGTGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAAt                                                  <  5:81594/149‑1 (MQ=255)
                                                                                                 aCTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCAGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACaaa                                                   >  1:324892/1‑148 (MQ=255)
                                                                                                  cTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATc                                                 <  6:57131/148‑1 (MQ=255)
                                                                                                   tAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTAt                                                                                         >  8:236854/1‑107 (MQ=255)
                                                                                                   tAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTa                                                       <  8:111408/141‑1 (MQ=255)
                                                                                                   tAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTa                                                       >  7:111408/1‑141 (MQ=255)
                                                                                                   tAAGGGCAGACATGAGCCCCTCCTGTATAGACGTAGATGAACAACTATGACTCCGACATCGCGGTATAACCGTCAGACGGGCCCGCCGTCACGAATAAAGAATATAACGAGCGCGAGAAGGCTGAAAGAACAAAACATTTACAAATCt                                                >  5:197555/1‑148 (MQ=255)
                                                                                                        gCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCg                                                                                       >  4:166170/1‑104 (MQ=255)
                                                                                                         cataCATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCAt                                         <  4:85239/146‑1 (MQ=255)
                                                                                                         cAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCATAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCAt                                         <  7:314757/149‑1 (MQ=255)
                                                                                                              aTGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACAttt                                                        <  4:61391/129‑1 (MQ=255)
                                                                                                              aTGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACAttt                                                        >  3:61391/1‑129 (MQ=255)
                                                                                                                 aGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCATGGACAAgat                                 >  1:30731/1‑147 (MQ=255)
                                                                                                                   ccccTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCATGGACAAgattt                               <  2:30731/149‑5 (MQ=255)
                                                                                                                     cctcctGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTaaaa                             >  7:110220/1‑149 (MQ=255)
                                                                                                                      ctcctGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAAc                            <  3:250023/149‑1 (MQ=255)
                                                                                                                              taGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTg                       >  1:187413/1‑146 (MQ=255)
                                                                                                                              taGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCa                     <  2:187413/148‑1 (MQ=255)
                                                                                                                                 tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTg                  >  3:209332/1‑148 (MQ=255)
                                                                                                                                 tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGAAGTg                  >  1:306132/1‑148 (MQ=255)
                                                                                                                                      gTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGaaaaa             <  3:9014/148‑1 (MQ=255)
                                                                                                                                       ttGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAaact           <  8:89080/149‑1 (MQ=255)
                                                                                                                                         gAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAaactaa         >  1:104986/1‑149 (MQ=255)
                                                                                                                                          aacaacTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAaactaac        >  1:236868/1‑149 (MQ=255)
                                                                                                                                                  tGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAAACTAACTTCGca  <  2:104986/147‑1 (MQ=255)
                                                                                                                                                  tGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAAACTAACTTCGca  <  8:110220/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAAACTAACTTCGCA  >  NC_002947/3741821‑3742113

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: