Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,971,996 T→A 6.0% L246Q (CTG→CAG)  wecC → UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,971,9960TA 6.0% 436.0 / 14.4 117L246Q (CTG→CAG) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (55/55);  new base A (2/5);  total (57/60)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.10e-01

Only 100 of 118 total aligned reads displayed.

AAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGAT  >  NC_000913/3971849‑3972138
                                                                                                                                                   |                                                                                                                                              
aaGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCt                                                                                                                                                >  1:755148/1‑148 (MQ=255)
    tGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGCTCAGGGGATTAACGTCTGGGAACTGATTCGCCTgg                                                                                                                                              <  2:278565/146‑1 (MQ=255)
    tGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTgg                                                                                                                                              >  1:278565/1‑146 (MQ=255)
        tgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATc                                                                                                                                        <  3:401108/148‑1 (MQ=255)
        tgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATc                                                                                                                                        <  2:1000904/148‑1 (MQ=255)
        tgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATc                                                                                                                                        <  1:33893/148‑1 (MQ=255)
            tcgtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCa                                                                                                                                    <  2:599109/148‑1 (MQ=255)
                              cggcggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttct                                                                                                                  <  2:755148/148‑1 (MQ=255)
                              cggcggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttct                                                                                                                  <  1:285257/148‑1 (MQ=255)
                                gcggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAAt                                                                                                                                         >  4:223638/1‑123 (MQ=255)
                                gcggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAAt                                                                                                                                         <  3:223638/123‑1 (MQ=255)
                                         tgtAAGCTCACCGAAAACGGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTgg                                                                                                                                              <  2:180624/109‑1 (MQ=255)
                                         tgtAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTgg                                                                                                                                              >  1:180624/1‑109 (MQ=255)
                                          gtAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATc                                                                                                                                        >  1:741414/1‑114 (MQ=255)
                                          gtAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATc                                                                                                                                        <  2:741414/114‑1 (MQ=255)
                                               cTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggc                                                                                                        >  2:567425/1‑141 (MQ=255)
                                               cTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggc                                                                                                        <  1:567425/141‑1 (MQ=255)
                                                 cACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGa                                                                                                                                           <  3:180899/104‑1 (MQ=255)
                                                 cACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGa                                                                                                                                           >  4:180899/1‑104 (MQ=255)
                                                   ccGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTgg                                                                                                                                              <  1:393475/99‑1 (MQ=255)
                                                   ccGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTgg                                                                                                                                              >  2:393475/1‑99 (MQ=255)
                                                    cGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTgg                                                                                                                                              <  4:1093756/98‑1 (MQ=255)
                                                    cGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTgg                                                                                                                                              >  3:1093756/1‑98 (MQ=255)
                                                          cAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTgg                                                                                                                                              >  1:704142/1‑92 (MQ=255)
                                                          cAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTgg                                                                                                                                              <  2:704142/92‑1 (MQ=255)
                                                          cAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAatat                                                                                                                     <  4:1022516/117‑1 (MQ=255)
                                                          cAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAatat                                                                                                                     >  3:1022516/1‑117 (MQ=255)
                                                               tCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCAt                                                                                 >  4:1131499/1‑148 (MQ=255)
                                                               tCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGaa                                                                                                                                          <  1:879807/91‑1 (MQ=255)
                                                               tCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGaa                                                                                                                                          >  2:879807/1‑91 (MQ=255)
                                                                    gATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTg                                                                            <  3:1131499/148‑1 (MQ=255)
                                                                         gAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCa                                                                                                                                    <  3:954413/87‑1 (MQ=255)
                                                                          aaTATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccctgg                                                                                                   >  4:945478/1‑119 (MQ=255)
                                                                          aaTATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccctgg                                                                                                   <  3:945478/119‑1 (MQ=255)
                                                                           atatCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCAtt                                                                                >  1:924192/1‑137 (MQ=255)
                                                                           atatCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCAtt                                                                                <  2:924192/137‑1 (MQ=255)
                                                                              tCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                     <  3:293037/81‑1 (MQ=255)
                                                                              tCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                     >  4:293037/1‑81 (MQ=255)
                                                                                   tttGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAtt                                                                                                                    >  4:151250/1‑93 (MQ=255)
                                                                                   tttGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAtt                                                                                                                    <  3:151250/93‑1 (MQ=255)
                                                                                   tttGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCa                                                                                                               >  3:214906/1‑98 (MQ=255)
                                                                                   tttGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCa                                                                                                               <  4:214906/98‑1 (MQ=255)
                                                                                       cTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGt                                                         <  4:718711/148‑1 (MQ=255)
                                                                                         aaTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAatat                                                                                                                     <  2:921450/86‑1 (MQ=255)
                                                                                         aaTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAatat                                                                                                                     >  1:921450/1‑86 (MQ=255)
                                                                                         aaTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCa                                                                                       >  3:904891/1‑116 (MQ=255)
                                                                                         aaTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCa                                                                                       <  4:904891/116‑1 (MQ=255)
                                                                                         aaTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGt                                                                   <  2:927966/136‑1 (MQ=255)
                                                                                         aaTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGt                                                                   >  1:927966/1‑136 (MQ=255)
                                                                                           tGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttcttc                                                                                                                <  3:583587/89‑1 (MQ=255)
                                                                                           tGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttcttc                                                                                                                >  4:583587/1‑89 (MQ=255)
                                                                                             aaTTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGt                                                                   >  1:1069733/1‑132 (MQ=255)
                                                                                             aaTTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGt                                                                   <  2:1069733/132‑1 (MQ=255)
                                                                                              aTTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTaa                                                                                                                        <  3:729985/78‑1 (MQ=255)
                                                                                                tGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCa                                                                                       >  3:901823/1‑109 (MQ=255)
                                                                                                tGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCa                                                                                       <  4:901823/109‑1 (MQ=255)
                                                                                                     cTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTc                                                                                        <  3:754730/103‑1 (MQ=255)
                                                                                                     cTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTc                                                                                        >  4:754730/1‑103 (MQ=255)
                                                                                                      tGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggcc                                                                                                       >  1:1119556/1‑87 (MQ=255)
                                                                                                      tGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggcc                                                                                                       <  2:1119556/87‑1 (MQ=255)
                                                                                                      tGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGt                                                                   >  4:110311/1‑123 (MQ=255)
                                                                                                      tGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGt                                                                   <  3:110311/123‑1 (MQ=255)
                                                                                                         tttGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAtt                                                                                                                    >  1:147211/1‑71 (MQ=255)
                                                                                                                    cAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGcac                                                    >  1:493205/1‑124 (MQ=255)
                                                                                                                    cAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGcac                                                    <  2:493205/124‑1 (MQ=255)
                                                                                                                      ggggATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccc                                                                                                      <  3:700926/72‑1 (MQ=255)
                                                                                                                      ggggATTAACGTCTGGGAACTGATTCGACTGGCGAATCGTCACCCACGCGTTAATATTCTTCAGcctggccc                                                                                                      >  4:700926/1‑72 (MQ=255)
                                                                                                                        ggATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccctgg                                                                                                   <  1:210399/73‑1 (MQ=255)
                                                                                                                        ggATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccctgg                                                                                                   >  2:210399/1‑73 (MQ=255)
                                                                                                                                 tCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttct                                                                                                                  >  2:853596/1‑49 (MQ=255)
                                                                                                                                   tGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGc                                                                                            >  3:435949/1‑69 (MQ=255)
                                                                                                                                   tGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGc                                                                                            <  4:435949/69‑1 (MQ=255)
                                                                                                                                    gggAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGGTGATCCGTGGttt                                                              <  1:846385/98‑1 (MQ=255)
                                                                                                                                      gAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAg          >  1:756617/1‑148 (MQ=255)
                                                                                                                                       aaCTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAgg                                     >  2:217580/1‑120 (MQ=255)
                                                                                                                                       aaCTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAgg                                     <  1:217580/120‑1 (MQ=255)
                                                                                                                                        aCTGATTCGCCTGGCGAATCGTCACCCTCGCGTTa                                                                                                                         <  1:1079571/35‑1 (MQ=255)
                                                                                                                                        aCTGATTCGCCTGGCGAATCGTCACCCTCGCGTTa                                                                                                                         >  2:1079571/1‑35 (MQ=255)
                                                                                                                                        aCTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCgg                                <  1:450443/124‑1 (MQ=255)
                                                                                                                                        aCTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCgg                                >  2:450443/1‑124 (MQ=255)
                                                                                                                                          tGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCg                                                                                           >  4:1002957/1‑63 (MQ=255)
                                                                                                                                          tGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCg                                                                                           <  3:1002957/63‑1 (MQ=255)
                                                                                                                                           gATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggc                                                                                                        >  2:786659/1‑49 (MQ=255)
                                                                                                                                           gATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggc                                                                                                        <  1:786659/49‑1 (MQ=255)
                                                                                                                                           gATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGtt                                                               <  2:1041792/90‑1 (MQ=255)
                                                                                                                                           gATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGtt                                                               >  1:1041792/1‑90 (MQ=255)
                                                                                                                                           gATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTAt                                                            >  1:922880/1‑93 (MQ=255)
                                                                                                                                           gATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTAt                                                            <  2:922880/93‑1 (MQ=255)
                                                                                                                                           gATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAccc                                             <  2:131643/108‑1 (MQ=255)
                                                                                                                                           gATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAccc                                             >  1:131643/1‑108 (MQ=255)
                                                                                                                                             ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCAc                                                                                      >  2:156744/1‑65 (MQ=255)
                                                                                                                                             ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCAc                                                                                      >  4:469376/1‑65 (MQ=255)
                                                                                                                                             ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCAc                                                                                      <  3:469376/65‑1 (MQ=255)
                                                                                                                                             ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTg                                                                            <  4:950215/75‑1 (MQ=255)
                                                                                                                                             ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGa            >  1:548504/1‑139 (MQ=255)
                                                                                                                                             ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGa            <  2:548504/139‑1 (MQ=255)
                                                                                                                                             ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGa   <  2:756617/148‑1 (MQ=255)
                                                                                                                                              tCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCca                                          >  4:671463/1‑108 (MQ=255)
                                                                                                                                              tCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCca                                          <  3:671463/108‑1 (MQ=255)
                                                                                                                                              tCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGAt  >  3:1026840/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGAT  >  NC_000913/3971849‑3972138

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: