Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 102,844 A→G Q204Q (CAA→CAG ddlB → D‑alanine‑D‑alanine ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913102,8440AG96.5% 527.4 / 3.8 174Q204Q (CAA→CAGddlBD‑alanine‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base A (1/5);  new base G (86/82);  total (87/87)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 180 total aligned reads displayed.

TGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAACCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTCACAAGA  >  NC_000913/102754‑102933
                                                                                          |                                                                                         
tGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGcc                                                                                         <  1:3241966/93‑1 (MQ=255)
  aagaagTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGt                                                                                       >  2:3150919/1‑93 (MQ=255)
   agaagTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTc                                                                                      >  2:1927281/1‑93 (MQ=255)
   agaagTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTc                                                                                      >  1:3477827/1‑93 (MQ=255)
    gaagTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACTGTCAATACGTATTCAGCCGTcc                                                                                     >  2:3249079/1‑93 (MQ=255)
         attgattgAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAAc                                                                                >  2:415208/1‑93 (MQ=255)
          ttgattgaAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAAcc                                                                               <  1:1361825/93‑1 (MQ=255)
          ttgattgaAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAAcc                                                                               >  2:3693458/1‑93 (MQ=255)
          ttgattgaAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGCCCGGAAcc                                                                               >  1:4124358/1‑93 (MQ=255)
             attgaAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTc                                                                            <  2:1532219/93‑1 (MQ=255)
              ttgaAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCt                                                                           <  1:1200707/93‑1 (MQ=255)
               tgaAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTa                                                                          <  2:3197737/93‑1 (MQ=255)
                gaAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAACCGTCCGGAACCTTCTAt                                                                         <  1:3094903/93‑1 (MQ=255)
                  aaaaTGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGa                                                                       >  1:1047030/1‑93 (MQ=255)
                   aaaTGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGAt                                                                      <  1:2437343/93‑1 (MQ=255)
                    aaTGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGAtt                                                                     >  1:3138936/1‑93 (MQ=255)
                    aaTGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGAtt                                                                     >  1:3005357/1‑93 (MQ=255)
                     aTGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTa                                                                    >  2:3967848/1‑93 (MQ=255)
                     aTGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTa                                                                    >  1:986508/1‑93 (MQ=255)
                        gCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGa                                                                 <  2:2764745/93‑1 (MQ=255)
                        gCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGa                                                                 <  1:2209615/93‑1 (MQ=255)
                        gCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGa                                                                 <  1:2983756/93‑1 (MQ=255)
                        gCTAAGGGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAACCGTCCGGAACCTTCTATGATTAtca                                                                 <  2:3292252/93‑3 (MQ=255)
                         cTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGCACCTTCTATGATTATGAg                                                                >  1:4405814/1‑93 (MQ=255)
                         cTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAg                                                                <  1:1257016/93‑1 (MQ=255)
                         cTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAg                                                                <  2:844165/93‑1 (MQ=255)
                          tAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAgg                                                               <  2:4138358/93‑1 (MQ=255)
                          tAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAgg                                                               >  1:1469470/1‑93 (MQ=255)
                           aaGTGGGCCGGAGTTCACGGTTGCGATCCTCGGTGAAGAAATTTTCCCGTCAATACGTATTCAACCGTCCGGAACCTTCTATGATTATGAGGc                                                              <  1:3998774/93‑1 (MQ=255)
                           aaGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGc                                                              >  1:4477940/1‑93 (MQ=255)
                             gTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGa                                                            <  2:1430264/93‑1 (MQ=255)
                              tcggCCGGAGTTCACGGATGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAACCGTCCGGGACCTTCCATGATTATGAGGCGaa                                                           >  2:2203361/3‑93 (MQ=255)
                              ggggCCGGAGTTCACGGTTGCTATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGaa                                                           <  2:3461865/92‑1 (MQ=255)
                                ggCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGt                                                         <  2:544393/93‑1 (MQ=255)
                                   cGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTAtt                                                      >  2:1280251/1‑92 (MQ=255)
                                    ggAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTAtct                                                     >  1:167401/1‑93 (MQ=255)
                                    ggAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTAtct                                                     <  1:680034/93‑1 (MQ=255)
                                     gAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTAtctc                                                    <  2:3393391/93‑1 (MQ=255)
                                     gAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTAtctc                                                    >  1:3301231/1‑93 (MQ=255)
                                     gAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTAtctc                                                    <  1:2562348/93‑1 (MQ=255)
                                      aGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGGACCTTCTATGATTATGAGGCGAAGTAtctct                                                   >  2:4712174/1‑93 (MQ=255)
                                      aGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTAtctct                                                   <  1:254318/93‑1 (MQ=255)
                                      aGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTAtctct                                                   <  2:209998/93‑1 (MQ=255)
                                       gTTCACGGTCGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTAtctctc                                                  >  2:2759835/1‑93 (MQ=255)
                                        ttCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATctctct                                                 <  1:3739870/93‑1 (MQ=255)
                                        ttCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATctctct                                                 >  2:4432159/1‑93 (MQ=255)
                                        ttCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATctctct                                                 <  2:4676378/93‑1 (MQ=255)
                                          cACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCtga                                               <  1:1892969/93‑1 (MQ=255)
                                             cggtGCGATAATCTCTGAACAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATATCTCcgatga                                            <  2:3761266/90‑1 (MQ=255)
                                              gTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATgag                                           >  2:22728/1‑93 (MQ=255)
                                              gTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATgag                                           >  2:655416/1‑93 (MQ=255)
                                                tGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGACGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGac                                         <  1:1844053/93‑1 (MQ=255)
                                                  cGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGacac                                       <  2:3639418/93‑1 (MQ=255)
                                                  cGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGacac                                       <  2:4375995/93‑1 (MQ=255)
                                                   gATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGAcaca                                      >  1:3776785/1‑93 (MQ=255)
                                                     tACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCGGTCCGGAACCTTCTATGATTATGAGGCGAAGCATCTCTCTGATGAGACACAGt                                    >  1:1611854/1‑93 (MQ=255)
                                                      aCTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTa                                   >  2:4069702/1‑93 (MQ=255)
                                                      aCTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTa                                   <  2:3480174/93‑1 (MQ=255)
                                                       cTCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTAt                                  <  2:2949322/93‑1 (MQ=255)
                                                        tCGGTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTAtt                                 >  1:1960070/1‑93 (MQ=255)
                                                        acGGTTAAGAAAGTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTAtt                                 <  2:905925/92‑1 (MQ=255)
                                                          ggTGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTc                               <  2:485024/93‑1 (MQ=255)
                                                            tGAAGAAATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTGTCTGATGAGACACAGTATTTCTg                             >  2:3388419/1‑93 (MQ=255)
                                                             gaagaaATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGc                            >  2:1796807/1‑93 (MQ=255)
                                                             gaagaaATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGc                            <  1:4740524/93‑1 (MQ=255)
                                                             gaagaaATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGc                            >  2:2402269/1‑93 (MQ=255)
                                                             gaagaaATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGc                            <  1:2142685/93‑1 (MQ=255)
                                                             gaagaaATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGc                            >  2:2685801/1‑93 (MQ=255)
                                                               agaaATTTTATCGTCAATACATATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGccc                          <  2:1349097/93‑1 (MQ=255)
                                                               agaaATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGccc                          >  1:4604674/1‑93 (MQ=255)
                                                               agaaATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGccc                          >  2:3750978/1‑93 (MQ=255)
                                                                gaaATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGcccc                         >  2:250759/1‑93 (MQ=255)
                                                                gaaATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGcccc                         <  2:4308196/93‑1 (MQ=255)
                                                                gaaATTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGcccc                         <  2:3726197/93‑1 (MQ=255)
                                                                   aTTTTACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTATCTGATGAGACACAGTATTTCTGCCCAGCa                      >  1:3770021/1‑93 (MQ=255)
                                                                    ttttACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAg                     <  1:2378239/93‑1 (MQ=255)
                                                                    ttttACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAg                     >  1:2923928/1‑93 (MQ=255)
                                                                    ttttACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAg                     >  1:2946541/1‑93 (MQ=255)
                                                                      ttACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGt                   >  1:2531796/1‑93 (MQ=255)
                                                                      ttACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGt                   <  1:3125995/93‑1 (MQ=255)
                                                                       tACCGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTc                  <  2:1998757/93‑1 (MQ=255)
                                                                        tccGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCt                 <  1:4142444/92‑1 (MQ=255)
                                                                         ccGTCAATACGTATTCAGCCGTACGGTACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCCGGCCTg                >  2:2895527/1‑93 (MQ=255)
                                                                          cGTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTgg               >  1:1816350/1‑93 (MQ=255)
                                                                           gTCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGa              >  1:1700883/1‑93 (MQ=255)
                                                                            tCAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGaa             <  2:1148193/93‑1 (MQ=255)
                                                                             cAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAg            >  2:968181/1‑93 (MQ=255)
                                                                             cAATACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAg            <  1:2223385/93‑1 (MQ=255)
                                                                                tACGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGt         >  1:2702873/1‑93 (MQ=255)
                                                                                 aCGTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTc        <  1:1429668/93‑1 (MQ=255)
                                                                                   gTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTcac      <  2:176978/93‑1 (MQ=255)
                                                                                   gTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTcac      >  2:4237098/1‑93 (MQ=255)
                                                                                   gTATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTcac      >  2:2488627/1‑93 (MQ=255)
                                                                                    tATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTcaca     <  1:1286100/93‑1 (MQ=255)
                                                                                    tATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTcaca     >  2:2852206/1‑93 (MQ=255)
                                                                                    tATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTcaca     >  1:4047523/1‑93 (MQ=255)
                                                                                    tATTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTcaca     <  2:2023761/93‑1 (MQ=255)
                                                                                     aTTCAGCCGTCCTGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTCACaa    <  2:1880600/93‑1 (MQ=255)
                                                                                     aTTCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCCGGGCTGGAAGCGCCACaa    >  1:3865447/1‑93 (MQ=255)
                                                                                       tCAGCCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTCACAaga  <  1:714723/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTCAATACGTATTCAACCGTCCGGAACCTTCTATGATTATGAGGCGAAGTATCTCTCTGATGAGACACAGTATTTCTGCCCCGCAGGTCTGGAAGCGTCACAAGA  >  NC_000913/102754‑102933

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