Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,181,198 G→A A35T (GCC→ACC)  yegV → putative sugar kinase YegV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,181,1980GA99.1% 1260.6 / NA 327A35T (GCC→ACC) yegVputative sugar kinase YegV
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (164/160);  minor base C (0/3);  total (164/163)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

GCGGCGCTCGATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGG  >  NC_000913/2181100‑2181295
                                                                                                  |                                                                                                 
gcggcgCTCGATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATAcc                                                                                                 >  1:3927995/1‑101 (MQ=255)
  ggcgCTCGATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCt                                                                                               >  1:3927999/1‑101 (MQ=255)
    cgcTCGATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTcc                                                                                             >  2:3928002/1‑101 (MQ=255)
     gcTCGATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTccc                                                                                            <  2:3927418/101‑1 (MQ=255)
     gcTCGATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTccc                                                                                            >  2:3928006/1‑101 (MQ=255)
     gcTCGATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTccc                                                                                            >  1:3928008/1‑101 (MQ=255)
      cTCGATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTcccc                                                                                           >  1:3928012/1‑101 (MQ=255)
      cTCGATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTcccc                                                                                           <  2:3927570/101‑1 (MQ=255)
         gATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTgg                                                                                        >  2:3928018/1‑101 (MQ=255)
         gATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTgg                                                                                        <  1:3927563/101‑1 (MQ=255)
         gATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTg                                                                                         >  1:3928020/1‑100 (MQ=255)
           ttACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTg                                                                                         >  1:3928026/1‑98 (MQ=255)
           ttACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCg                                                                                      >  2:3928024/1‑101 (MQ=255)
             acacacGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTg                                                                                    >  1:3928029/1‑101 (MQ=255)
               acacGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTggg                                                                                  >  2:3928033/1‑101 (MQ=255)
               acacGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTggg                                                                                  >  1:3928035/1‑101 (MQ=255)
                cacGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGt                                                                                 <  2:3927404/101‑1 (MQ=255)
                 acGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTg                                                                                <  1:3927473/101‑1 (MQ=255)
                 acGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTg                                                                                <  2:3927227/101‑1 (MQ=255)
                   gctgctGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTg                                                                                >  1:3928040/1‑99 (MQ=255)
                     tgctgcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGat                                                                            <  1:3927541/101‑1 (MQ=255)
                     tgctgcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGat                                                                            <  1:3927478/101‑1 (MQ=255)
                     tgctgcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGat                                                                            <  2:3927467/101‑1 (MQ=255)
                      gctgcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGata                                                                           >  2:3928045/1‑101 (MQ=255)
                      gctgcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGata                                                                           >  1:3928046/1‑101 (MQ=255)
                      gctgcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGata                                                                           <  2:3927526/101‑1 (MQ=255)
                      gctgcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGata                                                                           <  2:3927715/101‑1 (MQ=255)
                       ctgcCTGAATTAACCCCGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGatat                                                                          <  2:3927600/101‑1 (MQ=255)
                        tgcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTg                                                                                         <  1:3927581/85‑1 (MQ=255)
                        tgcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATc                                                                         <  1:3927598/101‑1 (MQ=255)
                        tgcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATc                                                                         >  1:3928051/1‑101 (MQ=255)
                         gcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCg                                                                        >  2:3928054/1‑101 (MQ=255)
                         gcCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCg                                                                        <  2:3927761/101‑1 (MQ=255)
                          ccTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGa                                                                       >  2:3928057/1‑101 (MQ=255)
                             gAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACt                                                                    <  1:3927649/101‑1 (MQ=255)
                              aaTTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTg                                                                   >  1:3928062/1‑101 (MQ=255)
                               aTTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGa                                                                  <  2:3927591/101‑1 (MQ=255)
                               aTTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGa                                                                  >  2:3928064/1‑101 (MQ=255)
                                ttAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATCCCCTCCCCTGGCGTGGGTGCGATATCGAACTGaa                                                                 <  1:3926935/101‑1 (MQ=255)
                                 tAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGaaa                                                                <  1:3927381/101‑1 (MQ=255)
                                    ccACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAAcag                                                             <  2:3927453/101‑1 (MQ=255)
                                    ccACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAAcag                                                             <  1:3927530/101‑1 (MQ=255)
                                     cACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAAcagc                                                            >  1:3928069/1‑101 (MQ=255)
                                       cgcgTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATCCCCTCCCCTGGCGTGGGTGCg                                                                              <  1:3927831/81‑1 (MQ=255)
                                       cgcgTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAAcagcag                                                          <  1:3927681/101‑1 (MQ=255)
                                       cgcgTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAAcagcag                                                          >  2:3928073/1‑101 (MQ=255)
                                          gTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGc                                                       <  2:3927362/101‑1 (MQ=255)
                                            cAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGt                                                     >  2:3928078/1‑101 (MQ=255)
                                            cAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGt                                                     >  1:3928077/1‑101 (MQ=255)
                                             aGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGtt                                                    <  2:3927515/101‑1 (MQ=255)
                                              gTCGGTGATGGTCGTCGGCTCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTa                                                   <  2:3927787/101‑1 (MQ=255)
                                               tCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTaa                                                  <  1:3927505/101‑1 (MQ=255)
                                               tCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTaa                                                  <  2:3927490/101‑1 (MQ=255)
                                                 ggTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACg                                                >  2:3928085/1‑101 (MQ=255)
                                                 ggTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACg                                                >  1:3928083/1‑101 (MQ=255)
                                                   tGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTcc                                                                                             >  2:3928091/1‑54 (MQ=255)
                                                   tGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGtt                                              <  2:3927504/101‑1 (MQ=255)
                                                    gATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTg                                             <  2:3927519/101‑1 (MQ=255)
                                                     aTGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTgg                                            >  2:3928094/1‑101 (MQ=255)
                                                      tGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTggc                                           >  1:3928095/1‑101 (MQ=255)
                                                       ggTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCg                                                                              >  2:3928099/1‑65 (MQ=255)
                                                          cgtcgGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCt                                       <  1:3927728/101‑1 (MQ=255)
                                                           gtcgGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTggc                                           >  2:3928106/1‑96 (MQ=255)
                                                           gtcgGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTg                                      >  1:3928104/1‑101 (MQ=255)
                                                             cggcgcggcgGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTgcg                                    <  1:3927323/101‑1 (MQ=255)
                                                                 gcggcgGTCATTGACGTGATCGCCGACGCTTATACCCTccc                                                                                            <  1:3927417/41‑1 (MQ=255)
                                                                 gcggcgGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACt                                <  1:3927617/101‑1 (MQ=255)
                                                                 gcggcgGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACt                                <  2:3927693/101‑1 (MQ=255)
                                                                 gcggcgGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACt                                >  2:3928110/1‑101 (MQ=255)
                                                                 gcggcgGTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACt                                >  1:3928112/1‑101 (MQ=255)
                                                                   ggcggTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGa                              >  2:3928115/1‑101 (MQ=255)
                                                                      ggTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAata                           <  1:3927432/101‑1 (MQ=255)
                                                                      ggTCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAata                           >  1:3928117/1‑101 (MQ=255)
                                                                        tCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATAtt                         >  2:3928122/1‑101 (MQ=255)
                                                                        tCATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATAtt                         <  1:3927263/101‑1 (MQ=255)
                                                                         cATTGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTg                        <  1:3927452/101‑1 (MQ=255)
                                                                           ttGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGcc                      <  1:3927446/101‑1 (MQ=255)
                                                                           ttGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGcc                      >  1:3928125/1‑101 (MQ=255)
                                                                            tGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCa                                  >  2:3928131/1‑88 (MQ=255)
                                                                            tGACGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCg                     >  1:3928128/1‑101 (MQ=255)
                                                                              aCGTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTg                   <  1:3927497/101‑1 (MQ=255)
                                                                                gTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGc                 >  2:3928136/1‑101 (MQ=255)
                                                                                gTGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGc                 <  1:3927743/101‑1 (MQ=255)
                                                                                 tGATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCg                >  1:3928138/1‑101 (MQ=255)
                                                                                  gATCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGt               <  2:3927896/101‑1 (MQ=255)
                                                                                   aTCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTg                                             >  1:3928144/1‑70 (MQ=255)
                                                                                   aTCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGtt              >  2:3928142/1‑101 (MQ=255)
                                                                                   aTCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGtt              <  2:3927550/101‑1 (MQ=255)
                                                                                   aTCGCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGtt              <  1:3927559/101‑1 (MQ=255)
                                                                                      gCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTaaa           <  2:3927661/101‑1 (MQ=255)
                                                                                      gCCGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTaaa           >  2:3928149/1‑101 (MQ=255)
                                                                                        cGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgc         <  1:3927625/101‑1 (MQ=255)
                                                                                        cGACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgc         <  2:3927639/101‑1 (MQ=255)
                                                                                         gACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgcg        >  1:3928154/1‑101 (MQ=255)
                                                                                         gACGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAgcg        >  2:3928153/1‑101 (MQ=255)
                                                                                           cGCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGcc      >  2:3928158/1‑101 (MQ=255)
                                                                                            gCTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGc                       <  2:3927425/83‑1 (MQ=255)
                                                                                             cTTATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTc    <  2:3927341/101‑1 (MQ=255)
                                                                                              ttATACCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCg   >  1:3928162/1‑101 (MQ=255)
                                                                                               tataCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCgg  <  2:3927620/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
GCGGCGCTCGATTACACACGCTGCTGCCTGAATTAACCACGCGTCAGTCGGTGATGGTCGTCGGCGCGGCGGTCATTGACGTGATCGCCGACGCTTATGCCCTCCCCTGGCGTGGGTGCGATATCGAACTGAAACAGCAGAGCGTTAACGTTGGCGGCTGCGCACTGAATATTGCCGTGGCGTTAAAGCGCCTCGG  >  NC_000913/2181100‑2181295

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