Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,992,862 A→G Q164R (CAG→CGG)  dapF → diaminopimelate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,992,8620AG99.7% 1111.0 / ‑5.4 356Q164R (CAG→CGG) dapFdiaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (171/184);  total (172/184)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 375 total aligned reads displayed.

ACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTT  >  NC_000913/3992772‑3992951
                                                                                          |                                                                                         
aCAAAGCGGAAAAGATCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCggg                                                                                         <  1:1004681/93‑1 (MQ=255)
aCAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCggg                                                                                         <  1:6787583/93‑1 (MQ=255)
aCAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCggg                                                                                         >  1:3155878/1‑93 (MQ=255)
  aaaGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGtg                                                                                       >  1:4629754/1‑92 (MQ=255)
   aaGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCg                                                                                      >  2:2802769/1‑93 (MQ=255)
    aGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCga                                                                                     <  2:7694557/93‑1 (MQ=255)
    aGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCga                                                                                     <  2:2156787/93‑1 (MQ=255)
     gCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgat                                                                                    <  1:5709078/93‑1 (MQ=255)
     gCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgat                                                                                    <  1:1467449/93‑1 (MQ=255)
     gCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgat                                                                                    >  1:6871438/1‑93 (MQ=255)
       ggAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgatga                                                                                  <  1:720079/93‑1 (MQ=255)
       ggAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgatga                                                                                  >  2:3176167/1‑93 (MQ=255)
        gAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgatgat                                                                                 <  1:5472067/93‑1 (MQ=255)
        gAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgatgat                                                                                 >  2:1520322/1‑93 (MQ=255)
         aaaaGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGatgatg                                                                                >  1:3967491/1‑93 (MQ=255)
         aaaaGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGatgatg                                                                                >  1:1293487/1‑93 (MQ=255)
          aaaGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGt                                                                               >  2:3396412/1‑93 (MQ=255)
          aaaGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGt                                                                               >  2:6353763/1‑93 (MQ=255)
           aaGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTc                                                                              >  1:1164055/1‑93 (MQ=255)
           aaGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTc                                                                              <  2:7487101/93‑1 (MQ=255)
             gACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGa                                                                            <  2:1237029/93‑1 (MQ=255)
             gACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGa                                                                            <  1:5077051/93‑1 (MQ=255)
              aCCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGAt                                                                           <  1:6260950/93‑1 (MQ=255)
               ccTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATa                                                                          <  2:5554608/93‑1 (MQ=255)
               ccTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATa                                                                          >  2:816474/1‑93 (MQ=255)
               ccTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGAAATTCCGCATTGCGTGATTCGGGTCGATGATGTCGATa                                                                          <  2:7051914/93‑1 (MQ=255)
                cTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATAc                                                                         <  2:6406486/93‑1 (MQ=255)
                 tataTTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATAcc                                                                        >  2:3812784/1‑93 (MQ=255)
                  atatTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACcg                                                                       >  2:6165364/1‑93 (MQ=255)
                   tattatGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACcgc                                                                      <  2:6003323/93‑1 (MQ=255)
                   tattatGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACcgc                                                                      >  2:170619/1‑93 (MQ=255)
                   tattatGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACcgc                                                                      >  2:252363/1‑93 (MQ=255)
                   tattatGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACcgc                                                                      <  2:3244357/93‑1 (MQ=255)
                    attatGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACcgcg                                                                     >  2:4137918/1‑93 (MQ=255)
                      tatGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCgcggc                                                                   >  1:5312929/1‑93 (MQ=255)
                        tGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGcggcgg                                                                 >  1:6645518/1‑93 (MQ=255)
                          ggcgcCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTa                                                               <  1:5408658/92‑1 (MQ=255)
                          cgcgcCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTa                                                               <  1:7532114/93‑1 (MQ=255)
                           gcgcCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAg                                                              <  1:2005290/93‑1 (MQ=255)
                             gccgccGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGaa                                                            <  1:3309320/93‑1 (MQ=255)
                              ccgccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGaag                                                           <  1:3756758/93‑2 (MQ=255)
                               cgccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGaagc                                                          >  2:3068698/1‑91 (MQ=255)
                                 ccgAGCGGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAgcgc                                                        >  1:963400/1‑93 (MQ=255)
                                 ccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGTCGGTAGAAgcgc                                                        <  2:5756917/93‑1 (MQ=255)
                                 ccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAgcgc                                                        >  1:1058033/1‑93 (MQ=255)
                                  cgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCt                                                       <  1:4850092/93‑1 (MQ=255)
                                   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                      >  1:6609119/1‑93 (MQ=255)
                                   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                      <  2:6959193/93‑1 (MQ=255)
                                   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                      >  1:5515403/1‑93 (MQ=255)
                                   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                      >  1:447179/1‑93 (MQ=255)
                                    aGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTg                                                     >  1:745954/1‑93 (MQ=255)
                                    aGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTg                                                     >  1:6064145/1‑93 (MQ=255)
                                     gCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTgg                                                    >  1:1575045/1‑93 (MQ=255)
                                       aGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTc                                                  >  1:2902241/1‑93 (MQ=255)
                                        gACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTcc                                                 >  2:2862279/1‑93 (MQ=255)
                                        gACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTcc                                                 <  2:7305956/93‑1 (MQ=255)
                                         aCAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCt                                                >  1:5426453/1‑93 (MQ=255)
                                          cAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTg                                               <  1:2770159/93‑1 (MQ=255)
                                           aaTCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGt                                              <  1:2514516/93‑1 (MQ=255)
                                           aaTCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGt                                              >  2:4025832/1‑93 (MQ=255)
                                               ttATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTg                                          >  1:7811586/1‑93 (MQ=255)
                                               ttATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTg                                          <  2:4949147/93‑1 (MQ=255)
                                                   gcggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGaaa                                      <  2:7848529/93‑1 (MQ=255)
                                                     ggcgTGGTTTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGc                                    <  1:377578/93‑1 (MQ=255)
                                                       cgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCa                                  <  2:2309786/93‑1 (MQ=255)
                                                       cgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCa                                  <  1:6973726/93‑1 (MQ=255)
                                                        gtggtgTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCAc                                 <  1:4343463/93‑1 (MQ=255)
                                                         tggtgTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACg                                >  1:409769/1‑93 (MQ=255)
                                                          ggtgTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGa                               <  2:5116973/93‑1 (MQ=255)
                                                           gtgtGGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAg                              <  1:5220493/93‑1 (MQ=255)
                                                            tgtCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGc                             >  1:365645/1‑93 (MQ=255)
                                                            tgtCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGc                             <  1:2390792/93‑1 (MQ=255)
                                                            tgtCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGc                             <  2:7550805/93‑1 (MQ=255)
                                                              tCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGt                           <  1:492226/93‑1 (MQ=255)
                                                              tCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGt                           <  2:4681940/93‑1 (MQ=255)
                                                                gATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGttt                         >  1:4216323/1‑93 (MQ=255)
                                                                gATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGttt                         >  1:7865167/1‑93 (MQ=255)
                                                                gATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGttt                         <  1:7599240/93‑1 (MQ=255)
                                                                gATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGttt                         >  2:4784542/1‑93 (MQ=255)
                                                                 aTGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGtttt                        <  2:4423113/93‑1 (MQ=255)
                                                                    ggAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCg                     <  1:560894/93‑1 (MQ=255)
                                                                     gAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTGCgg                    <  2:5316674/93‑1 (MQ=255)
                                                                     gAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCgg                    <  2:657617/93‑1 (MQ=255)
                                                                      aaaTCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGa                   <  1:3035339/93‑1 (MQ=255)
                                                                       aaTCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAg                  <  2:3651087/93‑1 (MQ=255)
                                                                       aaTCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTATGGAAAGCCACGAGCGTTTTCCGGAg                  >  2:4448307/1‑93 (MQ=255)
                                                                        aTCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgc                 >  1:6402291/1‑93 (MQ=255)
                                                                        aTCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgc                 >  2:5058944/1‑93 (MQ=255)
                                                                          ccGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgc               <  2:3508494/93‑1 (MQ=255)
                                                                           cGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgcg              >  2:2050972/1‑93 (MQ=255)
                                                                             cATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGcc            >  2:2385420/1‑93 (MQ=255)
                                                                               ttGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCaa          <  2:1379144/93‑1 (MQ=255)
                                                                                  cGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAatat       <  1:2143029/93‑1 (MQ=255)
                                                                                  cGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAatat       <  1:4551695/93‑1 (MQ=255)
                                                                                   gTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAatatg      >  1:4729572/1‑92 (MQ=255)
                                                                                   gTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATc      >  2:6865449/1‑93 (MQ=255)
                                                                                   gTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATc      <  2:2778699/93‑1 (MQ=255)
                                                                                     gATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCgg    >  1:5574030/1‑93 (MQ=255)
                                                                                       ttCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGtt  <  1:5951781/93‑1 (MQ=255)
                                                                                       ttCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGtt  <  1:1898571/93‑1 (MQ=255)
                                                                                          |                                                                                         
ACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTT  >  NC_000913/3992772‑3992951

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: ‑