Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 281,596 T→C L39P (CTG→CCG)  yagE → putative lyase/synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,5960TC87.4% 682.0 / 81.4 253L39P (CTG→CCG) yagEputative lyase/synthase
Reads supporting (aligned to +/- strand):  ref base T (19/13);  new base C (97/124);  total (116/137)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

Only 100 of 266 total aligned reads displayed.

CCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAG  >  NC_000913/281505‑281688
                                                                                           |                                                                                            
 cGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGa                                                                                            <  2:24128/93‑1 (MQ=255)
 cGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGa                                                                                            <  1:1125563/93‑1 (MQ=255)
 cGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGa                                                                                            <  1:5651193/93‑1 (MQ=255)
  gCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGAt                                                                                           >  1:1534533/1‑93 (MQ=255)
   cAGTCCGCGTTTTGCACGGGAATCATTCCCCCTGCCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATc                                                                                          <  1:86823/93‑1 (MQ=255)
   cAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATc                                                                                          <  2:2770164/93‑1 (MQ=255)
    aGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATcg                                                                                         >  1:185549/1‑93 (MQ=255)
      tCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATcgac                                                                                       >  2:5956469/1‑93 (MQ=255)
         gcgTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGAt                                                                                    >  1:953223/1‑93 (MQ=255)
          cgTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATc                                                                                   <  1:5758229/93‑1 (MQ=255)
          cgTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATc                                                                                   <  2:2967276/93‑1 (MQ=255)
           gttgttCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCt                                                                                  >  2:4277791/1‑93 (MQ=255)
           gttgttCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCt                                                                                  >  2:1142321/1‑93 (MQ=255)
            ttgttCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTg                                                                                 >  1:6003273/1‑93 (MQ=255)
            ttgttCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTg                                                                                 >  2:3901434/1‑93 (MQ=255)
            ttgttCACGGGAATCATTCCCCCTGCCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTg                                                                                 <  1:6496248/93‑1 (MQ=255)
              gttCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGAt                                                                               <  2:1475368/93‑1 (MQ=255)
              gttCACGGGAATCATCCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGAt                                                                               <  1:172761/93‑1 (MQ=255)
                tCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCa                                                                             >  2:3182843/1‑93 (MQ=255)
                 cACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCaa                                                                            >  2:3286280/1‑93 (MQ=255)
                    gggAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGc                                                                         <  2:5685751/93‑1 (MQ=255)
                     ggAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCa                                                                        <  2:1093215/93‑1 (MQ=255)
                     ggAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCa                                                                        >  1:6406573/1‑93 (MQ=255)
                     gcaaTCATTCCCCCTGTCCCCCCCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCa                                                                        <  1:1211242/91‑1 (MQ=255)
                      gAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAg                                                                       <  1:3476713/93‑1 (MQ=255)
                       aaTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAgg                                                                      <  1:2647608/93‑1 (MQ=255)
                       aaTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAgg                                                                      <  2:2026031/93‑1 (MQ=255)
                       aaTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAgg                                                                      <  1:6610915/93‑1 (MQ=255)
                        aTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGc                                                                     <  1:2706052/93‑1 (MQ=255)
                         tCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCg                                                                    >  1:3337505/1‑93 (MQ=255)
                         tCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCg                                                                    <  1:5473289/93‑1 (MQ=255)
                          cATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCTATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGt                                                                   >  2:951264/1‑93 (MQ=255)
                           aTTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGTTCTGATCAAAGCAGGCGtt                                                                  <  1:6485268/93‑1 (MQ=255)
                           aTTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGtt                                                                  >  2:3457489/1‑93 (MQ=255)
                            ttCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGCTCGACGATCTGATCAAAGCAGGCGATg                                                                 >  2:2839876/1‑93 (MQ=255)
                              cccccTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGAc                                                               >  1:4184710/1‑93 (MQ=255)
                               ccccTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACg                                                              <  2:1342363/93‑1 (MQ=255)
                                cccTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACgg                                                             >  1:3253147/1‑93 (MQ=255)
                                 ccTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGc                                                            <  2:5572201/93‑1 (MQ=255)
                                 ccTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGc                                                            <  2:1790330/93‑1 (MQ=255)
                                 ccTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGc                                                            >  2:4436215/1‑93 (MQ=255)
                                  cTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGcc                                                           >  1:2301307/1‑93 (MQ=255)
                                    gTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTg                                                         >  2:6441162/1‑93 (MQ=255)
                                     tctcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGt                                                        >  2:5019479/1‑93 (MQ=255)
                                       tcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGttt                                                      >  2:1053291/1‑92 (MQ=255)
                                       tcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGttc                                                      <  2:773523/93‑1 (MQ=255)
                                        ccaccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGttct                                                     >  1:4397988/1‑93 (MQ=255)
                                        ccaccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGttct                                                     >  1:3422604/1‑93 (MQ=255)
                                         caccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGttctt                                                    <  1:2821779/93‑1 (MQ=255)
                                          accaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGttcttc                                                   <  2:2261152/93‑1 (MQ=255)
                                           ccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTcc                                                  >  1:2466959/1‑93 (MQ=255)
                                           ccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTcc                                                  >  1:1772061/1‑93 (MQ=255)
                                            caTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCt                                                 <  1:2940972/93‑1 (MQ=255)
                                            caTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCt                                                 >  2:143755/1‑93 (MQ=255)
                                            caTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCt                                                 >  1:4745516/1‑93 (MQ=255)
                                             aTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTg                                                >  2:1908950/1‑93 (MQ=255)
                                              tttttACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTgg                                               >  2:4621878/1‑93 (MQ=255)
                                              tttttACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTgg                                               <  1:2593394/93‑1 (MQ=255)
                                                tttACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACTGCCTGTTCTTCCTGGGc                                             <  2:1380441/93‑1 (MQ=255)
                                                 ttACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCa                                            >  1:6518192/1‑93 (MQ=255)
                                                  tACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAg                                           >  2:5146497/1‑93 (MQ=255)
                                                  tACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGCTCTTCCTGGGCAg                                           >  2:5253821/1‑93 (MQ=255)
                                                      gccgACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGt                                       <  1:4818927/93‑1 (MQ=255)
                                                        cgACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTgg                                     <  1:4271325/93‑1 (MQ=255)
                                                             gCCAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGt                                >  2:1160365/1‑93 (MQ=255)
                                                              ccAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGtt                               >  2:6298825/1‑93 (MQ=255)
                                                              ccAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGtt                               >  2:5498952/1‑93 (MQ=255)
                                                               cAGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTtc                              <  2:2917156/93‑1 (MQ=255)
                                                                aGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTtct                             >  1:4249409/1‑93 (MQ=255)
                                                                aGCTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTtct                             >  2:3327739/1‑93 (MQ=255)
                                                                  cTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTcc                           <  2:4216951/93‑1 (MQ=255)
                                                                  cTCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGACCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTcc                           <  2:4741171/93‑1 (MQ=255)
                                                                   tCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTccc                          <  1:5362130/93‑1 (MQ=255)
                                                                   tCGATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTccc                          <  1:3932111/93‑1 (MQ=255)
                                                                     gATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAg                        >  2:2453037/1‑93 (MQ=255)
                                                                     gATAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAg                        >  1:2119362/1‑93 (MQ=255)
                                                                      aTAAGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGc                       <  1:5307834/93‑1 (MQ=255)
                                                                       tAAGCCGGACACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCt                      <  2:1943580/93‑1 (MQ=255)
                                                                         aGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCg                    <  1:3830968/93‑1 (MQ=255)
                                                                         aGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCg                    <  2:5831531/93‑1 (MQ=255)
                                                                         aGCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGCTCTCCCAGCTCg                    >  1:4776130/1‑93 (MQ=255)
                                                                          gCCGGGCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGCTCTCCCAGCTCgg                   >  2:1196064/1‑93 (MQ=255)
                                                                             gggCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgcgc                <  2:6551546/93‑1 (MQ=255)
                                                                             gggCACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgcgc                <  2:2521441/93‑1 (MQ=255)
                                                                                cACCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGa             <  2:663294/93‑1 (MQ=255)
                                                                                 cccgccgCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGaa            <  1:3592922/92‑1 (MQ=255)
                                                                                 aCCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGaa            <  1:2274282/93‑1 (MQ=255)
                                                                                 aCCGCCGCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGaa            >  1:1454617/1‑93 (MQ=255)
                                                                                  ccgccgCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAag           >  1:2023304/1‑93 (MQ=255)
                                                                                  ccgccgCGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAag           <  1:1283703/93‑1 (MQ=255)
                                                                                   cgccgcgccgATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAaga          <  1:6328720/93‑1 (MQ=255)
                                                                                   cgccgcgccgATCGACGATCTGATCAAAGCAGGCATTGACGTCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAaga          <  2:4062981/93‑1 (MQ=255)
                                                                                   cgccgcGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAaga          <  1:3006484/93‑1 (MQ=255)
                                                                                   ccccgcGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAaga          <  2:3527823/91‑1 (MQ=255)
                                                                                    gccgcGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAagag         <  2:5860324/93‑1 (MQ=255)
                                                                                    gccgcGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAagag         <  2:487373/93‑1 (MQ=255)
                                                                                    gccgcGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCGTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAagag         >  1:6146789/1‑93 (MQ=255)
                                                                                     ccgcGCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGc        <  1:2142771/93‑1 (MQ=255)
                                                                                      cgcgCCGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCg       <  2:1685503/93‑1 (MQ=255)
                                                                                       gcgcCGATCGACGATCTGATCAAAGCAGGCGTTGACGCCCTGTTCTTCCTGGGCAGCTGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGt      <  1:4671041/93‑1 (MQ=255)
                                                                                           |                                                                                            
CCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAG  >  NC_000913/281505‑281688

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: