Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,039,245 A→G 6.3% intergenic (‑265/+29) PP_5290 ← / ← pyrE transposase/orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,039,2450AG 6.3% 185.7 / 2.0 63intergenic (‑265/+29)PP_5290/pyrEtransposase/orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base G (2/2);  total (34/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 = 9.48e-01

GGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGATAC  >  NC_002947/6039101‑6039392
                                                                                                                                                |                                                                                                                                                    
ggTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTc                                                                                                                                                  >  6:194599/1‑149 (MQ=255)
ggTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTc                                                                                                                                                  >  7:65018/1‑149 (MQ=255)
  tGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGc                                                                                                                                                <  6:144171/149‑1 (MQ=255)
        ctcCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCtttt                                                                                                                                           <  7:106703/148‑1 (MQ=255)
           cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCTTCTTCAGAGGGTCGCCttttttt                                                                                                                                        <  8:254815/148‑3 (MQ=255)
           cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCAt                                                                                                                                       >  3:359806/1‑149 (MQ=255)
                   tGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGc                                                                                                                               >  5:139693/1‑149 (MQ=255)
                        tGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGa                                                                                                                          >  2:63505/1‑149 (MQ=255)
                             gAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGa                                                                                                                     <  3:46682/149‑1 (MQ=255)
                               aCAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGc                                                                                                                               <  6:315904/137‑1 (MQ=255)
                               aCAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGc                                                                                                                               >  5:315904/1‑137 (MQ=255)
                                cAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATcc                                                                                                                  <  8:65018/149‑1 (MQ=255)
                                 aGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATccc                                                                                                                 >  7:302585/1‑149 (MQ=255)
                                 aGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATccc                                                                                                                 >  6:362237/1‑149 (MQ=255)
                                   tCAGAATCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCATACCATACAAGCGCGTTTACCCGCGAAAGGGACAGATAAGTCGGCACAAAAAAACGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCATCGGACAGGATCAGATCCCg                                                                                                                >  6:233463/1‑148 (MQ=25)
                                        aGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTgc                                                                                                          <  5:86035/149‑1 (MQ=255)
                                                 tGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTgg                                                                                                                                 >  3:63861/1‑117 (MQ=255)
                                                 tGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTgg                                                                                                                                 <  4:63861/117‑1 (MQ=255)
                                                 tGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTa                                                                                                  <  5:362237/148‑1 (MQ=255)
                                                  ggCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCTTCTTCAGAGGGTCGCCttt                                                                                                                                            >  3:319513/1‑105 (MQ=21)
                                                  ggCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCTTCTTCAGAGGGTCGCCttt                                                                                                                                            <  4:319513/105‑1 (MQ=21)
                                                  ggCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATccc                                                                                                                 >  5:66191/1‑132 (MQ=255)
                                                  ggCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAgg                                                                                                >  6:44118/1‑149 (MQ=255)
                                                  ggCCCTCCTCCTACCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATccc                                                                                                                 <  6:66191/132‑1 (MQ=255)
                                                    cccTCCTCCTCCCTTCAAGTCCTACCATACAAGGGGTTTTACCCGCGAATGGGCCATATCAGGCGGCACAAAAAATGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGAt                                                                                                                         <  2:179948/122‑1 (MQ=37)
                                                    cccTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGAt                                                                                                                         >  1:179948/1‑122 (MQ=255)
                                                         ctcctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGa                                                                                          <  5:44118/148‑1 (MQ=255)
                                                                       tCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCt                                                                           <  5:194599/149‑1 (MQ=255)
                                                                         cTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGATAGCCGGCAGATGCTGt                                                                         <  2:217668/149‑1 (MQ=255)
                                                                         cTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGATAGCCGGCAGATGCTGt                                                                         >  1:217668/1‑149 (MQ=255)
                                                                          tACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACa                                                                                        >  7:194212/1‑133 (MQ=255)
                                                                          tACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACa                                                                                        <  8:194212/133‑1 (MQ=255)
                                                                           aCCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGttt                                                                       >  8:302558/1‑149 (MQ=255)
                                                                            ccATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCTGGCGGCACAAAAATAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGTAGATGCTGttt                                                                       <  5:233463/148‑1 (MQ=255)
                                                                             ccctttcgCGGGTAAACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGttt                                                                       >  7:120089/8‑147 (MQ=255)
                                                                             cATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCa                                                                     <  4:359806/149‑1 (MQ=255)
                                                                               tACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGc                                                                   >  8:345837/1‑149 (MQ=255)
                                                                                  aaGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCAGTTCTGCGCCCGGTAGGCCTCGACAGCCGGCAGCTGCTGTTTCAGCtgt                                                                >  8:278380/1‑149 (MQ=255)
                                                                                    gCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgtg                                                               >  8:221098/1‑148 (MQ=255)
                                                                                       ggTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATc                                                           <  8:302585/149‑1 (MQ=255)
                                                                                       ggTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATc                                                           <  7:302558/149‑1 (MQ=255)
                                                                                        gTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGa                                                                               <  8:107943/128‑1 (MQ=255)
                                                                                        gTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGa                                                                               >  7:107943/1‑128 (MQ=255)
                                                                                        gTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGc                                                                   <  4:201576/140‑1 (MQ=255)
                                                                                        gTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGc                                                                   >  3:201576/1‑140 (MQ=255)
                                                                                              ccGCGAAAGGGCCAGATCAGGCGTCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGc                                                     <  8:297930/148‑1 (MQ=255)
                                                                                              ccGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGc                                                     >  7:297930/1‑148 (MQ=255)
                                                                                                gcgAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAg                                                  >  5:157372/1‑149 (MQ=255)
                                                                                                gcgAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATACCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGGTTCAGATGTGGATCGTCAGCAAg                                                  >  3:278665/1‑149 (MQ=255)
                                                                                                    aaGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAAc                                              <  7:345837/149‑1 (MQ=255)
                                                                                                       ggCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGc                                                                                  >  6:143506/1‑110 (MQ=255)
                                                                                                       ggCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGc                                                                                  <  5:143506/110‑1 (MQ=255)
                                                                                                       ggCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTcc                                           <  7:221098/149‑1 (MQ=255)
                                                                                                        gCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTc                                            >  5:295095/1‑147 (MQ=255)
                                                                                                                 ggcggcACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGGCCGGTACGGCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTggg                                >  1:116727/1‑149 (MQ=255)
                                                                                                                  gcggcACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCa                                          <  6:359764/139‑1 (MQ=255)
                                                                                                                  gcggcACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCa                                          >  5:359764/1‑139 (MQ=255)
                                                                                                                       acaAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTc                                            >  1:174526/1‑132 (MQ=255)
                                                                                                                         aaaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGGCAGGAACTCCAGCCCCTGGGTCAGCGa                         >  8:285810/1‑148 (MQ=255)
                                                                                                                                ggCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGGTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGc                 >  7:132260/1‑149 (MQ=255)
                                                                                                                                 gCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCt                <  6:157372/149‑1 (MQ=255)
                                                                                                                                   gACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGa              <  8:120089/149‑1 (MQ=255)
                                                                                                                                        tctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGaccacc         >  6:165665/1‑149 (MQ=255)
                                                                                                                                           gaagaaGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCgg       >  2:296324/1‑148 (MQ=255)
                                                                                                                                            aagaagGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGa     >  5:349519/1‑149 (MQ=255)
                                                                                                                                            aagaagGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGa     >  3:164737/1‑149 (MQ=255)
                                                                                                                                                agGTCGCCTTTTTCATTCCTGGGCCAGGATCACATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCAACCTGGGTCAGCGAAACAATGCTGACCACCGGGATAc  <  1:63505/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGATAC  >  NC_002947/6039101‑6039392

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: