Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,240,000 T→C 7.3% intergenic (+51/+335) dadX → / ← cvrA alanine racemase, catabolic, PLP‑binding/putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,240,0000TC 7.3% 162.8 / 3.1 55intergenic (+51/+335)dadX/cvrAalanine racemase, catabolic, PLP‑binding/putative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (23/28);  new base C (2/2);  total (25/30)
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.02e-01

GGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTACATAAAACGCTTACGCTTCATTGTT  >  NC_000913/1239853‑1240144
                                                                                                                                                   |                                                                                                                                                
ggAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAatcatc                                                                                                                                                 <  8:42632/149‑1 (MQ=255)
          aTTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTg                                                                                                                                       <  1:155444/149‑1 (MQ=255)
                   gTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAAc                                                                                                                              >  7:82259/1‑149 (MQ=255)
                    tCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACt                                                                                                                             <  7:130228/149‑1 (MQ=255)
                    tCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACt                                                                                                                             <  2:261035/149‑1 (MQ=255)
                      gccgccgcTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTg                                                                                                                            >  4:95363/1‑148 (MQ=255)
                      gccgccgcTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTg                                                                                                                            >  6:393923/1‑148 (MQ=255)
                       ccgccgcTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTg                                                                                                                            >  8:87211/1‑147 (MQ=255)
                        cgccgcTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAg                                                                                                                         <  5:393923/149‑1 (MQ=255)
                         gccgcTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAgt                                                                                                                        <  2:239407/149‑1 (MQ=255)
                            gcTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAgttgt                                                                                                                     <  7:459683/149‑1 (MQ=255)
                             cTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAgttgtt                                                                                                                    <  7:87211/149‑1 (MQ=255)
                                ccGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTaa                                                                                                                  <  4:303024/148‑1 (MQ=255)
                                   gAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAatcat                                                                                                                                                  >  2:344030/1‑113 (MQ=255)
                                   gAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAatcat                                                                                                                                                  <  1:344030/113‑1 (MQ=255)
                                   gAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGGTAAtgt                                                                                                               >  3:339329/1‑148 (MQ=255)
                                        gTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACaga                                                                                                         <  3:95363/149‑1 (MQ=255)
                                            gCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAgttgt                                                                                                                     >  7:368552/1‑133 (MQ=255)
                                            gCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAgttgt                                                                                                                     <  8:368552/133‑1 (MQ=255)
                                              tATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATcatccatc                                                                                                                                             >  3:195443/1‑107 (MQ=255)
                                              tATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATcatccatc                                                                                                                                             <  4:195443/107‑1 (MQ=255)
                                              tATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCAtt                                                                                                   >  1:303326/1‑149 (MQ=255)
                                                  aGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGcc                                                                                                >  5:220823/1‑148 (MQ=255)
                                                   gTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCa                                                                                              <  2:230903/149‑1 (MQ=255)
                                                        tgtgCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATa                                                                                         <  6:99359/149‑1 (MQ=255)
                                                        tgtgCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATa                                                                                         <  5:303579/149‑1 (MQ=255)
                                                             gcgcTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTc                                                                                    <  2:431780/149‑1 (MQ=255)
                                                              cgcTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTTTTAATATGACAGAGCCATTGCCCATGATAGTGTcc                                                                                   <  8:82259/149‑1 (MQ=255)
                                                               gcTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCa                                                                                  <  3:296414/149‑1 (MQ=255)
                                                                      cTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAg                                                                           >  6:68513/1‑149 (MQ=255)
                                                                         cgcgTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGAt                                                                        <  8:398553/149‑1 (MQ=255)
                                                                          gcgTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAgg                                                                          <  5:150736/146‑1 (MQ=255)
                                                                          gcgTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAgg                                                                          >  6:150736/1‑146 (MQ=255)
                                                                                      tgACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCa                                                                                              <  5:207397/114‑1 (MQ=255)
                                                                                      tgACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCa                                                                                              >  6:207397/1‑114 (MQ=255)
                                                                                         cGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCCTTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTcccc                                                         <  4:339329/148‑1 (MQ=255)
                                                                                                 cTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTg                                                 >  6:313086/1‑148 (MQ=255)
                                                                                                  ttgttgTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATgg                                                                      >  8:30572/1‑126 (MQ=255)
                                                                                                  ttgttgTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATgg                                                                      <  7:30572/126‑1 (MQ=255)
                                                                                                     ttgtAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACt                                             >  4:18486/1‑148 (MQ=255)
                                                                                                      tgtAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGAc                                                                                                            <  8:317723/84‑1 (MQ=255)
                                                                                                      tgtAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGAc                                                                                                            >  7:317723/1‑84 (MQ=255)
                                                                                                           gCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGca                                      >  6:159652/1‑149 (MQ=255)
                                                                                                           gCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGca                                      >  5:61861/1‑149 (MQ=255)
                                                                                                            ccGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCAtt                                                                                >  2:85967/1‑106 (MQ=255)
                                                                                                            ccGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCAtt                                                                                <  1:85967/106‑1 (MQ=255)
                                                                                                             cGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCACCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGa                                                                     >  6:133496/1‑116 (MQ=255)
                                                                                                             cGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCACCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGa                                                                     >  8:428538/1‑116 (MQ=255)
                                                                                                             cGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCACCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGa                                                                     <  7:428538/116‑1 (MQ=255)
                                                                                                             cGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCACCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGa                                                                     <  5:133496/116‑1 (MQ=255)
                                                                                                                      ggCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTa                           <  6:220823/149‑1 (MQ=255)
                                                                                                                      ggCAACGTCTTCTGGGTGCAAAAAAATCATCCATCAGGCTGGTCAGAAAATGTAGTTGTTAATGTGACAGAGCCATTGACCATGATAGTGTCCATTAAAAGGATGGACACTATTTACCCGGAACATGAACTCAACGCACAGGCGttat                            >  6:178216/1‑146 (MQ=255)
                                                                                                                       gCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTa                           >  3:169130/1‑148 (MQ=255)
                                                                                                                               ttctGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTACATAAAacg                  >  8:62483/1‑149 (MQ=255)
                                                                                                                                 ctGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTACATAAAacgc                 <  5:159652/148‑1 (MQ=255)
                                                                                                                                 ctGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTACATAAAacg                  <  3:18486/147‑1 (MQ=255)
                                                                                                                                  tGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTACATAAAacgct                >  2:450601/1‑148 (MQ=255)
                                                                                                                                    ggTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTACATAAAacgcttac             >  1:157974/1‑149 (MQ=255)
                                                                                                                                    ggTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATCTGACAGAGCCATTGCCCATGATAGTGTCCATTTCAAGGATGGACACTATTTCCCCGGAACCTGACCTCACCGCACAGGCGTGCTACATAAAacgctta              <  7:118590/148‑1 (MQ=255)
                                                                                                                                              aaTCATCCATCCTGCTGGTCAGCAACTGTAGTTGTTAATTTGACAGAGCCATTGCCAATGATAGTGTACATTAAAATGATGTACAATATTTCCCCTGAACCTGAACTCACCGCACATGCGTTCTACATAAAACGCTTACGCTTCAttgt   <  5:178216/149‑1 (MQ=255)
                                                                                                                                               atcatcCATCCGGCTGGTCCTCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCAACGCACAGGCGTTCTACATAAAACGCTTACGCTTCAttgtt  <  8:366201/149‑1 (MQ=255)
                                                                                                                                               atcatcCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTACATAAAACGCTTACGCTTCAttgtt  <  8:372532/149‑1 (MQ=255)
                                                                                                                                               atcatcCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTACATAAAACGCTTACGCTTCAttgtt  >  1:391108/1‑149 (MQ=255)
                                                                                                                                               atcatcCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTACATAAAACGCTTACGCTTCAttgtt  <  4:169130/149‑1 (MQ=255)
                                                                                                                                                   tccatccGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTaaa                                                                             >  1:78713/1‑70 (MQ=255)
                                                                                                                                                   tccatccGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTaaa                                                                             <  2:78713/70‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGGAACCTGAACTCACCGCACAGGCGTTCTACATAAAACGCTTACGCTTCATTGTT  >  NC_000913/1239853‑1240144

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: